main.m

来自「罚函数法matlab程序.希望对大家有用.」· M 代码 · 共 21 行

%主函数main.m
clc
clear
a=input('请输入初始a：');
b=input('请输入初始b：');
e=input('请输入控制精度：');

x=b-((b-a)/(f1(b)-f1(a)))*f1(b);
count=1;
while abs(a-b)>e
    a=b;
    b=x;
    x=b-(b-a)*f1(b)/(f1(b)-f1(a));
    count=count+1;
end
fprintf('割线法的解：\n');
x
fprintf('牛顿迭代法迭代次数：\n');
count