rp=128;
ct=0:1:90;
fan10=11.25;
lamt=0;
ctf=135;
T1=ct/ctf;
if T1<=0.125
fani=1/(4+pi)*(pi*T1-0.25*sin(4*pi*T1));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
if T1>=0.875
fani=1/(4+PI)*(4+pi*T1-0.25*sin(4*pi*T1));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
if T1>0.125&&T1<0.875
fani=1/(4+pi)*(2+pi*T1-9/4*sin((PI+4*pi*T1)/3));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
figure
plot(x,y)
ct=0:1:90;
fan10=11.25;
lamt=0;
ctf=135;
T1=ct/ctf;
if T1<=0.125
fani=1/(4+pi)*(pi*T1-0.25*sin(4*pi*T1));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
if T1>=0.875
fani=1/(4+PI)*(4+pi*T1-0.25*sin(4*pi*T1));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
if T1>0.125&&T1<0.875
fani=1/(4+pi)*(2+pi*T1-9/4*sin((PI+4*pi*T1)/3));
fan=fani+fan10;
x=rp*sin(ct+fan+fan10-lamt)-C*sin(ct-fan10-lamt);
y=-rp*cos(ct+fan-fan10-lamt)+C*cos(ct-fan10-lamt);
end
figure
plot(x,y)