| Author |
Comment/Response |
dani
|
 06/21/03 1:29pm
In[23]:=
NDSolve[{
a'[z] == -0.1*a[z] - 149*10^6*2*a[z]*b[z]*Sin[d[z] + 10^5*z],
b'[z] == -0.1*b[z] - 149*10^6*2*a[z]^2*Sin[2*c[z] - d[z] - 10^5*z],
c'[z] == -149*10^6*b[z]*Cos[d[z] + 10^5*z],
d'[z] == -149*10^6*2*a[z]^2/b[z]*Cos[2*c[z] - d[z] - 10^5*z],
a[0] == 0.9, b[0] == 0.9, c[0] == 0.5, d[0] == 0.5}, {a, b, c, d}, {z,
10^-10, 5*10^-5}]
{{a -> InterpolatingFunction[{{0., 00005.}}, <>],
b -> InterpolatingFunction[{{0., 00005.}}, <>],
c -> InterpolatingFunction[{{0., 00005.}}, <>],
d -> InterpolatingFunction[{{0., 00005.}}, <>]}}
Plot[Evaluate[{a[z], b[z], c[z], d[z]} /. %], {z, 10^-10, 10^-5}]
From In[23]:=
NDSolve::"mxst": "Maximum number of \!\(1000\) steps reached at the point \
\!\(z\) == \!\(2.129961290733448`*^-8\)."
Out[23]=
\!\(\*
RowBox[{"{",
RowBox[{"{",
RowBox[{
RowBox[{"a", "->",
TagBox[\(InterpolatingFunction[{{1.`*^-10,
2.129961290733448`*^-8}}, "<>"]\),
False,
Editable->False]}], ",",
RowBox[{"b", "->",
TagBox[\(InterpolatingFunction[{{1.`*^-10,
2.129961290733448`*^-8}}, "<>"]\),
False,
Editable->False]}], ",",
RowBox[{"c", "->",
TagBox[\(InterpolatingFunction[{{1.`*^-10,
2.129961290733448`*^-8}}, "<>"]\),
False,
Editable->False]}], ",",
RowBox[{"d", "->",
TagBox[\(InterpolatingFunction[{{1.`*^-10,
2.129961290733448`*^-8}}, "<>"]\),
False,
Editable->False]}]}], "}"}], "}"}]\)
From In[23]:=
Syntax::"sntxf": "\"\!\(a ->\)\" cannot be followed by \
\"\!\(\(InterpolatingFunction[\(\(\(\({\({0., 00005.}\)}\)\), \
<>\)\)]\)\)\".\!\(\"\"\)"
From In[23]:=
{{a -> InterpolatingFunction[{{0., 00005.}}, <> ],
b -> InterpolatingFunction[{{0., 00005.}}, <>],
c -> InterpolatingFunction[{{0., 00005.}}, <>],
d -> InterpolatingFunction[{{0., 00005.}}, <>]}}
URL: , |
|
