1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
clear;
clc;
function [r,th]=rect2pol(x,y)
//rectangle to polar coordinate conversion
r=sqrt(x^2+y^2);
th = atan(y,x)*180/%pi;
endfunction
// linear cominations
// 1 part (1+a)
a=complex((-0.5),0.866);
b=1+a;
[r,th]=rect2pol(real(b),imag(b));
mprintf("a) magnitude = %f, Angle = %f\n",r,th);
// 2 part (a^2 - 1)
c=((a^2)-1)
[r,th]=rect2pol(real(c),imag(c));
mprintf(" b) magnitude = %f, Angle = %f\n",r,th);
// 3 part (a^2 + a)
d=((a^2)+a)
[r,th]=rect2pol(real(d),imag(d));
mprintf(" c) magnitude = %f, Angle = %f\n",r,th);
// 4 part (a^2 + a + 1)
e=((a^2)+a+1)
[r,th]=rect2pol(real(e),imag(e));
mprintf(" d) magnitude = %f, Angle = %f",r,th);
|
