文档介绍: ,用于离散数据的拟合
C/C++ code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
2,n);
 printf("Input xx:");
 scanf("%f",&xx);
 yy=y[20];
 for(i=n-1;i>=0;i--) yy=yy*(xx-x[i])+y[i];
 printf("NewtonInter(%f)=%f",xx,yy);
 getch();
}
36
37
38
,求解其次线性方程组
C/C++ code
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
#include<>
#include <>
#define N 20
int main()
{ int n,i,j,k;
 int mi,tmp,mx;
 float a[N][N],b[N],x[N];
 printf("\nInput n:");
 scanf("%d",&n);
 if(n>N)
 { printf("The input n should in(0,N)!\n");
    getch();
    return 1;
 }
 if(n<=0)
 { printf("The input n should in(0,N)!\n");
    getch();
    return 1;
 }
 printf("Now input a(i,j),i,j=0...%d:\n",n-1);
 for(i=0;i<n;i++)
 { for(j=0;j<n;j++)
    scanf("%f",&a[i][j]);}
 printf("Now input b(i),i,j=0...%d:\n",n-1);
 for(i=0;i<n;i++)
 scanf("%f",&b[i]);
 for(i=0;i<n-2;i++)
 { for(j=i+1,mi=i,mx=fabs(a[i][j]);j<n-1;j++)
    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;j<n;j++)
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
      { tmp