不得不说,这个方程是个变态,很难找到真实解的区间。经过几次尝试之后可以绘得了一个图clearr=linspace(-5e1,-2.5,1000);rr=3.99-0.32*(1+r).^(-1)-0.28498*(1+r).^(-2)-0.24899*(1+r).^(-3)-0.27558*(1+r).^(-4)-0.3132*(1+r).^(-5)-0.3547*(1+r).^(-6)-0.40035*(1+r).^(-7)-0.45036*(1+r).^(-8)-0.50506*(1+r).^(-9)-0.56463*(1+r).^(-10)-0.63936*(1+r).^(-11)-0.69947*(1+r).^(-12)-0.81354*(1+r).^(-12./r);plot(r,rr)holdonplot([-50,0],[6.73,6.73],'r')如图所示,两线的交点就是解。将图的精度取高一点,就可以直接从图中得到解了,或者再用数值解法得到更精确的解。
|