Метод стрельбы: ошибки

#include <cstdio>
#include <math.h>
#include <iostream>
 
using namespace std;
 
void Show(int);
 
void Show(int n) {
    cout << "Исходная задача:\n";
    cout << "x'(t) = 0.5*x(t) + t*t*y(t)+t*t-1";
    cout << "y'(t) = -0.5*y(t) + x(t)\n";
    cout << "При условиях:\n";
    cout << "х(0) = -1\n";
    cout << "х(1.5) = 0\n";
    cout << endl;
}
 
double f(double x, double y, double t) // x'
{
    return (0.5*x + t * t*y + t * t - 1);
}
 
double g(double x, double y, double t) // y'
{
    return -0.5 * y + x; 
}
 
void rk4(double a, double b, double h, double y0, double yb, double Y[])
{
    double k1[4]; //коэффиценты К
    double k2[4]; //коэффиценты К
    int i = 0;
    Y[0, 0] = y0;  
    Y[1, 0] = yb;
    for (double x = a; x < b - h; x += h)
    {
        k1[0] = f(x, Y[0, i], Y[1, i]);
        k2[0] = g(x, Y[0, i], Y[1, i]);
 
        k1[1] = f(x + h / 3, Y[0, i] + k1[0] * (h / 3), Y[1, i] + k2[0] * (h / 3));
        k2[1] = g(x + h / 3, Y[0, i] + k1[0] * (h / 3), Y[1, i] + k2[0] * (h / 3));
 
        k1[2] = f(x + h * 2 / 3, Y[0, i] - k1[1] * (h / 3), Y[1, i] - k2[1] * (h / 3));
        k2[2] = g(x + h * 2 / 3, Y[0, i] - k1[1] * (h / 3), Y[1, i] - k2[1] * (h / 3));
 
        k1[3] = f(x + h, Y[0, i] + k1[0] - k1[1] + k1[2], Y[1, i] + k1[0] - k1[1] + k1[2]);
        k2[3] = g(x + h, Y[0, i] + k1[0] - k1[1] + k1[2], Y[1, i] + k1[0] - k1[1] + k1[2]);
 
        Y[0, i + 1] = Y[0, i] + (h / 8) * (k1[0] + 3 * k1[1] + 3 * k1[2] + k1[3]);
        Y[1, i + 1] = Y[1, i] + (h / 8) * (k2[0] + 3 * k2[1] + 3 * k2[2] + k2[3]);
        i++;
    }
}
 
double shooting_method(double a, double b, double y0, double y1, double h, double eps)
{
    double nu1 = 1.0, nu2 = 0.8, nu3; // некоторое значение тангенса угла наклона касательной к решению в точке x = a
    int N = (int)((b - a) / h + 1); // количество элементов
    double[, ] Y1 = new double[2, N], Y2 = new double[2, N]; // массивы для хранения решений
    Y1 = rk4(a, b, h, y0, nu1, Y1); // решаем систему для nu1
    Y2 = rk4(a, b, h, y0, nu2, Y2); // решаем систему для nu2
    while (abs(Y2[0, N - 1] - y1) > eps) //
    {
        nu3 = nu2 - ((nu2 - nu1) / (Y2[0, N - 1] - Y1[0, N - 1])) * (Y2[0, N - 1] - y1); // вычисляем новое nu
        nu1 = nu2;
        nu2 = nu3;
        Y1 = Y2;
        Y2 = new double[2, N]; // без этого не работает ?!? о_О
        Y2 = rk4(a, b, h, y0, nu2, Y2); //вычисляем решение для нового nu
    }
    return Y2;
}
 
void main() {
    setlocale(0, "");
    int a = 0, b = 1.5, // отрезок [a; b] 
        n = 10;
    double h = (b - a) / n; // шаг
    double eps = 0.001; // точность
    Show(n);
    double x0 = -1; //y(a)
    double x1 = 0; //y(b)
    int N = (int)((b - a) / h + 1);
    double[, ] Y = shooting_method(a, b, x0, x1, h, eps);
    
    
}