Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1997
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1997

[Date Index] [Thread Index] [Author Index]

Search the Archive

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



  • Prev by Date: Re: groups
  • Next by Date: Re: Creating table with multiple lists
  • Previous by thread: Re: Changing PDE to ODE
  • Next by thread: A question !!!