ep6_p1.m

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

%  Ep6_p1: << 实验六 >>  个人住房抵押贷款 --- 月还款额
%  Designed by FGH

A0= 10000;
sdr= [6.12 6.39  6.66 7.20 7.56 7.56];
sdk= [0.5 1 3 5 6 20];
fdk= 1:9;
fdr= 0.135*(fdk-1)+6.12;

for k=10:20
   fdk(k)= k; 
   fdr(k)= 7.20;   
end

m(1)= 0;
s(1)= 10612.00;
r= fdr/12/100;
for k=2:20
   m(k)= ( A0*r(k)*(1+r(k))^(k*12) )/( ((1+r(k))^(k*12)) - 1 );
   s(k)= k*12*m(k);
end

tb= sprintf('         贷款%d元还款表\n\n',A0);
tb= [tb sprintf('  年   月  年利率  月利率  月还款额  本息总额\n\n')];
for k=1:20
   tb= [tb sprintf('  %2d  %3d  %6.4f  %6.4f  %8.4f  %8.2f\n',...
          				k,12*k,fdr(k),fdr(k)/12,m(k),s(k))];
end

tb

plot(sdk,sdr,fdk,fdr)

pause
close
clear all