Re: Differential eqn. with interpolated coeficients
- To: mathgroup at smc.vnet.net
- Subject: [mg29187] Re: [mg29158] Differential eqn. with interpolated coeficients
- From: Ben Jacobson <bjacobson at illumitech.com>
- Date: Sat, 2 Jun 2001 05:54:41 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Works fine for me. I would guess you have a syntax error.
In[2]:=
p[x_] = Interpolation[Table[{j, Sin[j]}, {j, 0, 6, .6}]][x]
Out[2]=
InterpolatingFunction[{{0., 6.}}, "<>"][x]
In[6]:=
q[x_] = Interpolation[Table[{j, j*Cos[j]}, {j, 0, 6, .6}]][x]
Out[6]=
InterpolatingFunction[{{0., 6.}}, "<>"][x]
In[9]:=
NDSolve[{y''[x] + p[x]y'[x] + q[x] == 0, y[0] == .5, y'[0] == .3}, y, {x, 0,
6}]
Out[9]=
{{y -> InterpolatingFunction[{{0., 6.}}, "<>"]}}
Ben Jacobson
Illumitech Inc.
At 04:15 AM 6/1/01 -0400, you wrote:
>I need to solve a ordinary second order differential equation, of the
>type y''[x] + P[x] y'[x] + Q[x] == 0, where the coefficients P[x] and
>Q[x] are an interpolating function (or a table of numerical data).
>
>I have tried with NDSolve, but it seems to admit only an interpolating
>function for Q[x], but not for P[x]. Does any one know what could I do?
>
>Regards,
>
>Daniel Arteaga
>
>