ep6_f3.m

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

%  Ep6_f3: << 实验六 >>  养老保险的月利率(%)
%  Designed by FGH

function  rate= ep6_f3(p,q,N,M)
%  (6.15)
pol(1)= 1;
pol(N+1)= -1-q/p;
pol(M+1)= q/p;
x= roots(pol);
n1= 1;
for n= 1:M
   if isreal(x(n)) & x(n)-1 > 0.00001
      r(n1)= x(n)-1;
      n1= n1+1;
   end
end
rate= 100*r;