Upssss...
at the beginning of the prog instead of z:=sqrt(2)/2; there is z:=1-sqrt((3-2*sqrt(2))/2);