#include?
#include?
#include?
using?namespace?std;
void?twoFF（）
{cout<<“所求方程為x*x*x-x*x-2*x+1=0\n“;
double?ab;
cout<<“請輸入求根區(qū)間：\n“;
cout<<“a=“;
cin>>a;
cout<<“b=“;
cin>>b;
double?cd;
cout<<“請輸入要求達到的精度：\n“;
cout<<“精度為“;
cin>>d;
cout<<“a=0\n“;
cout<<“b=1\n“;
int?n=0;
while（c>d?||?c<-d）
{if（（a*a*a-a*a-2*a+1）*（（a+b）*（a+b）*（a+b）/8-（a+b）*（a+b）/4-a-b+1）>0）{a=（a+b）/2;cout<<“a=“<else?{b=（a+b）/2;cout<<“b=“<n++;
c=a-b;}
cout<<“迭代次數(shù)為“<}

void?Newton（）
{cout<<“求方程x*（x+1）*（x+1）-1=0在x=0.4附近的根\n“;
double?x=0.4m;
double?cd;
cout<<“請輸入要求達到的精度：\n“;
cout<<“精度為“;
cin>>d;
cout<<“x=0.4\n“;
int?n=0;
while（c>d?||?c<-d）{m=x;
x?=?x?-（x*（x+1）*（x+1）-1）?/（3*x*x+4*x+1）;
cout<<“x=“<c=m-x;
n++;}
cout<<“迭代次數(shù)為“<}

int?n;
double?a[100][101];
double?m[100];
void?changerow（int?g）
{int?ij;
if（a[g-1][g-1]==0）
{int?k;
for（k=g;k{if（a[k][g-1]!=0）break;
}
for（j=0;j{m[j]=a[g-1][j];
a[g-1][j]=a[k][j];
a[k][j]=m[j];}
}
for（i=g;im[i]=a[i][g-1]/a[g-1][g-1];
for（i=g;ifor（j=g-1;ja[i][j]=a[i][j]-m[i]*a[g-1][j];

for（i=0;i{cout<<“\n“;
	for（j=0;j????printf（“????%2f“a[i][j]）;
????}
cout<<“\n“;
}
void?GSXU（）
{
printf（“輸入方程個數(shù):“）;
cin>>n;
printf（“按行輸入各個方程的系數(shù):\n“）;
int?ij;
for（i=0;ifor（j=0;jcin>>a[i][j];
for（i=0;i{cout<<“\n“;
	for（j=0;j????printf（“????%2f“a[i][j]）;}
cout<<“\n“;
for（i=1;i{changerow（i）;}

double?x[100];
int?r=n-1;
x[n-1]=a[n-1][n]/a[n-1][n-1];
double?p=0q;
for（i=n-2;i>=0;i--）
{
for（j=n-1;j>=r;j--）
{q=a[i][j]*x[j];
p=p+q;
}
x[i]=（a[i][n]-p）/a[i][i];
r--;
p=0;}
printf（“方程組解為“）;
cout<<“X=“<<“（“;
for（i=0;icout<cout<cout<<“）“<}

void?Jacobi（）
{double?a[4][4]?x1[100]x2[100]x3[100];
x1[0]=0;x2[0]=0;x3[0]=0;
printf（“按行輸入三元方程的各個系數(shù):\n“）;
int?ij;
for（i=0;i<3;i++）
for（j=0;j<4;j++）
cin>>a[i][j];
x1[1]=a[0][3]/a[0][0];
x2[1]=a[1][3]/a[1][1];
x3[1]=a[2][3]/a[2][2];
double?c=1;
for（i=1;c>0.000001;i++）
{
x1[i]=x1[1]-a[0][1]/a[0][0]*x2[i-1]-a[0][2]/a[0][0]*x3[i-1];
x2[i]=x2[1]-a[1][0]/a[1][1]*x1[i-1]-a[1][2]/a[1][1]*x3[i-1];
x3[i]=x3[1]-a[2][0]/a[2][2]*x1[i-1]-a[2][1]/a[2][2]*x2[i-1];
c=fabs（x1[i]-x1[i-1]）+fabs（x2[i]-x2[i-1]）+fabs（x3[i]-x3[i-1]）;
cout<}
}

double?x[100]y[100];
void?input（int?n）
{
int?i;
for（i=0;i{cout<<“x[“<scanf（“%lf“&x[i]）;
cout<<“y[“<scanf（“%lf“&y[i]）;}
}
void?Lagrange（）
{
int?m;
printf（“輸入點的個數(shù)：\n“）;
scanf（“%d“&m）;
double?v;
printf（“輸入所求X的值：\n“）;
scanf（“%lf“&v）;
input（m）;
int?ij;
double?ab=1cd=1ef=1gh=1pq=0;
for（i=0;i{??for（j=0;j???{a=x[i]-x[j];
????b=b*a;
	c=v-x[j];
	d=d*c;}
???for（j=i+1;j???{e=x[i]-x[j];
????f=f*e;
	g=v-x[j];
	h=h*g;}
p=d*h*y[i]/b/f;
cout<<“L“<