罚函数法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