solution of a ordinary differential equation using matlab