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

MathGroup Archive 2007

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

Search the Archive

RV: System of differential-algebraic equations

  • To: mathgroup at smc.vnet.net
  • Subject: [mg78372] RV: [mg78271] System of differential-algebraic equations
  • From: José Luis Gómez <jose.luis.gomez at itesm.mx>
  • Date: Fri, 29 Jun 2007 05:49:15 -0400 (EDT)

You are very kind to use your time for answering our question, even it was a
bad question (sorry for that!) I will check your suggestions. Thank you!

Jose


-----Mensaje original-----
De: DrMajorBob [mailto:drmajorbob at bigfoot.com]
Enviado el: Mi=E9rcoles, 27 de Junio de 2007 02:25 p.m.
Para: Jos=E9 Luis G=F3mez
Asunto: Re: [mg78271] System of differential-algebraic equations

Attached is a partial solution for the initial conditions.

I took only the first solution at one point, and it may not be a 
MEANINGFUL solution. You should look at that further with the underlying 
problem in mind.

After evaluating subFunction[m]...

m[t] and Subscript[m,t] are interchangeable, they both evaluate to m[t], 
and they both display in the subscripted form. I often find that 
convenient.

You can take out the first cell without affecting the rest of the 
notebook, except for the way the functions are displayed.

Bobby

On Wed, 27 Jun 2007 04:29:43 -0500, Jos=E9 Luis G=F3mez 
<jose.luis.gomez at itesm.mx> wrote:

> Dear Mathematica Group.
>
>
> A colleague has asked me help to solve a system of 8 algebraic and
> differential equations. The system is included below, at the end of this
> e-mil, in InputForm.
>
> Mathematica 6.0 NDSolve command replies with this message:
>
>
> NDSolve::icfail: Unable to find initial conditions which satisfy the
> residual function within specified tolerances.  Try giving initial
> conditions for both values and derivatives of the functions.
>
>
> Now, my colleague does not want to give initial conditions for the
> derivatives, because he does not have actual information about those 
> values.
> We fool around a bit in the documentation, play a little bit with
> AccuracyGoal, and PrecisionGoal, and with different methods specified by
> Method, but we were not able to obtain an answer.
>
>
> Does anyone have a suggestion for us? Can we avoid the use of initial  
> values
> for the derivatives?
>
>
> The system is included below. Thanks in advance for any advice.
>
>
> Jose Luis Gomez-Munoz
>
>
>
>
> NDSolve[{m[t]*x[t] + q[t]*x[t]^2 == 2.75,
>
>    (12.6/10^15)*m[t]^2 - (10.2/10^16)*m[t]*r[t] -
>
>      (20.4/10^16)*m[t]*u[t]*y[t] + (20.4/10^16)*m[t]*u[t]*x[t] -
>
>      (92.4084/10^15)*q[t] == 0,
>
>    r[t]*(y[t] - x[t]) + u[t]*(y[t] - x[t])^2 == 0.444,
>
>    (10.2/10^16)*r[t]^2 - (9.82/10^12)*r[t]*v[t] -
>
>      (389.9256/10^16)*u[t] == 0, v[t] + 2*w[t]*(y[t] - 140/10^16) ==
>
>     0, 0.5*Derivative[1][m][t]*x[t]^2 +
>
>      m[t]*x[t]*Derivative[1][x][t] + (1/3)*Derivative[1][q][t]*
>
>       x[t]^3 + q[t]*x[t]^2*Derivative[1][x][t] -
>
>      (25.2/10^15)*q[t]*x[t] == 0,
>
>    0.444*Derivative[1][x][t]*0.5*r[t]*(y[t] - x[t])^2 +
>
>      (1/3)*Derivative[1][u][t]*(y[t] - x[t])^3 +
>
>      (Derivative[1][y][t] - Derivative[1][x][t])*
>
>       (r[t]*(y[t] - x[t]) + u[t]*(y[t] - x[t])^2) -
>
>      (20.4/10^16)*u[t]*(y[t] - x[t]) == 0,
>
>    0.5*Derivative[1][v][t]*(y[t] - 140/10^16)^2 +
>
>      v[t]*(y[t] - 140/10^6)*Derivative[1][y][t] +
>
>      (1/3)*Derivative[1][w][t]*(y[t] - 140/10^6)^3 +
>
>      (19.64/10^12)*w[t]*(140/10^6 - y[t]) + w[t]*(y[t] - 140/10^6)^2*
>
>       y[t] == 0, m[0] == 0., q[0] == 2.342*^8, r[0] == 3.7843*^7,
>
>    u[0] == -1.385*^13, v[0] == 0., w[0] == 1.9856*^7,
>
>    x[0] == 1.7/10^6, y[0] == 4.42/10^6}, {m, q, r, u, v, w, x, y},
>
>   {t, 0, 3600}]
>
>



--
DrMajorBob at bigfoot.com



  • Prev by Date: Re: Suggestions for Maintaining "Object" State?
  • Next by Date: Re: [Mathematica 6] Integrate strange result
  • Previous by thread: a definite integral and a workaround
  • Next by thread: Iterpolation and Dates