Blame view

Giac_maj/giac-1.4.9/check/TP03.tst 2.74 KB
6663b6c9   adorian   projet complet av...
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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
  [],
  "Warning: some commands like subs might change arguments order",
  0,0,0,1,0,1e-10,10,[1,50,0,25],0,0,0,
  3,
  x,
  2,
  2,
  6,
  1,
  1,
  [[1,1]],
  [[3,0,0,0],[0,6,0,0],[0,0,18,0],[0,0,0,36]],
  "Done",
  0,
  "Done",
  matrix[[0,0,0,0],[0,0,0,0],[0,0,0,0],[0,0,0,0]],
  proc(A) 
    local m,u,v,i,j; 
    m:=0;  
    u:=0;  
    v:=0;  
    for i from 1 to (dim(A))[1] do for j from 1 to (dim(A))[2] do if ((abs(A[i,j]))>0) and ((m=0) or ((abs(A[i,j]))<m)) then  
      m:=abs(A[i,j]); 
      u:=i; 
      v:=j 
    fi ; 
     od;; 
     od;;  
    u,v;  
   
  end;,
  matrix[[6,0,7,9],[6,8,6,9],[3,2,4,6],[3,2,3,3]],
  3,2,
  3,2,
  matrix[[1,0,0,0],[0,1,0,0],[0,0,1,0],[0,0,0,1]],
  "Done",
  "Done",
  [[1,0,0,0],[-1,1,-2,-3],[0,0,1,0],[0,0,0,1]],
  matrix[[6,0,7,9],[-2,8,-10,-15],[1,2,0,0],[1,2,-1,-3]],
  proc(A,i,j) 
    local n,U,V,l; 
    n:=(dim(A))[1];  
    U:=identity(n);  
    V:=identity(n);  
    for l from 1 to n do V[l,i]:=-iquo(A[l,j],A[i,j]); 
     od;;  
    V[i,i]:=1;  
    for l from 1 to n do U[j,l]:=-iquo(A[i,l],A[i,j]); 
     od;;  
    U[j,j]:=1;  
    V*A*U;  
   
  end;,
  matrix[[6,1,7,9],[6,8,6,9],[3,1,4,6],[3,2,1,3]],
  matrix[[0,-1,-1,-3],[0,6,-2,-3],[3,1,1,0],[0,1,-3,-3]],
  proc(A) 
    k:=0;  
    n:=(dim(A))[1];  
    l:=[seq(0,k=(1 .. n))];  
    B:=A;  
    i,j:=minval(B);  
    while (([i,j])<>[0,0]) and ((nops(B))>1) do  
      B:=trans(B,i,j); 
      if [i,j]=[minval(B)] then  
        k:=k+1; 
        l[k]:=B[i,j]; 
        B:=delcols(B,j .. j); 
        B:=delrows(B,i .. i) 
      fi ; 
      i,j:=minval(B); 
     od;;  
    diag([seq(l[k],k=(1 .. (n-1))),B[1,1]]);  
   
  end;,
  matrix[[2,2,2],[6,12,6],[6,4,6]],
  [[2,0,0],[0,-2,0],[0,0,0]],
  "Done",
  [[2,2,2],[6,12,6],[6,4,12]],
  [[2,0,0],[0,-2,0],[0,0,18]],
  [[2,0,0],[0,2,0],[0,0,18]],[[2,0,0],[0,2,0],[0,0,18]],
  matrix[[4,0,0],[0,6,0],[0,0,8]],
  [[2,0,0],[0,4,0],[0,0,24]],
  [[4,0,0],[0,6,0],[0,0,8]],
   (i,j)->if (i-j)*(i-1)=0 then 1; 
  fi ,
  matrix[[1,1,1],[0,1,0],[0,0,1]],
  matrix[[4,6,8],[0,6,0],[0,0,8]],
  [[2,0,0],[0,8,0],[0,0,-12]],
  proc(A,i,j) 
    local n,U,l; 
    n:=(dim(A))[1];  
    U:=identity(n);  
    for l from 1 to n do U[j,l]:=-iquo(A[i,l],A[i,j]); 
     od;;  
    U[j,j]:=1;  
    A*U;  
   
  end;,
  proc(A) 
    k:=0;  
    n:=(dim(A))[1];  
    l:=[seq(0,k=(1 .. n))];  
    B:=A;  
    i,j:=minval(B);  
    while (([i,j])<>[0,0]) and ((nops(B))>1) do  
      B:=transC(B,i,j); 
      if [i,j]=[minval(B)] then  
        k:=k+1; 
        l[k]:=B[i,j]; 
        B:=delcols(B,j .. j); 
        B:=delrows(B,i .. i) 
      fi ; 
      i,j:=minval(B); 
     od;;  
    diag([seq(l[k],k=(1 .. (n-1))),B[1,1]]);  
   
  end;,
  proc(A) 
    n:=(dim(A))[1];  
    d:=Zequiv(A);  
    L:=[];  
    for i from 1 to n-1 do  
    T:=matrix(n+1-i,n+1-i,f); 
    d:=ZequivC(T*d); 
    L:=[op(L),d[1,1]]; 
    d:=delrows(delcols(d,1 .. 1),1 .. 1); 
     od;;  
    [op(L),d[1,1]];  
   
  end;,
  matrix[[2,2,2],[6,12,6],[6,4,12]],
  [2,-2,18],
  [[2,0,0],[0,2,0],[0,0,18]],
  [2,4,48],
  [[2,0,0],[0,4,0],[0,0,48]]