function m4=y1(t,u)
format long
m4=0;
n=9;
for k=1:n
m1=0;
for j=0:k+2
p=2+j+k;
q=(p*log(2))/t;
a=[-0.5,-5,85-q,450+10*q,-1000-200*q];subs(a);
y=roots(a);z=find(real(y)<0);x=y(z);
c1=x(2)*(20+x(1))/(20*(x(2)-x(1)));
c2=x(1)*(20+x(2))/(20*(x(1)-x(2)));
a1=(c1*exp(x(1)*u)+c2*exp(x(2)*u))/q;
m=(-1)^j*factorial(k+2)/(factorial(j)*factorial(k+2-j));
m1=m1+m*a1;
end
m2=(-1)^(n-k)*k^n/(factorial(k)*factorial(n-k))*log(2)/t;
m3=factorial(2*k+4)/(factorial(k+2)*factorial(k+1))*m1;
m4=m2*m3+m4;
end
m4;
function m=y2(T,u)
m=0;
for i=1:100
t(i)=i*T/100;
s=exp((-1)*0.05*t(i))*(y1(t(i),u)-y1(t(i-1),u));
m=m+s;
end
m;
MATLAB默认的工作目录为安装目录下的work目录中，所以你编写的自定义函数都要放在work目录中陵世，告塌并且，MATLAB调用函数不是以函数名为标志，而是以自定义函数保存的文件名标志，鉴于此尺友肢，你可以把函数文件保存为函数名以方便识别和调用。