Question

急！！！！！哪位大哥会matlab 帮忙写个程序，用牛顿迭代法求方程X^3-3*x-1=0在2附近的解

把完整的语句复制过来，要求计算结果保留四位有效数字，小弟刚接触矩阵实验室，一输命令就错，很是沮丧

Answer 1

syms x
x0=2;
f=x^3-3*x-1;
eps=1e-6;
maxcnt=1000;

fx=diff(f,x);
x1=x0;
cnt=1;
while cnt<=maxcnt
x2=x1-subs(f/fx,x,x1);
if abs(x1-x2)<eps
break;
end
[cnt,x1,x2] %迭代次数，迭代前，迭代后
x1=x2;
cnt=cnt+1;
end

subs(f,x,x2)

Answer 2

#include<math.h>
main()
{
float f,fd,x,x1;
x=0;
do
{
f=x*x*x+9.2*x*x+16.7*x+4;
fd=3*x*x+18.4*x+16.7;
x1=x-f/fd;
x=x1;
} while(fabs(f/fd)>1e-5);
printf("x=%.6f",x1);
}