Bivariate Limits
The limit command has been enhanced for the case of limits of bivariate rational functions with non-isolated singularities. Many such limits that could not be determined previously are now computable. If the limit exists in such a situation, it is either +∞ or −∞. Maple can also determine if the limit does not exist, and then returns undefined. In Maple 18, all the following limit calls would return unevaluated, but they can now be computed in Maple 2015.
f≔x yx+y:
limitf,x=0,y=0
undefined
g≔4 x yx−y2:
limitg,x=0,y=0
h≔x4−x2−y2x−y4:
limith,x=0,y=0
−∞
Let us plot these three functions in the neighborhood of the origin.
In the first example, f tends to +∞ on one side of the singularity y=−x and to −∞ on the other side (as shown in the following plot). Therefore, the limit at the origin does not exist.
pf1≔plot3df, x=−0.1..0.1,y=−0.1..−x−1e−10,axes=boxed,view=−10..10:
pf2≔plot3df, x=−0.1..0.1,y=−x+1e−10..0.1,axes=boxed,view=−10..10:
plots:-displaypf1,pf2
Now, consider the second example.
s≔plots:-spacecurvex,−x,−1,x=−0.1..0.1,color=red,thickness=3:pg≔plot3dg,x=−0.1..0.1,y=−0.1..0.1,axes=boxed,view=−10..100,numpoints=40000:plots:-displays,pg
g tends to +∞ close to the singularity y=x. However, along the anti-diagonal y=−x, the limit is finite:
evalg,y=−x
−1
Thus, g does not have a limit at the origin. In fact, any number≥−1 can occur, namely, as the limit along the ray y=a x for −1≤a<1:
evalg,y=a x
4⁢x2⁢a−a⁢x+x2
limit,x=0
4⁢aa−12
plot,a=−3..1.1,view=−2..10
In the last example, h tends to −∞ on both sides of the singularity y=x.
plot3dh,x=−0.1..0.1,y=−0.1..0.1,axes=boxed,view=−100000..1
However, in this case the limit along any ray y=a x with a≠1 is −∞ as well.
evalh,y=a x
−a2⁢x2+x4−x2−a⁢x+x4
−signum⁡a2+1a−14⁢∞
factor
simplify assuming a<1;simplify assuming a>1
You can prove that the limit exists and is −∞ for any curve approaching the origin by using the theory of Lagrange multipliers. The extremal values (maxima and minima) of the function h on the circle, for a fixed radius r, satisfy the condition that the gradient of the function and the gradient of the constraint equation of the circle are parallel:
C≔x2+y2−r2:
withVectorCalculus:
df≔normal~Jacobianh,x,y
df≔−2⁢2⁢x3⁢y−x2−y⁢x−2⁢y2x−y52⁢2⁢x4−2⁢x2−y⁢x−y2x−y5
dC≔JacobianC,x,y
dC≔2⁢x2⁢y
eq≔factornumernormaldf1,1dC1,2−df1,2dC1,1
eq≔−8⁢x3−x−y⁢x2+y2
Thus, the local maximum and minimum values of f on C occur when both C=0 and eq=0. For the bivariate limit, this means that it is sufficient to consider only those critical paths satisfying eq=0. However, you also need to consider the global suprema and infima, which may occur close to the singularity y=x. In the example, the factor x2+y2 of eq does not admit any real paths, so there is only one critical path given by x3−x−y=0 (or, equivalently, y=−x+x3).
normalevalh,y=−x+x3
−x2−1x2⁢x2−23
In order to certify the limit close to the singularity y=x as well, you cannot take the limit along the singularity. Instead, consider two curves that approach the singularity closely from the top and from the bottom, respectively:
c1≔y=x+x2
c1≔y=x2+x
c2≔y=x−x2
c2≔y=−x2+x
plotx,rhsc1,rhsc2,−x+x3,x=−0.5..0.5, legend=singular path,c1,c2,critical path
normalevalh,c1
−2⁢x+1x6
normalevalh,c2
2⁢x−1x6
See Also
updates/Maple17/BivariateLimits, limit, limit/multi - multidimensional limits
Download Help Document