function f=fun(x)
%f=[14.9*(1+x)^180+0.1*(1+x-(1+x)^180)/x,2682*(1+x)^179+0.1*((1+x)^180-180*x*(1+x)^179-1)/(x^2)]
f=[0.1*((1+x)^180-1)-15*x*(1+x)^180,18*(1+x)^179-15*((1+x)^180+x*180*(1+x)^179)];%在網上找的公式