分别用牛顿法和割线法求解方程 x^3-6x^2+9x-2=0在区间[3,4]上的近似根.要求满足精度|x*-xk|
题目
分别用牛顿法和割线法求解方程 x^3-6x^2+9x-2=0在区间[3,4]上的近似根.要求满足精度|x*-xk|
答案
%
clc; clear all;
global fnq dfnq
fnq = @(x) x^3 - 6*x^2 + 9*x - 2;
dfnq = @(x) 3*x^2 - 12*x + 9;
tol = (1/2)*10^-4;x0 = 3.5;gmax = 1e3;x01 = 3;x02 = 4;[k,xk,yk,piancha]=newtonqx(x0,tol,gmax);
fprintf('\n牛顿法:%.5f\n',xk);
[k,xk,yk,piancha]=gexian(x01,x02,tol,gmax);
fprintf('\n割线法:%.5f\n',xk);% 牛顿法
function [k,xk,yk,piancha]=newtonqx(x0,tol,gmax)
global fnq dfnqx(1)=x0;for i=1:gmax
x(i+1)=x(i)-fnq(x(i))/(dfnq(x(i)+eps));
piancha=abs(x(i+1)-x(i));i=i+1;
xk=x(i);yk=fnq(x(i));
[(i-1) xk yk piancha];
if(pianchagmaxdisp('超过最大迭代次数')k=i-1;xk=x(i);yk=fnq(x(i));
[i-1 xk yk piancha];return;end% 割线法function [k,xk,yk,piancha]=gexian(x01,x02,tol,gmax)
global fnq dfnq
x(1)=x01;x(2)=x02;
for i=2:gmax
u(i)=fnq(x(i))*(x(i)-x(i-1));
v(i)=fnq(x(i))-fnq(x(i-1));
x(i+1)=x(i)-u(i)/(v(i));
piancha=abs(x(i+1)-x(i));i=i+1;xk=x(i);yk=fnq(x(i));
if(piancha>
举一反三
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
- 英语修改病句!紧急紧急...速度!
- 一根木料截去100厘米后剩下的木料是原来的1.2倍段30厘米,剩下多少
- 如果不等式(2a-1)x1,化简丨1-2a丨+丨a-3丨.
- 关于初三“液体的压强”一道计算题!
- 如图,正方形abcd边长为6.菱形efgh的三个顶点e,g,h分别在正方形abcd的边ab,cd,
- 三个不同的有理数,既可以表示为1,a+b,a,又可以表示为0,b分之a,试求a的2012次方+b的2012次方的值.
- only when you realize the importance of foreign languages_____them well.
- 印度的水稻,小麦与降水,地形的关系
- 他肯定来参加你的生日聚会 翻译成英语
- 雷锋日记读书笔记600字丶急!
热门考点
- 发挥你的想象,续写下列句子
- 水结冰后继续降温,会怎样?有没有-10度的冰?
- 一个人在两面墙之间,大喊一声,同时听到回声,则两面墙的最短距离为_______
- 食堂买了160千克白菜比买的芹菜的2有2分之1倍多30千克食堂买了芹菜多少千克
- 观察感受新年的习俗,写作文2篇,600字
- 解方程:x2+x+1x2+1 + 2 x2+x+2x2+x+1= 19/6.
- 已知函数f(x)=Asin(2x+φ)(A>0,0
- 物体从静止开始,以2m 每二次方秒的加速度运动,到第六秒内的平均速度是多少?位移是多少?
- 设向量a,b的夹角为135°,且|a|=√2,|b|=2,c=a+xb(其中x∈R,)当|a+xb|取最小值时,求a+xb与b夹角的大小
- 《少年中国说》的朗诵文章