I am trying to solve the differential equation using solve_ivp, but I am not getting the right solution. However, I obtained the right solution using