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.
[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,
fieldlines3D(X,Y,Z,U,V,W,x0,y0,z0) plots a field line with the initial point P(x0,y0,z0).The matrices X,Y,Z,U,V,and W must all be the same size and contain corresponding position and field components.The function employs the 4th order Runge-Kutta...