Question: Discontinuities... I think

Hello, again,
basically, I have an implicit relation for lambda:

> eq1 := (sqrt(lambda*alpha)+alpha)*BesselJ(alpha, 2*sqrt(lambda*alpha))/(2*sqrt(lambda*alpha)) = BesselJ(1+alpha, 2*sqrt(lambda*alpha))/alpha;

So solutions of lambda, for a particular alpha are given at the intersections of the two functions:

> plot({subs(alpha = 2, (sqrt(lambda*alpha)+alpha)*BesselJ(alpha, 2*sqrt(lambda*alpha))/(2*sqrt(lambda*alpha))), subs(alpha = 2, BesselJ(1+alpha, 2*sqrt(lambda*alpha))/alpha)}, lambda = 0 .. 50, -1 .. 1, color = [red, blue])

I'm interested in the first intersection ( the smallest non zero lambda)... so for alpha=2:

> fsolve(subs(alpha = 2, eval(eq2)), 0 .. 10);

But I really want is to plot these smallest lambdas as a function fo alpha, so 

> eq2 := unapply(eq1, lambda);

> plot([t-> fsolve(subs(alpha = t, eval(eq2))) ], 0 .. 10, 10 .. 10);

... but then there are some sort of discontinuities between 2-3. I don't think that they should be there, have I missed out something that is leaving maple confused? Is there anyway of removing them?

Thanks k

Please Wait...