MathGroup Archive 2004

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

Search the Archive

Question on PDE

  • To: mathgroup at
  • Subject: [mg48696] Question on PDE
  • From: Lautaro Vergara <lvergara at>
  • Date: Thu, 10 Jun 2004 02:44:50 -0400 (EDT)
  • Sender: owner-wri-mathgroup at

Dear people,

I need to automatize a problem like the one shown below. I'll appreciate
very much any comments to this respect.

Let the following PDE that depends on two parameters, a and b

\!\(s1 := NDSolve[{
    2\ Pi\ \((t^2 +
      V1[x, t])\)\ D[U[x, t], t] + V2[x, t,a] \[Equal] 0, U[x,
                t0] \[Equal] \((\(-x\) + 4\
                lambda\ x^3)\)\ \[ExponentialE]\^\(\(-a\)\ x\^2\),
        U[0, t] \[Equal] 0,
           U[b, t] \[Equal] 0}, U, {x, 0, b}, {t, t0, tf}]\)



V1[x_,t_]:=Exp[a x^2] (2 a x U[x,t]+D[U[x,t],x])

\!\(V2[x_, t_] := \((2\ a\ \((1 + 2\
      a\ x\^2)\)\ U[x, t] + 4\ a\ x\ D[U[x, t], x] + D[D[U[x, t], x],

I need to evaluate the following function and found the values of a and
b that makes this function independent of them

s2a[x_,t_]:=Exp[a x^2] U[x,t]/.s1[[1]]

I have done by hand by making a DO of the variable a and then playing
with different values of b such that what is looked for happens.

Do you know a way of doing this more efficiently?

Thanks in advance,


  • Prev by Date: AW: Combining plots
  • Next by Date: Re: using iFFT on a Continuous Time Transfer Function
  • Previous by thread: AW: Combining plots
  • Next by thread: (N)Integrate help - convergence problem