main.m

罚函数法matlab程序.希望对大家有用.

%主函数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