Re: Changing PDE to ODE
- To: mathgroup at smc.vnet.net
- Subject: [mg7524] Re: [mg7509] Changing PDE to ODE
- From: "w.meeussen" <w.meeussen.vdmcc at vandemoortele.be>
- Date: Tue, 10 Jun 1997 23:17:59 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
hi Sergio,
***** your input ****
In[32]:=eq1 = a^2*u[0][x, y, z] +
Derivative[2, 0, 0][u[0]][x, y, z] == 0
Out[32]=
2 (2,0,0)
a u[0][x, y, z] + (u[0]) [x, y, z] == 0
In[33]:=eq2 = b^2*Derivative[0,1,0][u[1]][x,y,z] +
Derivative[0,2,0][u[1]][x,y,z]==0
Out[33]=
2 (0,1,0) (0,2,0)
b (u[1]) [x, y, z] + (u[1]) [x, y, z] == 0
***** definitions : *****
In[34]:=u[0][x,y,z]:=u[0][y];
Derivative[n_,0,0][u[0]][x,y,z]:=Derivative[n][u[0]][x]
In[36]:=u[1][x,y,z]:=u[1][y];
Derivative[0,n_,0][u[1]][x,y,z]:=Derivative[n][u[1]][y]
***** end definitions *****
and, without further adoo :
In[38]:=eq1
Out[38]=
2
a u[0][y] + (u[0])''[x] == 0
In[39]:=eq2
Out[39]=
2
b (u[1])'[y] + (u[1])''[y] == 0
It's what's in The Book under "LaplaceTransform", discussing "Heat
Conduction PDE".
greetings,
At 03:48 7-06-97 -0400, Sergio Rojas wrote:
> Hello guys,
>
> I have a set of PDE that may be solved using ODE techniques.
> In general, they look (after using FullForm) like :
>
>eq1 = a^2*u[0][x, y, z] + Derivative[2, 0, 0][u[0]][x, y, z] == 0
>
>eq2 = b^2*Derivative[0,1,0][u[1]][x,y,z] + Derivative[0,2,0][u[1]][x,y,z]==0
>
> By using Mathematica, how these equations can be
> transformed into ordinary differential equation
> that may looks like:
>
>a^2*u[0][x] + Derivative[2][u[0]][x] == 0
>
>b^2*Derivative[1][u[1]][y] + Derivative[2][u[1]][y] == 0
>
> One can take care of the Derivative by using:
>
>neweq2 = eq2 /. Derivative[0, n_, 0] -> Derivative[n]
>
> How to change [x,y,z] to [y] ?
>
>Sergio
>
>E-mail: sergio at scisun.sci.ccny.cuny.edu
>
>
>
>
>
Dr. Wouter L. J. MEEUSSEN
eu000949 at pophost.eunet.be
w.meeussen.vdmcc at vandemoortele.be