#6479: [with patch, needs work] desolve for 2nd order ODE with initial condition
gives wrong answer
-------------------------+--------------------------------------------------
Reporter: gmhossain | Owner: burcin
Type: defect | Status: needs_work
Priority: major | Milestone: sage-4.2.1
Component: calculus | Keywords:
Work_issues: | Author:
Reviewer: | Merged:
-------------------------+--------------------------------------------------
Comment(by wdj):
Replying to [comment:11 robert.marik]:
> Replying to [comment:10 wdj]:
...
> > A request: in your new functions desolve_rk4 and desolve_system_rk4
> > there is an option called endpoint, with default value 10. I would
> > prefer an option called endpoints with a default value of [0,10]
> > (or something), so that a range can be plotted other than from
> > ics[0] to endpoint. If it is too much hassle, fine (you can just add
> > plots together to get that anyway...).
>
> what about this:
>
> endpoints=a .... integrate from ics[0] to a
>
> endpoints=[a] .... integrate from ics[0] to a
>
> endpoints=[a,b] .... integrate from ics[0] to b, then integrate back
from ics[0] to a, reverse the second list and join both lists together
/without repeating the point (ics[0],ics[1])/. If ics[0] is bigger than b
or smaller than a, raise error.
>
> I think that this is possible and I can try within a week.
This sounds excellent - thanks!
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/6479#comment:12>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
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-trac?hl=en
-~----------~----~----~----~------~----~------~--~---
