my $z = sqrt( (1/(0.5 + 0.02*$x**2)*cos($x))**2
            + (1/(0.5 + 0.02*$y**2)*cos($y))**2 );