matlab fsolve传递变量,使用fsolve函数，程序提示我变量不够，如何进行参数的传递...

t0=[0;7;1;2;0];

a=0:0.25:5;

p1=zeros(1,20);

p2=zeros(1,20);

p3=zeros(1,20);

p4=zeros(1,20);

p5=zeros(1,20);

p6=zeros(1,20);

for i=1:length(a)

sol=fsolve(@(t)fun2_(t,a(i)),t0);

p1(1,i)=real(sol(1));

p2(1,i)=real(sol(2));

p3(1,i)=real(sol(3));

p4(1,i)=real(sol(4));

p5(1,i)=real(sol(5));

end

figure

plot(a,p5);

function w=fun2_(t,a)

w1=exp(17.86806-4.60517*t(2))+exp(27.63102+log(t(1))-4.60517*t(2))+exp(35.8973+2*log(t(1))-4.60517*t(2))+exp(42.73598+3*log(t(1))-4.60517*t(2))+exp(47.80167+4*log(t(1))-4.60517*t(2))+exp(46.48919+5*log(t(1))-4.60517*t(2))+exp(-0.27631-2.30259*t(2))+exp(-17.96016)+exp(-46.69643+2.30259*t(2))+exp(-76.81424+4.60517*t(2))+exp(10.22348-4.60517*t(2))+exp(10.22348-4.60517*t(2))+exp(19.52592+3*log(t(1))-2.30259*t(2))+exp(-11.28267+2*log(t(1)))+exp(-42.55177+log(t(1))+2.30259*t(2))-t(4);

w2=exp(log(t(1)))+exp(22.56533+log(t(1))-2.30259*t(2))+exp(27.63102+log(t(1))-4.60517*t(2))+2*exp(35.8973+2*log(t(1))-4.60517*t(2))+3*exp(42.73598+3*log(t(1))-4.60517*t(2))+4*exp(47.80167+4*log(t(1))-4.60517*t(2))+5*exp(46.48919+5*log(t(1))-4.60517*t(2))+3*exp(19.52592+3*log(t(1))-2.30259*t(2))+2*exp(-11.28267+2*log(t(1)))+exp(-42.55177+log(t(1))+2.30259*t(2))-t(5)-2*a;

w3=exp(log(t(3)))-t(5);

w4=exp(22.56533+log(t(1))-2.30259*t(2))+exp(-2.30259*t(2))-exp(-32.65066+2.30259*t(2))+2*exp(17.86806-4.60517*t(2))-2*exp(log(t(3)))+2*exp(27.63102+log(t(1))-4.60517*t(2))+2*exp(35.8973+2*log(t(1))-4.60517*t(2))+2*exp(42.73598+3*log(t(1))-4.60517*t(2))+2*exp(47.80167+4*log(t(1))-4.60517*t(2))+2*exp(46.48919+5*log(t(1))-4.60517*t(2))+exp(-0.27631-2.30259*t(2))-exp(-46.69643+2.30259*t(2))-2*exp(-76.81424+4.60517*t(2))+2*exp(10.22348-4.60517*t(2))+exp(19.52592+3*log(t(1))-2.30259*t(2))-exp(-42.55177+log(t(1))+2.30259*t(2))-2*exp(-48.21613+4.60517*t(2))-exp(-26.20342+2.30259*t(2));

w5=exp(17.86806-4.60517*t(2)+1.5*(-32.65066+2.30259*t(2))+0.25*log(t(3)))-10.^(-16.86);

w=[w1;w2;w3;w4;w5];

end