function[m,ss,pesai]=f(V,theta,kapa,delta,epsilon)
tic
m=theta+(V-theta)*exp(-kapa*delta);
ss=V*epsilon^2*exp(-kapa*delta)*(1-exp(-kapa*delta))/kapa+theta*epsilon^2/(2*kapa)*(1-exp(-kapa*delta))^2;
pesai=ss/m^2;
toc
在工作空间使用：
f(1,3,1,5,1)
只显示m的值，其结果与
m=f(1,3,1,5,1)相同
如果要接收所有返回值，则应使用：
[m,ss,pesai]=f(1,3,2,5,4)
结果显示：
m =
2.9999
ss =
11.9993
pesai =
1.3333
