The implicit equation did not work out very well, due to "poles," but here is the classic Weierstrass parameterization, as per Wikipedia.
The surface which results is the same as for the _Scherk2 long .js, .htm script, but lacks symmetry at two ends,
at least partly due to the lack of numpoints parameter in _FXGraph3d (?). So I prefer the long script, after processing, for the Scherk tower approximation.
The FXGraph3d version comes in two halves, due to the exclusion of zero, and replacement with .0001.
The .8 can be increased in value, if U and V count are increased.
The two halves can be networked and joined, offset, edge swept, and joined to a solid, for a super "pringle". A few Fillet tries did not work.

Scherk2WFirstHalf01
_FxGraph3D log((1+u*u+2*u*cos(v)) / (1+u*u-2*u*cos(v))); log((1+u*u-2*u*sin(v)) / (1+u*u+2*u*sin(v))); 2*atan((2*u*u*sin(2*v)) / (u*u*u*u - 1));0;0.8;0.0001;3.141592653589793;21;21;1;;;UV

Scherk2WSecondHalf01
_FxGraph3D log((1+u*u+2*u*cos(v)) / (1+u*u-2*u*cos(v))); log((1+u*u-2*u*sin(v)) / (1+u*u+2*u*sin(v))); 2*atan((2*u*u*sin(2*v)) / (u*u*u*u - 1));0;0.8;0.0001;-3.141592653589793;21;21;1;;;UV

- Brian