Opened 2 years ago
Last modified 2 years ago
#22849 new enhancement
Heaviside in numerical resolutions
Reported by: | mforets | Owned by: | |
---|---|---|---|
Priority: | major | Milestone: | sage-8.0 |
Component: | symbolics | Keywords: | heaviside, integrate |
Cc: | rws, kcrisman, egourgoulhon, tscrim | Merged in: | |
Authors: | Reviewers: | ||
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description (last modified by )
- numerical ode:
sage: y = var('y') sage: desolve_rk4(x,y,ics=[0,0],end_points=1,step=0.5) # solution is x^2/2 [[0, 0], [0.5, 0.125], [1.0, 0.4999999999999999]] sage: desolve_rk4(x*unit_step(1+x),y,ics=[0,0],end_points=1,step=0.5) # OK [[0, 0], [0.5, 0.125], [1.0, 0.4999999999999999]] sage: desolve_rk4(x*heaviside(1+x),y,ics=[0,0],end_points=1,step=0.5) # ?? [[0, 0]]
- integration:
sage: integral(unit_step(x), (x, 0, 1)) 1 sage: integral(heaviside(x), (x, 0, 1)) # ?? integrate(heaviside(x), x, 0, 1)
- custom numerical value at 0:
sage: heaviside(0) heaviside(0) sage: f(x) = heaviside(x, 1/2); f(0) # new 2nd argument (?) 1/2
Change History (5)
comment:1 Changed 2 years ago by
- Description modified (diff)
comment:2 Changed 2 years ago by
- Cc egourgoulhon tscrim added
comment:3 Changed 2 years ago by
comment:4 Changed 2 years ago by
I cannot even find online documentation on hstep
, nor a mention of the Heaviside function. I don't think we can rely on Maxima having implemented it, so we cannot provide DE or integral services with it.
comment:5 Changed 2 years ago by
This is now #22850.
Note: See
TracTickets for help on using
tickets.
This is due to Maxima because with Maxima 5.38.1: