Return to net CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

Annotation of researchv10dc/cmd/matlab/exec/net, revision 1.1.1.1

1.1       root        1: C = <
                      2: 1 2 15 . . .
                      3: 2 1 3 . . .
                      4: 3 2 4 11 . .
                      5: 4 3 5 . . .
                      6: 5 4 6 7 . .
                      7: 6 5 8 . . .
                      8: 7 5 9 30 . .
                      9: 8 6 9 10 11 .
                     10: 9 7 8 30 . .
                     11: 10 8 12 30 31 34
                     12: 11 3 8 12 13 .
                     13: 12 10 11 34 36 .
                     14: 13 11 14 . . .
                     15: 14 13 15 16 38 .
                     16: 15 1 14 . . .
                     17: 16 14 17 20 35 37
                     18: 17 16 18 . . .
                     19: 18 17 19 . . .
                     20: 19 18 20 . . .
                     21: 20 16 19 21 . .
                     22: 21 20 22 . . .
                     23: 22 21 23 . . .
                     24: 23 22 24 35 . .
                     25: 24 23 25 39 . .
                     26: 25 24 . . . .
                     27: 26 27 33 39 . .
                     28: 27 26 32 . . .
                     29: 28 29 32 . . .
                     30: 29 28 30 . . .
                     31: 30 7 9 10 29 .
                     32: 31 10 32 . . .
                     33: 32 27 28 31 34 .
                     34: 33 26 34 . . .
                     35: 34 10 12 32 33 35
                     36: 35 16 23 34 36 .
                     37: 36 12 35 38 . .
                     38: 37 16 38 . . .
                     39: 38 14 36 37 . .
                     40: 39 24 26 . . . 
                     41: >;
                     42: 
                     43: <n, m> = size(C);
                     44: A = 0*ones(n,n);
                     45: for i=1:n, for j=2:m, k=c(i,j); if k>0, a(i,k)=1;
                     46: check = norm(A-A',1), if check > 0, return
                     47: <X,D> = eig(A+eye);
                     48: D = diag(D);  D = D(n:-1:1)
                     49: X = X(:,n:-1:1);
                     50: <x(:,1)/sum(x(:,1)) x(:,2) x(:,3) x(:,19)>

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.