Hello,
> What happened is that f(x, y) was evaluated for x=0 and y=0,
> which caused a divison by 0, which was not handled gracefully
> by plot3d. Sorry for the bother; this is a bug in the plotting code.
> (Specifically a function named COERCE-FLOAT-FUN.)
This is also a bug to be fixed in 'gpdraw' (ok, ok, let's call it 'draw'
in the near future :), since it also calls 'coerce-float-fun'.
When the plotting program finds a pole in function f(x), say in x0, my
present idea is to re-evaluate f in some point in the neighborhood of
x0, somewhere in (x0-eps/k,x0+eps/k)-{x0}, for some k>1. Andrej, in
private e-mail (I must say he is an excellent bug-hunter), suggested to
compute some average.
In explicit functions in 2d, we could also remove this conflictive
point, but in 3d the problem is harder, since points must be sent to
gnuplot in groups of equal size separated by blank lines.
Any other ideas?
--
Mario Rodriguez Riotorto
www.biomates.net