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

				%main.m
				clear;
				e=input('请输入控制精度：');
				x0=input('请输入初始迭代值：');
				 
				x1=x0-f1(x0)/f2(x0);
				cout=1
				while abs(x1-x0)>e
				    
				        x0=x1;
				        x1=x0-f1(x0)/f2(x0);
				        cout=cout+1;
				     
				end;
				
				fprintf('牛顿法迭代次数：\n');
				cout
				fprintf('最优点是: \n');
				x1
				fprintf('方程的最小值是：\n');
				x=f(x1)