f(x,y) = y - x**2 / tan(y)
set contour base
set cntrparam levels discrete 0.0
set nosurface
set term table
set out 'curve.dat'
splot f(x,y)
set out
set term postscript
set output "curve.ps"      
plot 'curve.dat' w l