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