On Oct 9, 9:26 am, kcrisman <kcris...@gmail.com> wrote:
> On Oct 9, 7:46 am, "ma...@mendelu.cz" <ma...@mendelu.cz> wrote:
>
>
>
>
>
> > Hello,
>
> > trying to fix desolve_laplace as described
> > athttp://groups.google.cz/group/sage-support/browse_thread/thread/b6f6b...
>
> > It is continuation of tickethttp://trac.sagemath.org/sage_trac/ticket/6479
> > which has been (hope) solved.
>
> > The temporary code ishttp://user.mendelu.cz/marik/temp/desolvers.py
> > and it has been tested on equations like equations on the bottom of
> > this post.
>
> > I have still a problem: I do not know how to get sage representation
> > of 'at(....)
>
> > current desolve_laplace produce (via maxima) something like
> > "x*%e^x*('at('diff('f(x),x,1),x=0))-'f(0)*x*%e^x+'f(0)*%e^x"
>
> > If I do
> > A=maxima("'at('diff('f(x),x,1),x=0)")
> > A.sage()
>
> > I expect the corresponding Sage expression but I get errors.
>
> Ah, but there IS no corresponding Sage expression! That's the kicker.
> See trac #385, as well as #3914 which would also be solved by it. So
> this has been around a long time.
>
> Ticket #385is on my to-do list, and since you are interested, I can
> try to get around to it soon. Or if you have a good solution at hand,
See patch at #385. Can you test it with your new desolve_laplace and
desolve code? That would help expose any errors in it. Maybe Sage
will finally be able to do everything it should in an intro DE
course...
- kcrisman
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to
sage-devel-unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
