limit of a function is giving another limit
3 views (last 30 days)
Show older comments
I tried calculating the limit of the function below and matlab gave me the answer as another limit. Would i need to change the function so that i get a better answer? If I do need to change it please let me know how I could do it.
fun2 = 0.5.*(1./(y.^2)).*besselj(0,(2.*pi./lambda).*(1./y^0.5)).*log((-z_0.*y^0.5+sqrt((z_0.^2).*(y)+1))./(-z_1.*y^0.5+sqrt((z_1.^2).*y+1)));
limit(fun2,y,0)
ans =
limit((log(((y + 1)^(1/2) - y^(1/2))/((4*y + 1)^(1/2) - 2*y^(1/2)))*besselj(0, pi/(50*y^(1/2))))/y^2, y, 0)/2
0 Comments
Answers (2)
Sulaymon Eshkabilov
on 9 Jun 2021
Because Symbolic Math Toolbox could not solve your exercise and find an explcit solution formulation.
Steven Lord
on 10 Jun 2021
Let's look at your function.
syms y
z_0 = 2;
z_1 = 3;
lambda = 4;
fun2 = 0.5.*(1./(y.^2)).*besselj(0,(2.*pi./lambda).*(1./y^0.5)).*log((-z_0.*y^0.5+sqrt((z_0.^2).*(y)+1))./(-z_1.*y^0.5+sqrt((z_1.^2).*y+1)));
fplot(fun2, [-1, 1])
That doesn't look all that promising in the vicinity of x = 0 for the limit to exist. Let's zoom in a bit.
figure
fplot(fun2, [-0.001 0.001])
Yeah, that really looks to me like the limit you're trying to compute does not exist which means you're encountering the first of Cleve's Golden Rules of computation.
2 Comments
Steven Lord
on 10 Jun 2021
syms y
z_0 = 2;
z_1 = 3;
lambda = 4;
fun2 = 0.5.*(1./(y.^2)).*besselj(0,(2.*pi./lambda).*(1./y^0.5)).*log((-z_0.*y^0.5+sqrt((z_0.^2).*(y)+1))./(-z_1.*y^0.5+sqrt((z_1.^2).*y+1)))
As y goes to 0, goes to Inf. How does behave as z goes to Inf? It oscillates according to the DLMF. That oscillation is multiplied by which is also going to Inf. So the appearance of the plot makes sense to me looking at your equation. If you were to evaluate fun2 for smaller and smaller values of y:
yvec = 10.^(-1:-1:-10).';
vpa(subs(fun2, y, yvec))
This is in agreement with the plot as well.
See Also
Categories
Find more on Calculus in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!