Question

牛顿法 matlab 求 x^3-x-1=0 在 x=1.5附近的根

就是这样。。。 牛顿法 matlab 求 x^3-x-1=0 在 x=1.5附近的根
还有 ，用二分法求方程x^2-x-1=0 的正根 ，要求精确到小数点后一位 同样 matlab实现

Answer 1

　　%牛顿迭代法 解方程y=x.^3-x-1
　　x=1.5;
　　format long;
　　x1=x-func1_1(x)/func1_1_1(x);
　　if(abs(x1)<1.5)
　　delt=abs(x1-x);
　　else
　　delt=abs((x1-x)/x1);
　　end
　　while(delt>1e-6|abs(func1_1(x))>1e-6)
　　x=x1; x1=x-func1_1(x)/func1_1_1(x);
　　if(abs(x1)<1.5)
　　delt=abs(x1-x);
　　else
　　delt=abs((x1-x)/x1);
　　end
　　if func1_1(x1)==0
　　break
　　end
　　end

　　disp('解方程y=x.^3-x-1,牛顿迭代法结果')
　　x1

　　编辑函数

　　function y=func1_1(x)
　　y=x.^3-x-1;

　　function y=func1_1_1(x)
　　y=3*x.^2-1

　　% 对分法
　　delta=10e-5;
　　a=1;
　　b=2;
　　fa=func2_1(a);
　　fb=func2_1(b);
　　n=1;
　　while(1)
　　if(fa*fb>0)
　　break;
　　end
　　x=(a+b)/2;
　　fx=func2_1(x);
　　if(abs(fx)<delta)
　　break;
　　elseif(fa*fx<0)
　　b=x;fb=fx;
　　else
　　a=x;fa=fx;
　　end
　　if(b-a<eps)
　　break;
　　end
　　n=n+1;
　　end
　　disp('对分法结果');
　　x

　　% 调用roots函数求解结果并作图
　　x=linspace(-5,5);
　　y=x.^2-x-1;
　　plot(x,y)
　　p=[1 -1 -1];
　　disp('调用roots函数求解')
　　x=roots(p)

　　编辑函数

　　function y=func2_1(x)
　　y=x.^2-x-1;

　　function y=func2_1_1(x)
　　y=2*x-1