Reference: computeFlowTrajectory
computeFlowTrajectory(expr, var, r, t0, t1, δt) returns the curve in n-dimensional space traced out by a particle that has the velocity expr(var) at point var in space. expr is thus a real n-dimensional vector-valued function in var, a valid identifier for a n-dimensional position vector. r ∈ R^n is the initial position of the particle. t0 is the initial time of the simulation, and t1 is the final time. δt is the temporal resolution.
Example: a ≔ 2
b ≔ 3.001
vectorField ≔ createVectorField("❨1 + a⋅x^2⋅y − b⋅x−x, -a⋅x^2⋅y + b⋅x❩", "x, y", [0, 10, 0.25]^2)
drawVectorField("vectorField", "colour:#333333")
flow ≔ computeFlowTrajectory("❨1 + a⋅r_1^2⋅r_2 − b⋅r_1 − r_1, -a⋅r_1^2⋅r_2 + b⋅r_1❩", "r", ❨1, 4❩, 0, 100, 0.01)
drawLines("flow", "colour:gold")