fieldlines

Plots a field line with the initial point P(x0,y0) for a 2D-field.

Plots a field line with the initial point P(x0,y0) for a 2D-field given by an equaly spaced grid of positions and corresponding field components. The function employs the 4th order Runge-Kutta method.

Example:
[x,y] = meshgrid(-2:.2:2,-1:.15:1);
z = x .* exp(-x.^2 - y.^2);
[px,py] = gradient(z,.2,.15);
quiver(x,y,px,py), axis image,
hold on
fieldlines(x,y,px,py,-0.7,0.2)