numerical method to find solution of differential equation