As you pointed out, the equation has separated variables, so it must
work allmost in any version
Perhaps problems with keepfloat:true. Is you replace 0.5 by 1/2,
everything works.
See below
Robert
[ma...@um-bc107 /opt/sage]$ ./sage
----------------------------------------------------------------------
| Sage Version 4.2.1.alpha0, Release Date: 2009-11-09 |
| Type notebook() for the GUI, and license() for information. |
----------------------------------------------------------------------
**********************************************************************
* *
* Warning: this is a prerelease version, and it may be unstable. *
* *
**********************************************************************
Loading Sage library. Current Mercurial branch is: inequality
sage: (g,t) = var('g,t')
sage: y=function('y',t)
sage: a = (1/2)^2*pi
sage: DEiii = pi * (39/100*y + 1/2)^2* diff(y,t) + a * sqrt(2*g*y)
sage: DEiii
1/10000*(39*y(t) + 50)^2*pi*D[0](y)(t) + 1/4*pi*sqrt(g*y(t))*sqrt(2)
sage: desolve(DEiii, [y,t])
-1/12500*(1521*y(t)^2 + 6500*y(t) + 12500)*sqrt(g*y(t))*sqrt(2)/g == c
+ t
sage:
On 11 lis, 21:13, David Joyner <[email protected]> wrote:
> On Wed, Nov 11, 2009 at 3:10 PM, Jaasiel Ornelas <[email protected]> wrote:
>
> > I'm using:
>
> > Sage Version 4.1.1, Release Date: 2009-08-14
>
> Okay. The above was done with Sage 4.2, so may or may not work in your
> version.
>
>
>
> > We just started laplace in my class, so I did not know that laplace
> > transforms don't solve non linear ode's. Sorry :p
>
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URL: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
