Need help making Delaunay's Disk in a figure.

Essentially what the title says. If you're not sure what the disk is, I've uploaded a picture of it. I tried mapping it out this way:

figure

hold on

x=n*[cos(pi/2),cos(pi),cos(3*pi/2),cos(2*pi)];

y=n*[sin(pi/2),sin(pi),sin(3*pi/2),sin(2*pi)];

plot(x,y)

fill(x,y, [ ])

But realized that this would only map a square with those points as vertices. I'm not entirely sure what steps to take next, though. I'm sorry ahead of time if anything about the code is wrong, I'm in my first semester of learning matlab.

Answer 1

[Akshit Agarwal]

Hi Christophre,

Thats a interesting problem statement you got there. You can refer to these queries that have already been answered:

How to plot a disc with a given center c and radius r

Drawing a segment of a circle

How to draw a filled semi circle?(like the one in the image) Thanks,

Solutions provided in above links should be enough for you to write the code for your problem statement.

[Christophre Dennen]

Hey Akshit,

This is the script I ended up using to get it to work, for anyone who might find it interesting or helpful!

clc
figure
hold on
t=linspace(0,0.5*pi,70);
 
x=[0 cos(t) 0];
y=[0,sin(t),0];
 
fill(x,y,[rand,rand,rand])
fill(x/1.15,y/1.15,[rand,rand,rand])
fill(x/1.35,y/1.35,[rand,rand,rand])
fill(x/1.65,y/1.65,[rand,rand,rand])
fill(x/2.05,y/2.05,[rand,rand,rand])
fill(x/2.75,y/2.75,[rand,rand,rand])
fill(x/4.35,y/4.35,[rand,rand,rand])
 
fill(-x,y,[rand,rand,rand])
fill(-x/1.15, y/1.15,[rand,rand,rand])
fill(-x/1.35,y/1.35,[rand,rand,rand])
fill(-x/1.65,y/1.65,[rand,rand,rand])
fill(-x/2.05,y/2.05,[rand,rand,rand])
fill(-x/2.75,y/2.75,[rand,rand,rand])
fill(-x/4.35,y/4.35,[rand,rand,rand])
 
fill(x,-y,[rand,rand,rand])
fill(x/1.15,-y/1.15,[rand,rand,rand])
fill(x/1.35,-y/1.35,[rand,rand,rand])
fill(x/1.65,-y/1.65,[rand,rand,rand])
fill(x/2.05,-y/2.05,[rand,rand,rand])
fill(x/2.75,-y/2.75,[rand,rand,rand])
fill(x/4.35,-y/4.35,[rand,rand,rand])
 
fill(-x,-y,[rand,rand,rand])
fill(-x/1.15,-y/1.15,[rand,rand,rand])
fill(-x/1.35,-y/1.35,[rand,rand,rand])
fill(-x/1.65,-y/1.65,[rand,rand,rand])
fill(-x/2.05,-y/2.05,[rand,rand,rand])
fill(-x/2.75,-y/2.75,[rand,rand,rand])
fill(-x/4.35,-y/4.35,[rand,rand,rand])
 
saveas(gcf, 'disk','jpg')