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