matlab solve 参数,用solve如何求解带有一个变参数的方程组

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