Non-linear differential equation

I have this equation $$y\left(\dot y^2+1\right)=m + \Lambda y^3,$$ where $$\Lambda=1.1\cdot 10^{-52}$$ (Cosmological constant). I want to get the graph of the solution of this equation (2-parametric curve $$(\Lambda, const)$$. Thus, I have to fix one of these parameters, and another one should change in interval, for example, [0,10].

Any ideas on how I can do this in Python or Matlab? Thanks in advance

• To plot the solution one first needs to obtain the solution, no way around it. So solving the ODE is the first part of the question. Once the equation is solved, the solution will be $y(t,m,\Lambda)$, with some other parameter referring to the initial conditions (Big Bang?). Then you can plot the solution, which is the second part of the question; however it is not clear what kind of plot you want, a sequence of $y(t)$ for one parameter fixed and the other taking a set of values? Is $const$ the same as $m$? – Maxim Umansky Jun 27 at 14:54