table1.m

来自「几个计算」· M 代码 · 共 15 行

% function table1()
n=1;
v=0.30;
for i=0:1:60
    for j =0:0.1:6
        table(n,1)=1+(i/60);
        table(n,2)=j;
        st=exp(-v*j);
        table(n,3)=A(5,1)*(st^table(n,1)+(st-st^table(n,1))*A(8,1))+A(4,1)*((1-st^(1/table(n,1)))+(st^(1/table(n,1))-st)*A(8,1));
        table(n,4)=A(5,2)*(st^table(n,1)+(st-st^table(n,1))*A(8,2))+A(4,2)*((1-st^(1/table(n,1)))+(st^(1/table(n,1))-st)*A(8,2));
        table(n,5)=A(5,3)*(st^table(n,1)+(st-st^table(n,1))*A(8,3))+A(4,3)*((1-st^(1/table(n,1)))+(st^(1/table(n,1))-st)*A(8,3));
        n=n+1;    
    end
end
save table table;