x********g
发帖数: 595 | 1 When using NDSolve to solve partial differential equations, get the error
message:
"
NDSolve::deqn: Equation or list of equations expected instead of True in the
first argument
" |
|
l******g
发帖数: 3 | 2 在用NDSolve解DE的时候, 那本厚厚的手册里给出一些参数选项, 比如: AccuracyGoal,
PrecisionGoal, MaxSteps, StartingStepSize, MaxStepSize.
我在用NDSolve求解一个Time Dependent Defferential Equation, 需要对MaxSteps和
StartingStepSize选择合适的值才能给出正确的结果. 但是我不理解MaxSteps,
StartingStepSize的含义和他们之间的关系. 因此总是做不对. 恳请各位高手帮我,
不胜感激 |
|
c****l
发帖数: 88 | 3 Hi,
I am using Mathematica to numercially solve the following equations:
2 D[f0[r], r]/r + D[f0[r], r, r] == -2 A f1[r]/r^4 + 2 A D[f1[r], r]/r^3
-2 f1[r]/r^2 + 2 D[f1[r], r]/r + D[f1[r], r, r] == 2 A D[f0[r], r]/r^3
And my boundary condition is
f0[10] == 0.01, f1[10] == 0.01, (D[f0[r], r] /. r -> 1) == 10^-3 A,
(D[f1[r], r] /. r -> 1) == 0
When I take A to be small, say 1, everything is fine. However, if I take A
to be large, e.g. 100. Mathematica complains and gives me crazy results.
Here ... 阅读全帖 |
|
h***o
发帖数: 539 | 4 BBS水木清华站∶精华区
发信人: FangQ (奥萨马·本·拉登), 信区: MathTools
标 题: Mathematica函数及使用方法
发信站: BBS 水木清华站 (Fri Nov 20 21:39:10 1998)
Mathematica函数及使用方法
—————————————————————————————————————
八、数值函数
N[expr] 表达式的机器精度近似值
N[expr, n] 表达式的n位近似值，n为任意正整数
NSolve[lhs==rhs, var] 求方程数值解
NSolve[eqn, var, n] 求方程数值解，结果精度到n位
NDSolve[eqns, y, {x, xmin, xmax}]微分方程数值解
NDSolve[eqns, {y1,y2,...}, {x, xmin, xmax}]
|
|
x********g
发帖数: 595 | 5 如图：
用NDSolve解完之后，解是InterpolatingFunction的形式；从中提取数据的时候，
Table[Evaluate[Re[q1[t] /. s4[[1]]]], {t, 0, 10}]
使用s4[[1]]]] 表示 s4的第一个元素，成功给出了数据。
但是很奇怪的是，当提取第2,3,4个元素的数据的时候，出现错误，没有成功。真不知
道哪里有问题， 请大家指点一下：） |
|
l*********i
发帖数: 483 | 6 我用Sol=NDSolve数值求解了一组微分方程，解是两个点的位置关于时间的函数x1[t],
y1[t],z1[t],
x2[t],y2[t],z2[t].我用下面的ParametricPlot3D命令试着把两个点的轨迹plot在一个
图里，总是
提示有错：
\!\(ParametricPlot3D[\ {Evaluate[{x\_1[t], y\_1[t], z\_1[t]} /. Sol],
Evaluate[{x\_2[t], y\_2[t], z\_2[t]} /. Sol]}, {t, 0, T\_max}]\)
哪位大大帮我看看是哪里的问题。 |
|
Q******g
发帖数: 607 | 7 tried Mathematica NDSolve?
most time, it's good enough. |
|
l*********i
发帖数: 483 | 8 因为是要嵌入到一个monte carlo程序里去的，如果用NDsolve的话不够快。 |
|
x********g
发帖数: 595 | 9 如图：
用NDSolve解完之后，解是InterpolatingFunction的形式；从中提取数据的时候，
Table[Evaluate[Re[q1[t] /. s4[[1]]]], {t, 0, 10}]
使用s4[[1]]]] 表示 s4的第一个元素，成功给出了数据。
但是很奇怪的是，当提取第2,3,4个元素的数据的时候，出现错误，没有成功。真不知
道哪里有问题， 请大家指点一下：） |
|
