slope field

{{short description|Visual representation of solutions to a differential equation}}

File:Slope Field.png

A slope field (also called a direction field{{cite book |last1=Boyce |first1=William |title=Elementary differential equations and boundary value problems |date=2001 |publisher=Wiley |isbn=9780471319993 |page=3 |edition=7}}) is a graphical representation of the solutions to a first-order differential equation{{cite book|author=Vladimir A. Dobrushkin|title=Applied Differential Equations: The Primary Course|url=https://books.google.com/books?id=d-5MBgAAQBAJ&pg=PA13|year=2014|publisher=CRC Press|isbn=978-1-4987-2835-5|page=13}} of a scalar function. Solutions to a slope field are functions drawn as solid curves. A slope field shows the slope of a differential equation at certain vertical and horizontal intervals on the x-y plane, and can be used to determine the approximate tangent slope at a point on a curve, where the curve is some solution to the differential equation.

Definition

=Standard case=

The slope field can be defined for the following type of differential equations

:y' = f(x, y),

which can be interpreted geometrically as giving the slope of the tangent to the graph of the differential equation's solution (integral curve) at each point (x, y) as a function of the point coordinates.{{cite book|author1=Andrei D. Polyanin|author2=Alexander V. Manzhirov|title=Handbook of Mathematics for Engineers and Scientists|url=https://books.google.com/books?id=ge6nk9W0BCcC&pg=PA453|year=2006|publisher=CRC Press|isbn=978-1-58488-502-3|page=453}}

It can be viewed as a creative way to plot a real-valued function of two real variables f(x,y) as a planar picture. Specifically, for a given pair x,y, a vector with the components [1, f(x,y)] is drawn at the point x,y on the x,y-plane. Sometimes, the vector [1, f(x,y)] is normalized to make the plot better looking for a human eye. A set of pairs x,y making a rectangular grid is typically used for the drawing.

An isocline (a series of lines with the same slope) is often used to supplement the slope field. In an equation of the form y'=f(x,y), the isocline is a line in the x,y-plane obtained by setting f(x,y) equal to a constant.

=General case of a system of differential equations=

Given a system of differential equations,

:\begin{align}

\frac{dx_1}{dt}&=f_1(t, x_1, x_2, \ldots, x_n) \\

\frac{dx_2}{dt}&=f_2(t, x_1, x_2, \ldots, x_n) \\

&\;\;\vdots \\

\frac{dx_n}{dt}&=f_n(t, x_1, x_2, \ldots, x_n)

\end{align}

the slope field is an array of slope marks in the phase space (in any number of dimensions depending on the number of relevant variables; for example, two in the case of a first-order linear ODE, as seen to the right). Each slope mark is centered at a point (t,x_1,x_2,\ldots,x_n) and is parallel to the vector

:\begin{pmatrix} 1 \\ f_1(t,x_1,x_2,\ldots,x_n) \\ f_2(t,x_1,x_2,\ldots,x_n) \\ \vdots \\ f_n(t,x_1,x_2,\ldots,x_n) \end{pmatrix}.

The number, position, and length of the slope marks can be arbitrary. The positions are usually chosen such that the points (t,x_1,x_2,\ldots,x_n) make a uniform grid. The standard case, described above, represents n=1. The general case of the slope field for systems of differential equations is not easy to visualize for n>2.

General application

With computers, complicated slope fields can be quickly made without tedium, and so an only recently practical application is to use them merely to get the feel for what a solution should be before an explicit general solution is sought. Of course, computers can also just solve for one, if it exists.

If there is no explicit general solution, computers can use slope fields (even if they aren’t shown) to numerically find graphical solutions. Examples of such routines are Euler's method, or better, the Runge–Kutta methods.

Software for plotting slope fields

Different software packages can plot slope fields.

=Direction field code in [[GNU Octave]]/[[MATLAB]] =

funn = @(x, y)y-x; % function f(x, y) = y-x

[x, y] = meshgrid(-5:0.5:5); % intervals for x and y

slopes = funn(x, y); % matrix of slope values

dy = slopes ./ sqrt(1 + slopes.^2); % normalize the line element...

dx = ones(length(dy)) ./ sqrt(1 + slopes.^2); % ...magnitudes for dy and dx

h = quiver(x, y, dx, dy, 0.5); % plot the direction field

set(h, "maxheadsize", 0.1); % alter head size

= Example code for [[Maxima (software)|Maxima]] =

/* field for y'=xy (click on a point to get an integral curve). Plotdf requires Xmaxima */

plotdf( x*y, [x,-2,2], [y,-2,2]);

= Example code for [[Mathematica]] =

(* field for y'=xy *)

VectorPlot[{1,x*y},{x,-2,2},{y,-2,2}]

= Example code for [[SageMath]]<ref>{{Cite web|url=https://doc.sagemath.org/html/en/reference/plotting/sage/plot/plot_field.html|title = Plotting fields — Sage 9.4 Reference Manual: 2D Graphics}}</ref> =

var('x,y')

plot_slope_field(x*y, (x,0,-5), (y,-5,5))

Examples

Image:Slope_field_1.svg|Slope field

Image:Slope_field_with_integral_curves_1.svg|Integral curves

image:Isocline_3.png|Isoclines (blue), slope field (black), and some solution curves (red)

See also

References

{{Reflist}}

  • Blanchard, Paul; Devaney, Robert L.; and Hall, Glen R. (2002). Differential Equations (2nd ed.). Brooks/Cole: Thompson Learning. {{ISBN|0-534-38514-1}}