Advertisement
Guest User

Untitled

a guest
Dec 6th, 2016
77
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C 3.38 KB | None | 0 0
  1. #include <stdio.h>///
  2. #include "math.h"///
  3. #define N 6///
  4. main()//
  5. {///
  6.       double x[N]={0.180, 0.185, 0.190, 0.195, 0.200, 0.205};///
  7.       double y[N]={5.61543, 5.46693, 5.32634, 5.19304, 5.06649, 4.94619};///
  8.       double X[4];///
  9.       double c, s1, sum;//
  10.       double A[4][4], B[4];//
  11.       int i,j,k,s;//
  12.       for(k=1; k<=3; k++)//
  13.       {///
  14.                if(k==1)printf("lin\n");
  15.                if(k==2)printf("kvadratnaya\n");
  16.                if(k==3)printf("kub\n");
  17.             //  
  18.                for(i=0; i<=k; i++)//
  19.                {
  20.                         for(j=0; j<=k; j++)//
  21.                         {///
  22.                                  A[i][j]=0;//
  23.                                  for(s=0;s<N;s++) A[i][j]=A[i][j]+pow(x[s], i+j);//
  24.                                  }//
  25.                                  B[i]=0;//
  26.                                  for(s=0;s<N;s++) B[i]=B[i]+y[s]*pow(x[s], i);//
  27.                                  }//
  28.                                  for(s=0; s<k; s++)//
  29.                                  {//
  30.                                           for(i=s+1; i<=k;i++)//
  31.                                           {//
  32.                                                      c=A[i][s]/A[s][s];//
  33.                                                      B[i]=B[i]-(c*B[s]);//
  34.                                                      for(j=s+1; j<=k; j++) A[i][j]=A[i][j]-c*A[s][j];//
  35.                                                      }//
  36.                                                      }//
  37.                                           X[k]=B[k]/A[k][k];//
  38.                                           for(i=k-1; i>=0; i--)//////////////////////
  39.                                           {//
  40.                                                      c=B[i];//
  41.                                                      for(j=i+1; j<=k; j++){//
  42.                                                                 c=c-A[i][j]*X[j];}// скобку добавить
  43.                                                                 X[i]=c/A[i][i];//
  44.                                                                 }//
  45.                                                                 sum=0;//
  46.                                                                 for(s=0;s<N; s++)//
  47.                                                                              {s1=0;//
  48.                                                                              for(i=0; i<=k;i++)
  49.                                                                               s1=s1+X[i]*pow(x[s],i);//
  50.                                                                              s1=(s1-y[s])*(s1-y[s]);//
  51.                                                                              sum=sum+s1;//
  52.                                                                              }//
  53.                                           sum=sqrt(sum/(N-1));//
  54.                                           printf("sred %f\n", sum);//
  55.                                          
  56.                                           for(i=0; i<=k; i++)
  57.                                           printf("X[%d]= %f\n",i, X[i]);
  58.                                           int a;
  59.                                           scanf("%d", &a);
  60.                                           }return 0;}
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement