Нужна помощь в проверки правельности алгоритма

class RUNGE_KUTTA {
public:
    int N;
    double t;
    double *K;
    double *KK,*K1,*K2,*K3,*K4;
    RUNGE_KUTTA(int Na) {
        N=Na;
        if(N<1) {
            return;
        }
        K = new double[N];
        KK = new double[N];
        K1 = new double[N];
        K2 = new double[N];
        K3 = new double[N];
        K4 = new double[N];
    }
    ~RUNGE_KUTTA() {
        delete[] K;
        delete[] KK;
        delete[] K1;
        delete[] K2;
        delete[] K3;
        delete[] K4;
    }
    void SET_TERM(double t0, double *K0) {
        t = t0;
        for(int i = 0; i < N; i++) {
            K[i] = K0[i];
        }
    }
    double GET_TIME() {
        return t;
    }
    void F(double t, double *K, double *dY);
    void STEP_IN(double dt) {
        if(dt < 0) {
            return;
        }
        F(t, K, K1);
        for(int i=0; i < N; i++) {
            KK[i] = K[i] + K1[i] * (dt/2.0);
        }
        F(t + dt/2.0, KK, K2);
        for(int i=0; i<N; i++) {
            KK[i] = K[i] + K2[i] * (dt/2.0);
        }
        F(t + dt/ 2.0, KK, K3);
        for(int i=0; i<N; i++) {
            KK[i] = K[i] + K3[i] * (dt);
        }
        F(t + dt, KK, K4);
        for(int i=0; i<N; i++) {
            K[i] = K[i] + dt / 6.0 * (K1[i] + 2.0 * K2[i] + 2.0 * K3[i] + K4[i]);
        }
        t = t + dt;
    }
};

void F(double t, double *K, double dY[]) { //y''(t) + y(t) = 0 (y'' = 2-ая производная.
        dY[0] = K[1];
        dY[1] = -K[0];
    }

int main() {
    RUNGE_KUTTA rk4(2);
    double K0[2] = {0, 1};
    rk4.SET_TERM(0, K0);
    for(int i = 0; i < 10; i++) {
        cout << rk4.GET_TIME() << " " << rk4.K[0] /*Искомое решение*/ << " " << rk4.K[1] /* i - ая производная*/ << endl;
        rk4.STEP_IN(0.01);
    }
    return 0;
}