![planeplotter graph planeplotter graph](https://keptenkurk.files.wordpress.com/2014/10/ppsnip4.jpg)
NOTE: If this page is misbehaving, it may be that your browser is still using an old version of one of the support files. Thanks to many users for helpful feedback, especially to Larry Friesen at Butler Community College, who suggested many improvements to this page, and (with his colleagues and students) tested it extensively. In addition, your Web browser's "zoom" features will work with this page (often not the case with Java), so you can increase the display size if you need to do a demonstration for a large class. (In particular, I never know if a particular classroom's computer will be able to run it.) This page is written in JavaScript (not Java), so it does not have those compatibility and security issues. I had used his Java applet for many years in my Calculus 2 and DE classes, but began to have issues with unstable Java installations and nuisance security alerts. I created this page as a replacement for the very nice JODEapplet by Marek Rychlik. The value of k defaults to 0.01 if omitted.Ībout this page: Acknowledgments, release history, basic troubleshooting delta(x,k) approximates the Dirac Delta function using the Gaussian function exp(-(x/k) 2)/(|k|√π).step(x,a) is the unit step function, equivalent to when(x>a,1,0).This can be used to create piecewise-defined functions, such as when(x>0,x^2,y). if(condition, true-value, false-value) or when(condition, true-value, false-value).In addition to many standard functions (and some exotic ones), the following functions can be helpful for some DE modeling problems:.For example, I once solved an IVP for which the constant was C=2e+1, and realized that I got the correct plot when I entered … +(1+2e)e^(-x), while … +(2e+1)e^(-x) produced the wrong graph, because 2e+1 is interpreted 2×10 1=20. e can be entered as a simple constant (as can pi), but take care, because e also is used in exponential notation.For numerical input (such as the coordinates for an initial value), fractions and complex expressions are allowed for example, you can enter 5/3 instead of 1.6666666667, or sqrt(2), or pi/2.Absolute values can be entered as either |x| or abs(x).However, some alternate notations are also accepted (for example, arctan). The inverse trig functions should be entered as atan, asin, and acos (and similarly for inverse hyperbolic trig).Exponentiation (like 7 x) can be entered as either 7^x or power(7,x), andĮ x can be entered as e^x or exp(x).All functions must have parentheses-for example, use sin(x) rather than sin x, and ln(|x|) rather than ln|x|.Closing parentheses are not optional (unlike, say, on TI-84 graphing calculators).x cos(x) or x y (There must be a space between x and cos, or x and y.).Multiplication is implied in expressions like:.
Planeplotter graph how to#
For example, for the DE dy/dx = - y/x (a circle), here are solution curves for RK4 with h=0.05 without switching (left) and with switching (right): How to enter expressionsįor the most part, expressions are entered using standard mathematical notation, with a few caveats: Specifically: If, at any point, | dy/dx| > 3 (i.e., if the tangent lines get too steep), the method switches the roles of x and y. It affects all of the numerical methods for ODEs except for RKF (it has no effect on solutions for systems). The "switching" option next to the choice of method is an adaptation that produces better solution plots in some cases. Runge-Kutta-Fehlberg (RKF) BETA - only working for ODEs.Runge-Kutta (both RK4 and the "3/8 rule").(In the case of non-autonomous systems-that is, where either dx/dt or dy/dt depends on t-the direction field shown is for t = 0.)īy specifying initial values, users can see approximate solution curves, with several choices for the solution method (click links to read more at Wikipedia): (Note: A limited number of alternative variables can be chosen, to make it easier to adapt to different applications or textbook conventions.)įor ODEs, a slope field is displayed for systems, a direction field is shown. Users enter a first-order ODE in the form dy/dx = f( x, y), or a system in the form dx/dt = f( t, x, y) and dy/dt = g( t, x, y).