Re: Maple: proc() and plot of implicit function

From: Carl Devore (devore_at_math.udel.edu)
Date: 11/30/04


Date: Tue, 30 Nov 2004 11:27:34 -0500

On 30 Nov 2004, Colin Rowat wrote:
> I define the function using conditionals in a proc() statement,
> resolving the implicit component with evalf(solve()).

There is a Maple-specific newsgroup comp.soft-sys.math.maple.

Your problem is that your solve's return 3 solutions, only one of which is
real. You always select the 2nd solution. There is no reliable
significance to the order that Maple returns solutions. Sometimes the 2nd
solution is not real. Instead of using evalf(solve(...)), use
fsolve(...). This will return the one real solution. In the next line,
where you select the 2nd solution with sol(k1)[2], remove the [2].

Also, fsolve is much faster than solve.