ep6_f2.m

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

%  Ep6_f2: << 实验六 >>  个人住房抵押贷款 --- 小李夫妇的贷款问题
%  Designed by FGH

%  本函数的返回值rate为小李夫妇实际享受到的贷款年利率
%  函数参数: A0为贷款总额, m为每次还款额, y为贷款年数，c为每年还款次数。
function  rate= ep6_f2(A0,m,y,c)
k= y*c;
p(1)= A0;
p(2)= -m-A0;
p(k+1)= m;
x= roots(p);
r1= (x-1)*c;
n1= 1;
for n= 1:k
   if r1(n) > 0
      r(n1)= r1(n);
      n1= n1+1;
   end
end
rate= r;