抛物线法matlab实验,数值分析：牛顿、弦歌、抛物线法的matlab实现

%m文件1

clc;clear;close;

fangcheng = 2;%指示对那个方程进行计算

%牛顿法

Endflag = 1e-8;

ComputeNum=0;%计算次数

xk=1.5;

tstart = cputime;

fprintf('牛顿法的计算结果为:\n计算次数 计算结果\n');

while 1

fxk=fun(xk,fangcheng,0);

fpxk=fun(xk,fangcheng,1);

if abs(fpxk) < Endflag

disp('出错！');

end

xkp1=xk-fxk/fpxk;

ComputeNum=ComputeNum+1;

hold on

graph1 = plot(ComputeNum,xkp1,'*');

fprintf(' %d %12.10f\n',ComputeNum,xkp1);

if abs(xkp1-xk) < Endflag

break;

else

xk=xkp1;

end

if ComputeNum >100 %防止死循环

disp('发散');

break;

end

end

tend = cputime;

disp('牛顿法的计算次数为:');

ComputeNum

disp('牛顿法的计算时间为:&