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
|
function t=trace(a)
// trace - computes the trace of a matrix
// Copyright INRIA
select type(a)
case 1 then
[m,n]=size(a)
if m<>n then error('the argument of trace must be a square matrix');end
t=sum(diag(a))
case 2 then
[m,n]=size(a)
if m<>n then error('the argument of trace must be a square matrix');end
t=sum(diag(a))
//-compat next case retained for list /tlist compatibility
case 15 then
if a(1)=='r' then
[m,n]=size(a)
if m<>n then error('the argument of trace must be a square matrix');end
t=sum(diag(a))
else
error('invalid argument for trace');
end
case 16 then
if a(1)=='r' then
[m,n]=size(a)
if m<>n then error('the argument of trace must be a square matrix');end
t=sum(diag(a))
else
error('invalid argument for trace');
end
case 5 then
[m,n]=size(a)
if m<>n then error('the argument of trace must be a square matrix');end
t=sum(diag(a))
else
error('invalid argument for trace');
end
|
