function A = zRotate(phi, points)
    %n = size(points, 2)
    A = rot(phi) * points;
