Mathematical Principles for Scientific Computing and Visualization
http : // www.farinhansford.com/books/scv
by Gerald Farin and Dianne Hansford

This notebook illustrates how you plot a vectorfield.
Figure 9.1 was created with this notebook.

p(t) is population at a given time
p'(t) = c * p(t) describes population growth.
c is constant the depends on the particular population.
If we know the population at some time t_0, an initial value of p_0, then what is the population at some time t? 

Below we have plotted the slope p'(t) for pairs (t,p). The circles represent initial values. These initial values determine functions that follow the vectorfield.

Initial value ODE

See Chapter 9 : Computing Dynamic Processes for details.

Updated July 2008

(* Set the output directory *)
SetDirectory[ToFileName[Extract["FileName"/.NotebookInformation[EvaluationNotebook[]],{1},FrontEnd`FileName]]];

(* Load the vectorfield package *)
Needs["VectorFieldPlots`"];

(* Plot a vectorfield *)
vecfield=VectorFieldPlot[{1,Sin[2 y]-0.4` (x^2)},{x,-2,2},{y,-2,2},PlotPoints->30,Frame->True];


startpoint1=ListPlot[{{-1,1}},Axes->False,PlotStyle->PointSize[0.02`]];

startpoint2=ListPlot[{{-1,-1.5`}},Axes->False,PlotStyle->PointSize[0.02`]];

pic=Show[{vecfield,startpoint1,startpoint2}]

Export["slopefield_problem.eps",pic];