clearclcticn=1000;h=linspace(0,0.4,n);i=25;R=0.16;r=0.019;r0=0.005;d=0.175;sd0=h/r0;sd1=sd0/i;ac=0.115;bc=0.065;l=0.1166;cc=sd0;fun=@(aa)((bc)^2)*((1-cos(aa))^2)+(ac+r*cos(cc)-bc*sin(aa))^2+(r*sin(cc))^2-l^2aa=solve(fun,0.0);