function question_24
clear;clc
b=0.8 ;
a=pi*5/180;
R0=30:1:50;
l=82.0;
T=300.0;
L=250.0;
C=100.0;
K=30.0;
c=200.0;
e=10.0;
f=5.0;
g=5.0;
i=10.0;
j=10.0;
k=150.0;
x0=[16 2 -1 0.1 100 60 0.6];
f=@(R) @(x) [b*x(5)*sin(x(7))-x(5)*cos(x(7))+b*x(6)-c*sin(a);
b*x(5)*cos(x(7))+x(5)*sin(x(7))+x(6)-c*cos(a);
sin(x(7))-(R-x(1))/R;
b*x(5)*R+x(6)*(b*R-sqrt((2*L-l*(sin(x(3))+sin(x(2))))^2+l^2*(cos(x(3))-cos(x(2)))^2))+e*cos(a)*sqrt((2*L-l*(sin(x(3))+sin(x(2))))^2+l^2*(cos(x(3))-cos(x(2)))^2)-0.5*f*l*sin(a+x(4)+x(2))-i*(l*sin(a+x(4)+x(2))+0.5*T*sin(a+x(4)))+g*(sqrt((2*L-l*(sin(x(3))+sin(x(2))))^2+l^2*(cos(x(3))-cos(x(2)))^2)*cos(a)+0.5*l*sin(x(3)-a-x(4)))+j*(sqrt((2*L-l*(sin(x(3))+sin(x(2))))^2+l^2*(cos(x(3))-cos(x(2)))^2)*cos(a)+l*sin(x(3)-a-x(4))-0.5*T*sin(a+x(4)))+k*((L-C*tan(x(4)+a))*cos(x(4)+a)-l*sin(a+x(4)+x(2))-T*sin(a+x(4)));;
K*180*x(2)/pi-b*x(5)*(l*sin(x(7)+x(2)+x(4))-R)+x(5)*l*cos(x(7)+x(2)+x(4));
-K*180*x(3)/pi-b*x(6)*(l*cos(x(3)-x(4))-R)+x(6)*l*sin(x(3)-x(4));
tan(x(4))-l*(cos(x(3))-cos(x(2)))/(2*L-l*(sin(x(3))+sin(x(2))));];
res=arrayfun(@(R) fsolve(f(R),x0),R0,'UniformOutput',false);
res{1}
res{2}
res{end}
end