| 1 | initial version |
The following code solves your equation numerically (with 4th order Runge-Kutta method) for y(0)=0.5 and plots the solution:
x,y=var('x,y')
mynumericalsolution=desolve_rk4(tan(y) * (y^3 - y), dvar=y, ivar=x, ics=[0,0.5],end_points=10, step=0.1)
myplot=line(mynumericalsolution,axes_labels=[x,y])
myplot.show()
I will edit if I can find an exact solution with SageMath.
| 2 | No.2 Revision |
The following code solves your equation numerically (with 4th order Runge-Kutta method) for y(0)=0.5 and plots the solution:
x,y=var('x,y')
mynumericalsolution=desolve_rk4(tan(y) * (y^3 - y), dvar=y, ivar=x, ics=[0,0.5],end_points=10, step=0.1)
myplot=line(mynumericalsolution,axes_labels=[x,y])
myplot.show()
I will edit if I can find an exact solution with SageMath.think SageMath cannot solve it symbolically.
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.