
clf

experiment=1;

t=[0:0.01:2*pi 0];

m1=20; m2=4;

xt=cos(t)/m1;
yt=sin(t)/m2;
xt1=0.75+1*cos(t)/m1;
yt1=1*sin(t)/m2;




%axis equal
axis off

hold on

switch experiment
    case 1
        
fill(xt,yt,'m-');
fill(xt1,yt1,'m-');
plot(xt1,yt1,'k-');
plot(xt,yt,'k-');

    case 2
     
        fill(xt,yt,'m-');
fill(xt1,yt1,'m-');
plot(xt1,yt1,'k-');
plot(xt,yt,'k-');
        
end
hold off

