Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2006
*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 2006

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

Search the Archive

Initial conditions for DAE

  • To: mathgroup at smc.vnet.net
  • Subject: [mg65788] Initial conditions for DAE
  • From: Maarten <maarten at stack.nl>
  • Date: Mon, 17 Apr 2006 02:28:45 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

Dear fellow users,

I would like to use NDSolve to solve a DAE with a single algebraic variable. The algebraic variable has to be solved numerically, so SolveDelayed->True. The problem is that it has a root of multiplicity >1. When starting the integration I know the correct value for the algebraic variable and I would like to use this value (eg. 0.9) as a starting point for the numerical root finding that occurs during initialisation. However, when I specify an additional equality of the form:

AEvar[0]==0.9

in addition to the algebraic equation

f(AEvar[t])==0

NDSolve reports the initial conditions are inconsistent (they are overdetermined probably are initial conditions for all the differential variables are also defined). As far as I can see there are no options to pass NDSolve my initial guess, while this is a common feature in other DAE solving programs I have worked with? Is there a way of accomplishing this without resorting to FindRoot functions in my DAE problem formulation?

Best Regards,

Maarten Nauta


  • Prev by Date: Re: Problem with limiits
  • Next by Date: Re: importing a mathematica plot in a latex file
  • Previous by thread: Re: Define new constant as Pi is
  • Next by thread: No success with a package from the Wolfram Website