matlab吧 关注：291,850贴子：1,658,283
- 3回复贴，共1页
不感兴趣
开通SVIP免广告
clc; clear;
f = @(x) x.^3+1./x.^2+1./(7-x.^2-2./x) - 2;
df = @(x) (2*x - 2./x.^2)./(2./x + x.^2 - 7).^2 + 3*x.^2 - 2./x.^3;
x0 = 0.5;
x1 = x0 - f(x0)/df(x0);
c = 1;
while abs(x1-x0) > 1e-10 && c < 200
x0 = x1;
x1 = x0 - f(x0)/df(x0);
c = c + 1;
end
y0 = 2 - x0^3 - 1/x0^2;
disp([x0, y0]);
牛顿迭代法求得数值解：
x=-0.6740 y=0.1051
f = @(x) x.^3+1./x.^2+1./(7-x.^2-2./x) - 2;
df = @(x) (2*x - 2./x.^2)./(2./x + x.^2 - 7).^2 + 3*x.^2 - 2./x.^3;
x0 = 0.5;
x1 = x0 - f(x0)/df(x0);
c = 1;
while abs(x1-x0) > 1e-10 && c < 200
x0 = x1;
x1 = x0 - f(x0)/df(x0);
c = c + 1;
end
y0 = 2 - x0^3 - 1/x0^2;
disp([x0, y0]);
牛顿迭代法求得数值解：
x=-0.6740 y=0.1051
