Re: System of differential-algebraic equations

*To*: mathgroup at smc.vnet.net*Subject*: [mg78299] Re: System of differential-algebraic equations*From*: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>*Date*: Thu, 28 Jun 2007 04:15:56 -0400 (EDT)*Organization*: Uni Leipzig*References*: <f5tbnq$24q$1@smc.vnet.net>*Reply-to*: kuska at informatik.uni-leipzig.de

Hi, if you don't have initial conditions that satisfy the algebraic constrains, than there is *no* solution and you are done -- no Mathematica, no NDSolve[] You can try FindInstance[] to find valid initial conditions but the input is than random ... You can try to set boundary conditions instead of initial conditions but this will change the character of the problem. Regards Jens Jos=E9 Luis G=F3mez 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 thi= s > 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 va= lues. > We fool around a bit in the documentation, play a little bit with > AccuracyGoal, and PrecisionGoal, and with different methods specified b= y > Method, but we were not able to obtain an answer. > > > > Does anyone have a suggestion for us? Can we avoid the use of initial v= alues > 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}] >