       NDSolve

• To: mathgroup at smc.vnet.net
• Subject: [mg15118] NDSolve
• From: Paula Chammas <chammas at dbag.bln.daimlerbenz.com>
• Date: Sat, 12 Dec 1998 03:59:16 -0500
• Sender: owner-wri-mathgroup at wolfram.com

```Dear all,

I am having a problem with NDSolve.

I am trying to solve the following problem:

Tin[t_] := - 2 t + 12   /;  t >= 5 && t  <= 6 Tin[t_] :=  2    /;  t >=
0 && t  <= 5 Tin[t_] := 0   /;  t >= 6 && t  <= 10

Fn[t_] :=  0.785 t  /; t >= 0 && t  <= 2 Fn[t_] :=  1.57 t  /; t >= 2 &&
t  <= 5 Fn[t_] :=  0.785 t  + 5.495 /; t >= 5 && t  <= 7 Fn[t_] :=  0
/; t >= 7 && t  <= 10

If [ (w[t] ! = v[t]) && ( ( (5/6) Tin[t] - (9/6) w[t]  ) > Fn[t] ),
NDSolve[
{ w'[t] + 2 w[t] - Tin[t] + Sign[ w[t] - v[t] ] (2/3)
Fn[t] == 0,
5 v'[t]  + v[t] -    Sign[ w[t] - v[t] ] (2/3)
Fn[t] == 0,
w == 0, v ==0, {w, v}, {t, 0, 10}
],
NDSolve [
{ 6 w'[t] -Tin[t] - 3w[t] == 0, w ==0, w, {t, 0,
10 }
]

There is a discontinuity at about t =4.08--that is when w[t] = v[t] and
when the switching  should happen. I am having problems imposing the
constraints!

Any ideas?

Thanks
Cheers
Paula

--------------------------------------
Paula Chammas
DaimlerChrysler
Research and Technology
Alt-Moabit 96A           Tel: +49 30 39982 219 10559 Berlin
fax: +49 30 39982 107

paula.chammas at daimlerchrysler.com
_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/ Worry is like a rocking
chair - it gives you something to do but it doesn't get you anywhere.
-Dorothy Galyean
---------------------------------------

```

• Prev by Date: Re: help!!
• Next by Date: Non-linear ODE
• Previous by thread: Re: rician random number
• Next by thread: Non-linear ODE