Прошу выручить и дописать код С++

#include <iostream> // для ввода-вывода (объекты std::cout, std::cin, манипулятор std::endl)
#include <cmath> // В стандартную математическую библиотеку языка Си
using namespace std; // подключение стандартного пространства имён
float y_an(float x) // объявление массива чисел типа float (числа с плавающей точкой )
{
    return -sin(x);
}
float p(float x)
{
    return 0;
}
float q(float x)
{
    return 1;
}
float f1(float x)
{
    return -cos(-x);
}
int main()
{
    setlocale(LC_ALL, ""); // установка русской локали
    float x[1000],y[1000],A[1000],B[1000],C[1000],F[1000],aa[1000],bb[1000];
    float h,xx;
    int n;
    int i;
    // определяем шаг сетки, по условию h = 10^-2 = 1/100
    h=0.01;
    // количество шагов сетки = pi/2 деленное на шаг сетки +1 т.к. округляем до int в меньшую сторону.
    n =(M_PI/2)/h + 1;
    // определяем значения x на сетке
     for (i=0; i < n; i++ )
        x[i]=h*i;
    for (i=0;i < n; i++)
    {
     // параметры СЛАУ
     A[i]=1/(h*h) - p(x[i])/(2*h);
     C[i]=1/(h*h) + p(x[i])/(2*h);
     B[i]=-2/(h*h)+q(x[i]);
     F[i]=f1(x[i]);
    }
    // задаем граничные параметры
     B[0] = -h-1; C[0]=1; F[0]=0;
     B[n] = 1+h; A[n]=-1;F[n]=0;
     aa[0] = -C[0]/B[0]; bb[0]= F[0]/B[0];
    // определяем коэфициенты для определенного шага сетки
    for (i=1;i < n; i++)
    {
     aa[i] = -C[i]/(A[i]*aa[i-1] + B[i]);
     bb[i] = (F[i] - A[i]*bb[i-1])/(A[i]*aa[i-1] + B[i]);
    }
     y[n] = (F[n] - bb[n-1]*A[n])/(B[n] + aa[n-1]*A[n]);
    // решаем уравнение для каждого шага сетки
    for (i=n-1;i >= 0; i--)
        y[i] = aa[i]*y[i+1] + bb[i];
    // вывод результата
    for (i=0;i < n; i++)
        cout <<  "Шаг № " << i << ": x=" << x[i] << " res=" << y_an(x[i]) << endl;
     xx=0;
    system("pause");
    return 0;
}