1function rotZ = rotationZ(angle) 2 3cos_angle = cos(angle); 4sin_angle = sin(angle); 5 6rotZ = [ cos_angle -sin_angle 0; sin_angle cos_angle 0; 0 0 1 ]; 7 8end