Решение краевой задачи разностным методом

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
 
double Pi = 3.14;
int main(void) {
    double a=0;
    double b=Pi;
    double h;
    int n, N;
    printf ("type N\n");    
    scanf ("%d", &N);
    h=(b-a)/N;
    n=N+1;
    double X[n];
    double Y[n];
    double A[n][n]; 
    int i;
    int j;
    for (i=0; i<n; i++) //заполнение массива аргументов
        X[i]=a+i*h;
    for(i=0; i<n; i++){
    for(j=0; j<n; j++){
        A[i][j]=0;
    }
 }
for(i=1; i<n-1; i++){ //заполнение матрицы слау
    A[i][i+1]=1;
    A[i][i]=2*h*h-2;
    A[i][i-1]=1;
}
    A[0][0]=-1;
    A[n-1][n-1]=-1; 
    double l;
    int k;
    for(k=0; k<n; k++) //метод гаусса
    {
        for(i=k+1; i<n; i++)
        {
            l=A[i][k]/A[k][k];
            A[i][n-1]=A[i][n-1]-l*A[k][n-1];
            for(j=k; j<(n); ++j)
            {
                A[i][j]-=l*A[k][j];
            }   
        }
    }
    double f, sum;
    for(k=n-1; k>=0; --k) {
        f=0;
        for(j=k; j<n; ++j)
        {
            sum=A[k][j]*Y[j];
            f+=sum;
        }
        Y[k]=(A[k][n-1]-f)/A[k][k];
    }
    for (i=0; i<n; ++i) //вывод данных на экран
    {
        printf ("X[%d]=%lf\n", i-1, X[i]);
        printf ("Y[%d]=%lf\n", i-1, Y[i]);
    }
    return 0;
}