Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <stdio.h>///
- #include "math.h"///
- #define N 6///
- main()//
- {///
- double x[N]={0.180, 0.185, 0.190, 0.195, 0.200, 0.205};///
- double y[N]={5.61543, 5.46693, 5.32634, 5.19304, 5.06649, 4.94619};///
- double X[4];///
- double c, s1, sum;//
- double A[4][4], B[4];//
- int i,j,k,s;//
- for(k=1; k<=3; k++)//
- {///
- if(k==1)printf("lin\n");
- if(k==2)printf("kvadratnaya\n");
- if(k==3)printf("kub\n");
- //
- for(i=0; i<=k; i++)//
- {
- for(j=0; j<=k; j++)//
- {///
- A[i][j]=0;//
- for(s=0;s<N;s++) A[i][j]=A[i][j]+pow(x[s], i+j);//
- }//
- B[i]=0;//
- for(s=0;s<N;s++) B[i]=B[i]+y[s]*pow(x[s], i);//
- }//
- for(s=0; s<k; s++)//
- {//
- for(i=s+1; i<=k;i++)//
- {//
- c=A[i][s]/A[s][s];//
- B[i]=B[i]-(c*B[s]);//
- for(j=s+1; j<=k; j++) A[i][j]=A[i][j]-c*A[s][j];//
- }//
- }//
- X[k]=B[k]/A[k][k];//
- for(i=k-1; i>=0; i--)//////////////////////
- {//
- c=B[i];//
- for(j=i+1; j<=k; j++){//
- c=c-A[i][j]*X[j];}// скобку добавить
- X[i]=c/A[i][i];//
- }//
- sum=0;//
- for(s=0;s<N; s++)//
- {s1=0;//
- for(i=0; i<=k;i++)
- s1=s1+X[i]*pow(x[s],i);//
- s1=(s1-y[s])*(s1-y[s]);//
- sum=sum+s1;//
- }//
- sum=sqrt(sum/(N-1));//
- printf("sred %f\n", sum);//
- for(i=0; i<=k; i++)
- printf("X[%d]= %f\n",i, X[i]);
- int a;
- scanf("%d", &a);
- }return 0;}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement