十大解法c语言Word格式文档下载.docx
《十大解法c语言Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《十大解法c语言Word格式文档下载.docx（12页珍藏版）》请在冰豆网上搜索。
i,n;
x[20],y[20],xx,yy;
printf（"
Input
n:
"
）;
scanf（"
%d"
&
n）;
if（n>
=20）
{printf（"
Error!
The
value
of
n
must
in
（0,20）."
getch（）;
1;
if（n<
=0）
x[%d]:
i）;
%f"
x[i]）;
\n"
y[%d]:
y[i]）;
xx:
xx）;
yy=lagrange（x,y,xx,n）;
x=%f,y=%f\n"
xx,yy）;
2.牛顿插值多项式，用于离散数据的拟合
void
difference（float
*y,int
n）
*f;
k,i;
f=（float
for（k=1;
k<
=n;
k++）
f[0]=y[k];
k;
f[i+1]=（f[i]-y[i]）/（x[k]-x[i]）;
y[k]=f[k];
return;
difference（x,（float
*）y,n）;
yy=y[20];
for（i=n-1;
i>
=0;
i--）
yy=yy*（xx-x[i]）+y[i];
NewtonInter（%f）=%f"
3.高斯列主元消去法,求解其次线性方程组
code#include<
math.h>
#define
N
20
n,i,j,k;
mi,tmp,mx;
a[N][N],b[N],x[N];
\nInput
N）
input
should
in（0,N）!
=0）
Now
a（i,j）,i,j=0...%d:
n-1）;
n;
a[i][j]）;
b（i）,i,j=0...%d:
b[i]）;
n-2;
for（j=i+1,mi=i,mx=fabs（a[i][j]）;
n-1;
if（fabs（a[j][i]）>
mx）
mi=j;
mx=fabs（a[j][i]）;
if（i<
mi）
tmp=b[i];
b[i]=b[mi];
b[mi]=tmp;
for（j=i;
tmp=a[i][j];
a[i][j]=a[mi][j];
a[mi][j]=tmp;
for（j=i+1;
tmp=-a[j][i]/a[i][i];
b[j]+=b[i]*tmp;
for（k=i;
a[j][k]+=a[i][k]*tmp;
x[n-1]=b[n-1]/a[n-1][n-1];
for（i=n-2;
i--）
x[i]=b[i];
x[i]-=a[i][j]*x[j];
x[i]/=a[i][i];
Answer:
\n
x[%d]=%f\n"
i,x[i]）;
0;
#include<
NUMBER
Esc
0x1b
Enter
0x0d
A[NUMBER][NUMBER+1]
ark;
flag,n;
exchange（int
r,int
k）;
max（int
message（）;
{
x[NUMBER];
r,k,i,j;
char
celect;
clrscr（）;
\n\nUse
Gauss."
\n\n1.Jie
please
press
Enter."
\n\n2.Exit
Esc."
celect=getch（）;
if（celect==Esc）
exit（0）;
\n\n
n="
\n\nInput
matrix
A
and
B:
for（i=1;
a%d1--a%d%d
b%d:
i,i,n,i）;
for（j=1;
=n+1;
j++）
A[i][j]）;
k++）
ark=max（k）;
if（ark==0）
\n\nIt'
s
wrong!
else
if（flag!
=k）
exchange（flag,k）;
for（i=k+1;
for（j=k+1;
A[i][j]=A[i][j]-A[k][j]*A[i][k]/A[k][k];
x[n]=A[n][n+1]/A[n][n];
for（
k=n-1;
k>
=1;
k--）
me=0;
me=me+A[k][j]*x[j];
x[k]=（A[k][n+1]-me）/A[k][k];
\n\nx%d=%f"
k）
i;
A[0][i]=A[r][i];
A[r][i]=A[k][i];
A[k][i]=A[0][i];
temp=0;
for（i=k;
if（fabs（A[i][k]）>
temp）
temp=fabs（A[i][k]）;
flag=i;
temp;
message（）
Go
on
Exit
Esc!
switch（getch（））
case
Enter:
main（）;
Esc:
default:
error!
4.龙贝格求积公式，求解定积分
f（x）
（sin（x）/x）
MAX
20
a
2
b
4
e
0.00001
LBG（float
p,float
q,int
sum=0,h=（q-p）/n;
for
（i=1;
sum+=f（p+i*h）;
sum+=（f（p）+f（q））/2;
return（h*sum）;
n=N,m=0;
T[MAX+1][2];
T[0][1]=LBG（a,b,n）;
n*=2;
for（m=1;
m<
MAX;
m++）
m;
T[i][0]=T[i][1];
=m;
T[i][1]=T[i-1][1]+（T[i-1][1]-T[i-1][0]）/（pow（2,2*m）-1）;
if（（T[m-1][1]<
T[m][1]+e）&
&
（T[m-1][1]>
T[m][1]-e））
Answer=%f\n"
T[m][1]）;
;
code5.牛顿迭代公式，求方程的近似解
100
PS
1e-5
TA
Newton（float
（*f）（float）,float（*f1）（float）,float
x0
）
x1,d=0;
k=0;
do
x1=
x0-f（x0）/f1（x0）;
if（（k++>
N）||（fabs（f1（x1））<
PS））
\nFailed!
exit（）;
d=（fabs（x1）<
1?
x1-x0:
（x1-x0）/x1）;
x0=x1;
x（%d）=%f\n"
k,x0）;
while（（fabs（d））>
PS&
fabs（f（x1））>
TA）
x1;
f（float
x）
x*x*x+x*x-3*x-3;
f1（float
3.0*x*x+2*x-3;
f（float）;
f1（float）;
x0,y0;
x0:
x0）;
x（0）=%f\n"
x0）;
y0=Newton（f,f1,x0）;
\nThe
root
is
x=%f\n"
y0）;
6.
牛顿-科特斯求积公式，求定积分
NC（a,h,n,r,f）
（*a）[];
h;
n,f;
*r;
nn,i;
ds;
1000||n<
2）
if
（f）
Faild!
Check
1<
n<
1000!
n）;
return（-1）;
if（n==2）
*r=0.5*（（*a）[0]+（*a）[1]）*（h）;
return（0）;
（n-4==0）
*r=0;
*r=*r+0.375*（h）*（（*a）[n-4]+3*（*a）[n-3]+3*（*a）[n-2]+（*a）[n-1]）;
if（n/2-（n-1）/2<
nn=n;
else
nn=n-3;
ds=（*a）[0]-（*a）[nn-1];
for（i=2;
=nn;
i=i+2）
ds=ds+4*（*a）[i-1]+2*（*a）[i];
*r=ds*（h）/3;
nn）
h,r;
n,ntf,f;
a[16];
the
x[i]（16）:
=15;
a[i]）;
h=0.2;
f=0;
ntf=NC（a,h,n,&
r,f）;
if（ntf==0）
\nR=%f\n"
r）;
Wrong!
Return
code=%d\n"
ntf）;
7.雅克比迭代，求解方程近似解
0.00001
i,j,k;
t;
a[N][N],b[N][N],c[N],g[N],x[N],h[N];
dim
0<
N!
a[i,j],i,j=0…%d:
c[i],i=0…%d:
c[i]）;
b[i][j]=-a[i][j]/a[i][i];
g[i]=c[i]/a[i][i];
h[j]=g[j];
for（k=0;
if（j==k）
continue;
h[j]+=b[j][k]*x[k];
t=0;
if（t<
fabs（h[j]-x[j]））
t=fabs（h[j]-x[j]）;
x[j]=h[j];
e）
x_i=\n"
i++）
after
%d
repeat
return\n"
MAX）;
8.秦九昭算法
qin（float
a[],int
n,float
r=0;
for（i=n;
r=r*x+a[i];
r;
a[50],x,r=0;
n,i;
frequency:
while（n<
1）;
value:
x）;
r=qin（a,n,x）;
9.幂法
3
x[n]={0,0,1};
a[n][n]={{2,3,2},{10,3,4},{3,6,1}};
y[n];
xm,oxm;
oxm=0;
N;
for（j=0
升级会员 