| Author |
Comment/Response |
Martin
|
 12/11/06 5:55pm
I have some complicated deformed dipolar magnetic field in cartesian components, and I need to show some field lines. I'll also have to extract coordinates of lines to a list. How should I do that ?
I can use the code below (with properly defined functions fx, fy and fz), but it's unreliable (it frequently hit the singularity at the center) and it frequently gives many unpredictable curves. Any idea ? Please, I need help ! :-(
FieldCurve = NDSolve[{
x'[t] == fx[ x[t], y[t], z[t] ],
y'[t] == fy[ x[t], y[t], z[t] ],
z'[t] == fz[ x[t], y[t], z[t] ],
x[0] == ... number ...,
y[0] == ... number ...,
z[0] == ... number ...
}, {x, y, z}, {t, 0, 100}, MaxSteps -> 10000]
Graph = ParametricPlot3D[
Evaluate[{x[t], y[t], z[t]}/. FieldCurve],
{t, 0, 100}, PlotPoints -> 1000]
URL: , |
|
