Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
1995
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 1995

[Date Index] [Thread Index] [Author Index]

Search the Archive

DirectionFieldPlot[]

  • To: mathgroup at christensen.cybernetics.net
  • Subject: [mg411] DirectionFieldPlot[]
  • From: Xah Y Lee <xyl10060 at fhda.edu>
  • Date: Mon, 16 Jan 1995 21:17:30 -0800 (PST)

MathGroupers,

I am taking my first Differential Equations course and encountered Direction 
Field Plot.

I wrote a mma program to do the plot. Included below.

Execute the following and you'll see the plot.

---------------------------------
(*variables used*)
Clear[xMin, xMax, xStep, yMin, yMax, yStep, ll, fun]
Clear[gridList, slopeList, destList, GP1, g1]

{xMin,xMax,xStep} = {-3,3,.2};   (*specify the grid points*)
{yMin,yMax,yStep} = {-3,3,.2};
ll = 1 (*ll = line length*);      (* specify the slop length to be drawn*)
fun = (#1^2-#2)&;                 (* specify the function to be ploted*)


gridList = N@Flatten[ Table[{x,y}, 
  {x,xMin,xMax,xStep},{y,yMin,yMax,yStep}], 1];
slopeList = fun[Sequence@@#]& /@ gridList;
destList = gridList + ll * ( {Cos@#,Sin@#}& /@ (ArcTan /@ slopeList) );


GP1 = N @ Line /@ Transpose[{gridList, destList}];
g1 = Graphics[{Hue[0.00], GP1 }, Axes->True, AspectRatio->Automatic];
Show[ g1 ]

(*GP stand for GraphicPrimitive. g1 is an graphic object that can be used 
with Show[]*)

---------------------------------
I am going to modify this into a package that contain a function 
DirectionFieldPlot[] that can be used as easily as Plot, with online 
documentations. I may upload this to MathSource sometimes in the future.

I wrote another program that do Trochoid animations. MathSource item number:
0207-234.

 Xah Lee      xyl10060 at tiptoe.fhda.edu    74631.731 at compuserve.com
 A stud at Foothill College.
 Mountain View, CA
 Quote of the day: Can you tell me how to get to Sesame Street?



  • Prev by Date: Shortening Polynomials
  • Next by Date: Fourier transforming TIFF (2D) images
  • Previous by thread: Re: Shortening Polynomials
  • Next by thread: Fourier transforming TIFF (2D) images