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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
% gauss.m: computes the same thing as gauss_demo.c
format short g
format compact
A = zeros (4,4) ;
for i = 0:3
for j = 0:3
A(i+1,j+1) = (i+1) + 1i * (2-j) ;
end
end
A (1,1) = 0 ;
A
sum(A, 'all')
A = A^2
sum (A, 'all')
D = eye (4) ;
C = D .* (A*A.')
fprintf ('C=D*A\n') ;
D = diag (diag (A)) ;
C = D*A
fprintf ('C=A*D\n') ;
C = A*D
C = ones (4,4) * (1 - 2i) ;
C = C + A*A.'
B = ones (4,4) * (1 - 2i) ;
C = C + B*A
C = C + A*B
C = (1-2i) + A
C = A * (1-2i)
C = A.' * (1-2i)
C = (1-2i) * A.'
R = real (C)
R = real (C')
S = zeros (4,4) ;
for i = 0:3
for j = 0:3
S(i+1,j+1) = i - j ;
end
end
C
R = real (C) + 1 + S
|
