Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2009

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

Search the Archive

Re: Solving nonlinear coupled differential equations in Mathematica

  • To: mathgroup at smc.vnet.net
  • Subject: [mg95297] Re: Solving nonlinear coupled differential equations in Mathematica
  • From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
  • Date: Wed, 14 Jan 2009 05:55:49 -0500 (EST)
  • References: <gkgpvp$48j$1@smc.vnet.net>

Hi,

your equations are not a differential equation for
z[x] because

ode = {y''[x] == y[x] + z'[x],
   z[x]^3 == y'[x] + 3};

ode /. Solve[D[#, x] & /@
     Last[ode], z'[x]][[1]]

gives:

{Derivative[2][y][x] == y[x] + Derivative[2][y][x]/(3*z[x]^2),
  z[x]^3 == 3 + Derivative[1][y][x]}

that does not include z'[x] any more and more over
z[x] is complete undetermined.

Setting

ode1 = ode /. Solve[D[#, x] & /@
      Last[ode], z'[x]][[1]];

yodes= First /@ (ode1 /. Solve[Last[ode1], z[x]])

gives three possible equations for y[x] and

NDSolve[
   {#, y[1] == 2 , y[0] == 3},
   y[x], {x, 0, 1}] & /@ yodes

will solve it with some numerical error messages
because I'm not able to to take the right solution
for z[x]

Regards
   Jens

SK wrote:
> Here is a set of equations I would like to solve
> 
> y''[x]==y+z'[x]
> z[x]^3=y'[x]+3
> 
> With boundary conditions of y[1]==2 and y[0]==3
> 
> Using NDSolve on these equations, Mathematica says the order of the
> equations is 3 and it has only 2 initial conditions. But the order of
> this system of equations is 2 as far as I see (since order is defined
> as the highest derivative)
> 
> When I do try to put another boundary condition in like z[0]==0
> Mathematica spits out that it cant solve for the derivatives and is
> using a mass matrix method (error:ntdvmm) and then it says that it has
> significant errors (error:berr) and will return the best solution
> found.
> 
> Any help on this matter will be greatly appreciated.
> Thanks
> 


  • Prev by Date: Re: Problems using MathLink
  • Next by Date: Re: File Corrupted
  • Previous by thread: Re: Problems using MathLink
  • Next by thread: Re: Solving nonlinear coupled differential equations in Mathematica