%% d apres un code original de Juergen Gilg
%% http://melusine.eu.org/lab/bpst/pst-solides3d/surfaces/ex_6.pst

-4 4 setxrange
-4 4 setyrange
60 setxunit

15 -20 0 SetCamPos
0 0 0 SetCamView

-10 10 setxrange3d
-3 3 setyrange3d

GetCamPos setlightsrc
1.5 setlightintensity

2 setlinejoin 
20 setfontsize
setTimes
.01 setlinewidth

/aretescachees false def

/klein {
3 dict begin
   /v exch def
   /u exch def
#rpn# (3*(1+Sin(v))+2*(1-Cos(v)/2)*Cos(u))*Cos(v)
#rpn# -2*(1-Cos(v)/2)*Sin(u)
#rpn# (4+2*(1-Cos(v)/2)*Cos(u))*Sin(v)
end
} def

.1 setlinewidth
gris

0 pi 2 mul 0 pi 2 mul [50 dup] {klein} newsurfaceparametree
dup videsolid
dup [.5 .6] solidputinouthuecolors
%dup (.8 setgray) inputcolors
drawsolid**