Re: ODE with in a piece wise plot

• To: mathgroup at smc.vnet.net
• Subject: [mg104243] Re: ODE with in a piece wise plot
• From: Hannes Kessler <HannesKessler at hushmail.com>
• Date: Sun, 25 Oct 2009 01:04:43 -0400 (EDT)
• References: <hbu7ps\$7he\$1@smc.vnet.net>

```Hello Becky,

I don't understand what you mean with "put in a ODE where the (x-1)
is".
Nevertheless, here is an example how one can use your function in an
ODE:

E1 = If[-1 <= x <= 0, 1,   If[3 <= x <= 5, 2, If[-1 <= x <= 3=
, (x -
1)]]]

(* Equivalent but better looking *)
E1=PiecewiseExpand[E1]

(* Function definition using E1*)
f[xx_] := E1 /. x :> xx

Plot[f[x], {x, -5, 5}]

(* Use it as rhs in an ODE *)
equation=y'[x] + y[x] == f[x]
solution = NDSolve[{equation, y[-5] == 1}, y, {x, -5, 5}]

(* The solution function and its plot *)
yfunc = y /. First[solution]
Plot[yfunc[x], {x, -5, 5}]

By the way, you can enter a piecewise function definition similar to
the one obtained by PiecewiseExpand by entering <escape>pw<esacpe>
followed by one or more <ctrl>enter. Or use Piecewise[..] (look into
the Mathematica help for Piecewise).

Best regards,
Hannes Kessler

On 24 Okt., 08:45, "Becky" <noslow... at comcast.net> wrote:
> Can this be done?
>
> First of all, a basic example and then what I want to change.
>
> E1=If[-1<=x<=0,1,If[3<=x<=5,2,If[-1<=x<=3,(x-1)]]]
> Plot[E1,{x,-5,5}]
>
> The purpose of this was to join all three areas into one plot.
>
> Now!  I want to put in a ODE where the (x-1) is.
>
> Is this possible, and if so, would you mind giving me an example?

```

• Prev by Date: Re: Re: which values of m satisfies the inequality
• Next by Date: Re: How to print/display symbolic expressions from inside a Module
• Previous by thread: ODE with in a piece wise plot
• Next by thread: Re: ODE with in a piece wise plot