как написать этот алгоритм? - C++ - Ответ 2466470

#include <stdio.h>
#include <math.h>
void main()
{
    const int MAX=101;
    const int n=100;
    double h=1./n;
    double x=1.5;
    double M_PI=3.1415926;
    double alfa[MAX],beta[MAX],y[MAX],A[MAX],B[MAX],C[MAX],F[MAX];
     C[0]=1; B[0]=0; F[0]=0;
     C[n]=1; F[n]=0; A[n]=0;
    for(int i=1;i<=n-1;i++)
    {   x=i*h; 
        A[i]=1./(h*h);
        B[i]=A[i];
        C[i]=2./(h*h);
        F[i]=4.*M_PI*M_PI*sin(2.*M_PI*x);
    }
    alfa[1]=B[0]/C[0];
    beta[1]=F[0]/C[0];
//  printf("alfa[%i]=%0.4lf | beta[%i]=%0.4lf\n",1,alfa[1],1,beta[1]);
    for(int j=1;j<=n-1;j++)
  
    {
        alfa[j+1]=B[j]/(C[j]-alfa[j]*A[j]);
        beta[j+1]=(F[j]+A[j]*beta[j])/(C[j]-alfa[j]*A[j]);
//      printf("alfa[%i]=%0.4lf | beta[%i]=%0.4lf\n",j+1,alfa[j+1],j+1,beta[j+1]);
    }
    y[n]=(F[n]+A[n]*alfa[n])/(C[n]-alfa[n]*A[n]);
//    printf("y[%i]=%0.4lf\n",n-1,y[n-1]);
    for(int k=n-1; k>=0;k--)
    
    {
     y[k]=alfa[k+1]*y[k+1]+beta[k+1];
     printf("%8.4lf",y[k]);
    }
 printf("\n");
    for(int k=n-1;k>=0;k--) printf("%8.4lf",sin(2*M_PI*k*h));
 
}