z1 = 1.5 + 3i;
z2 = 1.5 + 1.5i;
disp(z1 + z2);    % 3.0 + 4.5i
disp(z1 - z2);    % 0.0 + 1.5i
disp(z1 * z2);    % -2.25 + 6.75i
disp(z1 / z2);    % 1.5 + 0.5i
disp(-z1);        % -1.5 - 3i
disp(z1');        % 1.5 - 3i
disp(abs(z1));    % 3.3541 = sqrt(z1*z1')
disp(z1 ^ z2);    % -1.10248 - 0.38306i
disp( exp(z1) );  % -4.43684 + 0.63246i
disp( imag(z1) ); % 3
disp( real(z2) ); % 1.5
%...