ep6_f4.m

此代码展示了Euler计算方法与差分方法

%  Ep6_f4: << 实验六 >>  人寿保险的年利率(%)
%  Designed by FGH

%  函数参数: M为投保人寿命。
function  rate= ep6_f4(M)
N0= 40;
N= 60-N0;
M= M-N0;
F(M+1)= 0;
F(1)= 0;
r= 1/100;
while F(M+1) < 20000
   r= r+0.0001;
   for k= 1:M
      if k <= N
         p= 1540;
      else
         p= 0;
      end
      if mod(k,5) == 0
         q= 200*k+3000;
      else
         q= 0;
      end
      F(k+1)= (1+r)*F(k) + p - q;
   end
end   

rate= 100*r;