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

MathGroup Archive 2002

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

Search the Archive

DSolve solution validation

  • To: mathgroup at smc.vnet.net
  • Subject: [mg33910] DSolve solution validation
  • From: Vladimir Bondarenko <vvb at mail.strace.net>
  • Date: Tue, 23 Apr 2002 07:13:19 -0400 (EDT)
  • Reply-to: Vladimir Bondarenko <vvb at mail.strace.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Sometimes, it is easy to double check the DSolve's output.

However, an attempt to solve

                             DSolve[y'[z] == Abs[y[z]], y[z], z]
yields

{{y[z] -> InverseFunction[(-(Log[2*Sqrt[Im[K$3541]^2 + Re[K$3541]^2] +
(2*(-Im[K$3541]^2 + Im[K$3541]*Im[#1] - Re[K$3541]^2 + Re[K$3541]*Re[#1]))/
Sqrt[Im[K$3541]^2 - 2*Im[K$3541]*Im[#1] + Im[#1]^2 + Re[K$3541]^2 -
2*Re[K$3541]*Re[#1] + Re[#1]^2]]/Sqrt[Im[K$3541]^2 - 2*Im[K$3541]*Im[#1] +
Im[#1]^2 + Re[K$3541]^2 - 2*Re[K$3541]*Re[#1] + Re[#1]^2]) + Log[2*Sqrt[Im[#1]^2 +
Re[#1]^2] + (2*((-Im[K$3541])*Im[#1] + Im[#1]^2 - Re[K$3541]*Re[#1] + Re[#1]^2))/
Sqrt[Im[K$3541]^2 - 2*Im[K$3541]*Im[#1] + Im[#1]^2 + Re[K$3541]^2 -
2*Re[K$3541]*Re[#1] + Re[#1]^2]]/Sqrt[Im[K$3541]^2 - 2*Im[K$3541]*Im[#1] +Im[#1]^2 +
Re[K$3541]^2 - 2*Re[K$3541]*Re[#1] + Re[#1]^2])*(-K$3541 + #1) & ][z + C[1]]}}

If I try to use D[] to check it, I get an expression with ByteCount of
737608. The attempt to Simplify this huge expression for the derivative
gave no answer after 1 hour at CPU = Athlon FX 1600+ / RAM = 512 Mb.

Question # 1:   Is the above shown solution correct? (It involves 2 free constants,
                K$3541 and C[1], but after simplification K$3541 might disappear?
                (Also, this ODE is not linear, so the answer to it might have 2 free
                constants?)

By hand I found

                z Sign[y[z]] - Log[y[z]] = C[1]

Question # 2:   Is this implicit solution correct? (I am not 100% sure)


Question # 3:   Looks like I've calculated a collection of ODEs about which
                I suspect that the corresponding DSolve's solutions are invalid.
                What are the possible ways to double check these solutions?
                Any modules written in Mathematica? Anything else?


Vladimir Bondarenko




  • Prev by Date: Pattern matching
  • Next by Date: Re: Row vs. Column Vectors (or Matrices)
  • Previous by thread: Re: Re: Pattern matching
  • Next by thread: Re: DSolve solution validation