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