for i=1:n
mα=mα(1:i);
c=c(1:i);b=b(1:i);
W=W(1:i);μ=μ (1:i);φ=φ(1:i);
α=α(1:i);Q=Q(1:i);e=e(1:i);
Fs1=Fs1+sum(cb+(W-μb)tan(φ)/mα)/(sum(Wsinα)+sum(Q*e/R))
Fs(i)=Fs1;
end这个也运行不了啊