I Use the function NDSolve of mathematica to solve the …