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 function
Return
Double[]
approximations