Skip to content

Instantly share code, notes, and snippets.

@fusion809

fusion809/desolve-sle-attempt.sage

Created July 24, 2019 13:44
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save fusion809/b3e397f714e211ea7a18dcb9724f57ce to your computer and use it in GitHub Desktop.
Save fusion809/b3e397f714e211ea7a18dcb9724f57ce to your computer and use it in GitHub Desktop.
Download ZIP
Raw
desolve-sle-attempt.sage
y=function("y")
var("lambda_")
E=diff(y(x),x,2)-x*y(x)==lambda_*y(x)
S=desolve(E,y(x),ivar=x, contrib_ode=True, ics=[0,0,+Infinity,0]);
@fusion809
Copy link
Author

Output:

---------------------------------------------------------------------------
NotImplementedError                       Traceback (most recent call last)
<ipython-input-3-773d704ea403> in <module>()
      2 var("lambda_")
      3 E=diff(y(x),x,Integer(2))-x*y(x)==lambda_*y(x)
----> 4 S=desolve(E,y(x),ivar=x, contrib_ode=True, ics=[Integer(0),Integer(0),+Infinity,Integer(0)]);

/usr/lib/python2.7/site-packages/sage/calculus/desolvers.pyc in desolve(de, dvar, ics, ivar, show_method, contrib_ode, algorithm)
    604             if not show_method:
    605                 maxima_method=P("method")
--> 606             raise NotImplementedError("Unable to use initial condition for this equation (%s)."%(str(maxima_method).strip()))
    607         if len(ics) == 2:
    608             tempic=(ivar==ics[0])._maxima_().str()

NotImplementedError: Unable to use initial condition for this equation (odelin).

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment