%Vector Field Plotter function vectorField(tmax,xmin,xmax,n) dx = (xmax-xmin)/(n+1); dt = (tmax)/(n+1); [t,x] = meshgrid(0:dt:tmax,xmin:dx:xmax); vt = ones(n+2); vx = f(x,t); plt=quiver(t,x,vt,vx,1,'.') set(plt,'Color','red')