牛顿法7_2.cpp

来自「牛顿法求解非线性方程组」· C++ 代码 · 共 47 行

// 牛顿法7_2.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"
#include <iostream>
#include <math.h>
#include <conio.h>

using namespace std;
int _tmain(int argc, _TCHAR* argv[])
{
	float max(float s,float n);
	int i,j;
	float x[2],k,l,m,s,n;

	cout<<"【牛顿法求解非线性方程组】"<<endl<<endl;	
	cout<<" 请输入X1=";
	cin>>x[1];	
	cout<<" 请输入X2=";
	cin>>x[2];
	cout<<endl<<"计算结果："<<endl;
	for(j=0;j<10;j++)
	{
		l=18*x[0]*x[0]*x[0]+6*x[1]*x[1]*x[1]-6*x[0]*x[1]*x[1];
		k=18*x[0]*x[0]*x[0]-6*x[0]*x[1]*x[1]+6*x[0]*x[0]*x[1]*x[1]-2*x[0]*x[0]*x[0]*x[0]-2*x[0];
		m=
		6*x[0]*x[0]*x[0]*x[0]+3*x[1]*x[1]*x[0]*x[0]+3*x[1]*x[1]*x[1]*x[1]-3*x[0];
		s=k/l;
		n=m/l;
		x[0]=x[0]-s;
		x[1]=x[1]-n;
		cout<<"x1="<<x[1]<<"   x2="<<x[2]<<endl;
	}
	cout<<endl<<"按任意键继续...";
	getch();

	return 0;
}

float max(float s,float n)
{
	float m;
	if(s<0) s=-1*s;
	if(n<0) n=-1*n;
	m=s>n?s:n;
	return(m);
}