Неизвестная ошибка

#include <conio.h>
#include <cmath>
#include <cstdlib>
#include <iostream>
 
void func (double *y, double *ys, double t)
 
{ // функция вычисления правых частей уравнений
 
ys[0] = y[1]; // ys[1]-первая производная; ys[2]-вторая и т.д.
 
ys[1] = y[2]; // t-независимый аргумент
 
ys[2] = 5 + t * t - y[0] - 3.0 * y[1] - 2.0 * y[2];
    //*ys=2.5+0.5*t*t-0.5*y[1]-0.5*y[1]*y[1]*y[1]-1.5*y[2];
 
}
 
void runge_k (
 
void f (double *y, double *ys, double x), double *y, int n, double tn, double tk,int m, double eps) 
{
 
 double *k1,*k2,*k3,*k4,*Yh;
 
 k1 = new double[n];
 k2 = new double[n];
 k3 = new double[n];
 k4 = new double[n];
 
 Yh = new double[n+1];
 
    double h=(tk - tn)/m;
    double* ys=y+n+1;
    double X;
    for(X=tn;X<tk;X+=h){
        f(y,ys,X);
        
        for( int i=0;i<n;i++) {k1[i]=h*y[i+1];}
        for( int i=0;i<n;i++) {Yh[i]=y[i]+0.5*k1[i];}
        f(Yh,Yh+n+1,X+h*0.5);
        
        for( int i=0;i<n;i++) {k2[i]=h*Yh[i+1];}
        for( int i=0;i<n;i++) {Yh[i]=y[i]+0.5*k2[i];}
        f(Yh,Yh+n+1,X+h*0.5);
        
        for( int i=0;i<n;i++) {k3[i]=h*Yh[i+1];}
        for( int i=0;i<n;i++) {Yh[i]=y[i]+k3[i];}
        f(Yh,Yh+n+1,X+h);
        
        for( int i=0;i<n;i++) {k4[i]=h*Yh[i+1];}
        f(Yh,Yh+n+1,X+h);
        
        for( int i=0;i<n;i++) y[i]+=(k1[i]+2*k2[i]+2*k3[i]+k4[i])/6.0;
    }
 /*   h=tk-X;
    f(y,ys,X);
    
    for( int i=0;i<n;i++) k1[i]=h*y[i+1];
    for( int i=0;i<n;i++) Yh[i]=y[i]+0.5*k1[i];         
    f(Yh,Yh+n+1,X+h*0.5);
    
    for( int i=0;i<n;i++) k2[i]=h*Yh[i+1];
    for( int i=0;i<n;i++) Yh[i]=y[i]+0.5*k2[i];
    f(Yh,Yh+n+1,X+h*0.5);
    
    for( int i=0;i<n;i++) k3[i]=h*Yh[i+1];
    for( int i=0;i<n;i++) Yh[i]=y[i]+k3[i];
    f(Yh,Yh+n+1,X+h);
    
    for( int i=0;i<n;i++) k4[i]=h*Yh[i+1];
    
    for( int i=0;i<n;i++) y[i]+=1./6.*(k1[i]+2*k2[i]+2*k3[i]+k4[i]);*/
}
 
 
int _tmain(int argc, _TCHAR* argv[])                            
{
 
double y[3], xs, xe;
 
int i;
 
y[0] = 1.0; y[1] = 0.1; y[2] = 0.0; // Начальные условия
 
runge_k (func, y, 1000, 0, 2, 10, 0.001);
 
    return 0;
}