Parameters
double y0
Initial value
double start
Start Time
double end
End Time
int N
Size of output array(the larger, the finer)
Func<double, double, double> f
ode model
Return
Double[]
approximation with size N