two-dimensional implicit plotting
implicitplot(expr, x=a..b, y=c(x)..d(x), options)
implicitplot(ineq, x=a..b, y=c(x)..d(x), options)
implicitplot(expr, x, y, options)
implicitplot(f, a..b, c..d, options)
implicitplot([expr1,expr2,t], x=a..b, y=c(x)..d(x), options)
expression or equation depending on x and y
inequality depending on x and y
equation containing procedures or operators representing a function of 2 variables
equations or expressions in x, y and t, polynomial in t
expressions that evaluate to real constants for a fixed x value
(optional) as described in the Options section, or any plot options; see plot/options
This section describes options to implicitplot.
Any additional arguments are interpreted as options which are specified as equations of the form option = value. Most options available for the plot command can also be used with the implicitplot command. For more information, see plot/options. The options that result in specific behavior when used with implicitplot are described below.
The option grid = [m, n] where m and n are positive integers (larger than 1) specifies that the points used initially to plot the 2-D curve lie on an m by n grid of equally spaced points in the ranges a..b and c..d respectively. By default a 26 by 26 grid is used, but refinements are applied and more than 262=676 points will be used typically. In contrast, the numpoints = k option controls the total number of points, so the grid will have isqrt(k+3)+1 points in each direction. If an exact value for grid is desired the default refinements should be turned off by explicitly setting the refinement options to 0: gridrefine = 0 and crossingrefine = 0.
The option gridrefine = l where l is a non-negative integer (default 2) specifies the number of recursive subdivisions to perform on any cells in the initial grid that have been determined to contain part of the curve to be plotted. This allows for a more efficient refinement of the curve, as it only computes function values where necessary. One problem is that for some curves, and some choices of the initial grid, it may not appear as though a cell contains any part of the curve, but implicitplot has some degree of intelligence in this respect, as it uses a quadratic interpolation to detect if some part of the curve may be in a cell before rejecting it. The subdivisions are performed recursively, so if the initial grid is [26,26], setting gridrefine=1 will have the highest grid level as [51,51], setting gridrefine=2 will have the highest grid level as [101,101], etc. The highest grid level for gridrefine=l will be [2^l*m,2^l*n].
The option crossingrefine = i where i is a non-negative integer (default 1) provides a facility for improving the interpolation used to determine the location where the curve crosses a cell boundary. Normally a simple linear interpolation is used, but setting crossingrefine to i>0 performs i-1 steps of bisection, and uses a quadratic interpolant to determine the cell boundary crossing. This is most useful for high degree curves where a linear interpolation does not provide a good estimate of the crossing point.
The option rational = true/false (default false), when set to true, tells implicitplot to use a rational interpolant to determine the location where the curve crosses a cell boundary. With this option implicitplot is better able to interpolate across a singular curve (for example, this provides a more accurate interpolation for the singular curve of 1/(x^2+y^2-1)). The rational interpolation is more expensive than standard interpolation as it always requires an additional function evaluation for each crossing point.
The option signchange = true/false (default true), when set to true, implicitly plots sign changes of the function. This includes, for example, the line y=0 for the input x/y. If this is not desired, and only the zeros of the input function are of interest, then signchange should be set to false. Note that when signchange is set to false, rational is automatically set to true, as it is not generally possible to detect sign changes that are not zeros unless a rational interpolant is used.
If the first argument is an expression or a procedure f, the option filledregions = true specifies that the regions defined by f = 0 should be colored. The option coloring = [c1, c2] allows the default colors to be changed, and causes region f < 0 to be given color c1 and region f > 0 to be given color c2. Note that the color of the curve defined by f = 0 defaults to black, but this can be changed with the color option.
If the first argument is an equation, such as f=g, the option filledregions = true specifies that the regions defined by f-g should be colored.
If the first argument is an inequality, and filledregions is specified, then only the region that satisfies the inequality is shaded. Just as for filledregions=false, if the inequality is strict, then the lines for f = 0 are plotted as dotted lines instead of solid lines.
The option resolution = r applies resolution-based smoothing of the output curve to reduce the size of the output data structure, but retain the visual features of the curve up to a resolution of r x r pixels. By default, this value is 0, which means that no smoothing is to be applied. Note that this option cannot be used in combination with filledregions. Specifying resolution = 1000, when used in combination with crossingrefine, can reduce the size of a high detail plot by as much as a factor of 10.
The option resultant = fast/stable applies only to the fifth calling sequence, and specifies the algorithm used to numerically compute the resultant of expr1 and expr2 with respect to t. The default is fast, which uses Euclid's algorithm to compute the numerical resultant for each evaluation. The option stable uses a full Gauss elimination of the Sylvester matrix to compute the resultant for each evaluation. Use of resultant = stable is only really helpful when the input, as polynomials in t, have numerical instability when the resultant is computed using Euclid's algorithm.
The option factor = true/false when set to true tells implicitplot to perform symbolic preprocessing to make the plotting easier, and may be the only way, for simple functions, of plotting the zeros of a function that is greater or equal to zero over the plotting region. Note that this option cannot be used for the resultant form (fifth calling sequence), operator input (no way to factor), or for inequalities. Also note that this option may disable the use of adaptive ranges (only when multiple factors are present) as the zeros are computed independently for each factor.
The default value of the factor option is false if the resultant form (fifth calling sequence) is selected, for operator input, and for inequalities. For any other input, if the factor option is not given, then factorization is attempted. However, in some cases, there is a risk that the factorization may lose some solutions due to a real product being factored into nonreal factors. Usually, Maple can determine that this is not the case. If that fails, factoring is aborted and the original expression used. If the option factor = true is explicitly given, then this test is not done and the factorization is used.
The option rangeasview = true/false (default false) is simply a convenience option that, when set to true, tells implicitplot to use the specified region (a..b,c..d) as the view option (see view in plot/options). Note that when the vertical range is dependent on x, in the form c(x)..d(x), then the vertical range is ignored for this option.
The option outlines = true/false (default false), when set to true, displays the mesh over which the sampling takes place. Note that it is only reasonable to display the mesh when using a fairly coarse grid or a low refinement.
The option adaptranges = true/false when set to true tells implicitplot that if a plot is produced where only a portion of the specified range contains curves, then the plot should be automatically recomputed with smaller ranges, thus providing greater resolution in the output plot. Note that the ranges are determined from the curves produced in the initial plot, so if the initial resolution (as specified by grid and gridrefine) is poor, then parts of the plot may be missed. This option is only effective when plotting a single function, and is ignored when multiple functions are being plotted. In the same manner it is ignored when factor = true and multiple factors are present. This option is false by default when filledregions=true, and otherwise it is true by default.
The implicitplot command computes the two-dimensional plot of an implicitly defined curve. By default, the curve is computed in Cartesian coordinates.
In the first calling sequence, implicitplot(expr, x=a..b, y=c(x)..d(x)), the equation expr must have components that are expressions in the names x and y. The expr parameter can also be an expression instead of an equation, in which case the equation expr = 0 is plotted. The range a..b must evaluate to real constants, and the range c(x)..d(x) must evaluate to real constants for any fixed value of x.
In the second calling sequence, implicitplot(ineq, x=a..b, y=c(x)..d(x)), the inequality can depend only upon x and y. If the inequality is strict, the border of the inequality is plotted as a dotted line, while for nonstrict it is plotted as a solid line. There is also a difference in the behavior of the filledregions option (described in the Options section).
The ranges x=a..b, y=c(x)..d(x) in the first and second calling sequence are optional. If they are omitted, Maple attempts to infer reasonable values for the domain of the plot. You can include a range for both x and y in the calling sequence, or neither, but you cannot select a range for only one variable.
If you omit both ranges, you can still include (both of) the variable names x and y to indicate which should be used for the horizontal and which for the vertical axis. This is demonstrated in the third calling sequence. This can be used with inequalities, too. It cannot be used if x or y do not occur in expr and x or y is also the name of a keyword plot option, so that Maple can distinguish such options.
If you don't specify the variable names x and y, they will be taken from expr. The variable used for the horizontal axis will be the one that comes first in Maple's default set ordering; for global variable names of the same length, this means it is the one that comes first in lexicographic ordering.
In the fourth calling sequence, implicitplot(f, a..b, c..d), the assumption is made that the equation f consists only of procedures or operators taking no more than two arguments. The f parameter can also be a procedure or operator instead of an equation, in which case the equation f = 0 is plotted. Operator notation must be used, that is, the procedure name is given without parameters specified, and the ranges must be given simply in the form a..b and c..d, rather than as equations.
In the fifth calling sequence, implicitplot([expr1,expr2,t], x=a..b(x), y=c..d(x)), the expr1 and expr2 must be equations or expressions in the names x,y and t, and must be polynomial with respect to t. This form of call to implicitplot is equivalent to the call implicitplot(resultant(expr1,expr2,t), x=a..b(x), y=c..d(x)), except that the resultant is computed more efficiently avoiding expression growth and excessive round-off error.
Because the implicitplot command samples the function being plotted and builds the final image from the sample, it does not, by default, detect discontinuities in the function. Instead, the function is interpolated across the discontinuities. To change this behavior, see the rational and signchange options (described in the Options section).
For similar reasons, a sample-based method cannot be used to detect isolated points which should be plotted.
If expr or f is a set or list (for the first four calling sequences), then its members are plotted together. If it is a list, then particular option values can also be given as lists, with elements corresponding to elements of expr or f. The options that can take lists as values are: color, coords, grid, legend, linestyle, numpoints, style, symbol, symbolsize, thickness, and transparency.
Note that implicitplot with the default settings does not do well for problems where the solution curve is very close to a region in which the expression being plotted becomes undefined or complex valued (for example expressions containing ln, where the argument can become negative). For problems of this class, the ideal approach is to increase the initial grid until points appear in the region, then increase gridrefine until the points describe a smooth curve.
A similar class of problematic cases is those where the only real value of the expression is 0. The only way such expressions can be dealt with is if they factor into a real and a nonreal part. An example is included below.
Plot a circle.
Use the operator form of the command.
p:= proc(x,y) if x<y then x^2+y^2-2 else 2*x+y-3 end if end proc:
Use a different coordinate system.
Use a set or a list to combine plots.
Color the regions with blue and green.
Implicit plot of 4 concentric circles benefits a lot from refinements
fn ≔ x2+y2−1⁢x2+y2−0.73⁢x2+y2−0.5⁢x2+y2−0.31:
Default plot is very smooth
With refinements disabled, it is very rough
Use of gridrefine=2 to get a finer grid (up to [101,101])
Use of crossingrefine=2 to more accurately determine crossings
Comparison of size for high detail plots using resolution
p1 ≔ implicitplot⁡fn,resolution=0:
p2 ≔ implicitplot⁡fn,resolution=1000:
Plot of a singular line using default linear interpolation
Improvement using rational interpolation
Intersection of a unit sphere and a cylinder projected into the x-y plane (two intersecting ellipses).
Plotting of a curve that is strictly ≥0
Plotting of the region where x2+y2<1 (note the dotted border)
Plotting a function where the vertical range is restricted by the horizontal variable
For the following example, we want Maple to select the domain for us, but we want to specify a non-default ordering of the variables. This can be done using the third calling sequence.
A polar plot with outlines to indicate sampling grid
In the following example, factoring is necessary because the original expression is nonreal everywhere except where it is zero. The protections outlined in the note above make it so that Maple does not factor by default, so factoring has to be turned on explicitly.
expr ≔ x3−4⁢I⁢x−2⁢x2−8⁢I⁢y−4⁢I⁢y2−2⁢y⁢x2+I⁢x⁢y2+x⁢y2+2⁢I⁢y3−2⁢y2−4⁢I⁢x2+2⁢I⁢y⁢x2−2⁢y3+I⁢x3−4⁢x+8+16⁢I+8⁢y
The command to create the plot from the Plotting Guide is
Download Help Document