Ask Your Question

solve ode on given interval

asked 2017-06-10 07:24:58 -0500

soking gravatar image

updated 2017-06-13 04:19:21 -0500


I would like to tell Sage to solve ode on a given interval that I will specify. My final goal is to use this method to solve ode that contains 'piecewise' function like 2*diff(u,x) + u == f(x), where

f = piecewise([((-1,0),0), ((0,1),x), ((1,2),2-x), (RealSet.unbounded_above_closed(2),0)])


edit retag flag offensive close merge delete


did you try the approach suggested in your previous post, to write $f$ as a sum of unit_step functions? for instance, here it would be unit_step(t)*unit_step(1-t)*t + unit_step(t-1)*unit_step(2-t)*(-t+2). perhaps a method to_unit_step of piecewise functions could be useful to have..

mforets gravatar imagemforets ( 2017-06-11 11:13:59 -0500 )edit

No, I didn't think about the problem that way, I was rather think of a way to make Sage solve the ode on a specified interval. But, thanks, it's still something that works (at least for this problem).

soking gravatar imagesoking ( 2017-06-13 03:39:13 -0500 )edit

1 answer

Sort by ยป oldest newest most voted

answered 2017-06-11 15:53:21 -0500

calc314 gravatar image

Here is a numerical solution using that suggested from your previous post.

var('u x')
ode = -1/2*(u)+f
edit flag offensive delete link more

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower


Asked: 2017-06-10 07:24:58 -0500

Seen: 133 times

Last updated: Jun 13