clear all;
a0=150;b=1;n=15;
%a1=500;b1=4.5;n1=15;
%a2=500;b2=45;n2=20;
x0=1;
[x1,fv1,ef1,out1]=fzero(@rate,x0,[],a0,n,b,1);%對第一題小張夫婦買方月利率問題進行求解
a1=500;b1=4.5;n1=15;
[x2,fv2,ef2,out2]=fzero(@rate,x0,[],a1,n1,b1,1);%對二問第一家銀行的利率求解
a2=500;b2=45;n2=20;
[x3,fv3,ef3,out3]=fzero(@rate,x0,[],a2,n2,b2,-1);%對二問第二家銀行的利率求解
format long;% 顯示格式
disp('第一題小張夫婦買方月利率問題進行求解')
x1%輸出其中x1匙小張夫婦買房的月利率
disp('第二問第一家銀行的利率求解')
disp('第一家月利率')
x2%x2是第一家銀行的月利率
disp('第一家年利率')
X2=x2*12
disp('第二問第二家銀行的利率求解')
disp('第二家月利率')
X3=x3/12
disp('第二家年利率')
x3%x3市第二家銀行的年利率
if x2*12<x3
    disp('diyi jia youhui ')
else
    disp('di er jia ')
end
%y=inline('1/r-1/(r*(1+r)^180)-150','r');
%r=-1:0.01:1;
%y_char=vectorize(y);
%Y=feval(y_char,r);
%clf,plot(r,Y,'r');hold on,plot(r,zeros(size(r)),'k') ;
%xlabel('r');ylabel('y(r)'),hold off
%xlabel('r');ylabel('y(r)'),hold off
%clf,plot(r,Y,'r');hold on,plot(r,zeros(size(r)),'k') ;
%xlabel('r');ylabel('y(r)'),hold off
%clf,plot(r,Y,'r');hold on,plot(r,zeros(size(r)),'k') ;
%xlabel('r');ylabel('y(r)'),hold off
%clf,plot(r,Y,'r');hold on,plot(r,zeros(size(r)),'k') ;
%xlabel('r');ylabel('y(r)'),hold off