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

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

Search the Archive

Iterate initial conditions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg48775] Iterate initial conditions
  • From: "Michael Hohendorf" <se0085 at uni-essen.de>
  • Date: Wed, 16 Jun 2004 04:54:27 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

Hi,

I have a system of differential equations, wich are results in two  
functions f[x], g[x].
Now I have to iterate the initial conditons  
(Tau_11[0],Tau_22[0],Tau_33[0]) until the functions f[x] and g[x] are  
equal to other functions Sigma_m[x] and Sigma_u[x].
Can someone help me to solve my problem?

Michael

solution = NDSolve[{
Tau_11'[x] == (2 a_T[x] v'[x] η[x] - Exp[Epsilon[x] λ[x](Ï?au_11[x]  
+ Tau_22[x] + Tau_33[x])/η[x]] Ï?au_11[x] + 2 a_T[x] v'[x] η[x] λ[x] (1 -  
ξ)) 1/(a_T[x]v[x]λ[x]),

Tau_22'[x] == (2 a_T[x] η[x] (v[x] s'[x]/s[x])- Exp[Epsilon[x]  
λ[x](Ï?au_11[x] + Tau_22[x] + Tau_33[x])/η[x]] Ï?au_22[x] + 2 a_T[x] η[x]  
λ[x] (v[x] s'[x]/s[x]) (1 - ξ)) 1/(a_T[x]v[x]λ[x]),

Tau_33'[x] == (2 a_T[x] η[x] (v[x] r'[x]/r[x])- Exp[Epsilon[x]  
λ[x](Ï?au_11[x] + Tau_22[x] + Tau_33[x])/η[x]] Ï?au_33[x] + 2 a_T[x] η[x]  
λ[x] (v[x] r'[x]/r[x])(1 - ξ)) 1/(a_T[x]v[x]λ[x]),

Tau_12'[x] == (a_T[x] λ[x] v'[x] r'[x] Tau_12[x] (1 - ξ) + a_T[x] λ[x]  
Tau_12[x] (v[x] s'[x]/s[x]) (1 - ξ) - Exp[Epsilon[x] λ[x](Ï?au_11[x]  
+ Tau_22[x] + Tau_33[x])/η[x]] Ï?au_12[x]) 1/(a_T[x]v[x]λ[x]),

Tau_11[0]==a , Tau_22[0]==b , Tau_33[0]==c , Tau_12[0]==9810},
{Tau_11[x],Tau_22[x],Tau_33[x],Tau_12[x]},{x,0,0.42}]

Tau_11[x_]=Tau_11[x]/.solution[[1]]
Tau_22[x_]=Tau_11[x]/.solution[[1]]
Tau_33[x_]=Tau_11[x]/.solution[[1]]
Tau_12[x_]=Tau_11[x]/.solution[[1]]

f[x_]=Tau_11[x]-Tau_22[x]
g[x_]=Tau_22[x]-Tau_33[x]

Iterate initial conditions until

f[x]==Sigma_m[x]
g[x]==Sigma_u[x]


  • Prev by Date: Complexes, Reals, FullSimplify
  • Next by Date: Re: Iterate initial conditions
  • Previous by thread: RE: Complexes, Reals, FullSimplify
  • Next by thread: Re: Iterate initial conditions