> with(linalg): Warning, new definition for norm Warning, new definition for trace > A:=matrix(3,3); A := array(1 .. 3, 1 .. 3, []) > evalm(det(A) * A^(-1)); [A[2, 2] A[3, 3] - A[2, 3] A[3, 2] , -A[1, 2] A[3, 3] + A[1, 3] A[3, 2] , A[1, 2] A[2, 3] - A[1, 3] A[2, 2]] [-A[2, 1] A[3, 3] + A[2, 3] A[3, 1] , A[1, 1] A[3, 3] - A[1, 3] A[3, 1] , -A[1, 1] A[2, 3] + A[1, 3] A[2, 1]] [A[2, 1] A[3, 2] - A[2, 2] A[3, 1] , -A[1, 1] A[3, 2] + A[1, 2] A[3, 1] , A[1, 1] A[2, 2] - A[1, 2] A[2, 1]] > B := matrix(4,4); B := array(1 .. 4, 1 .. 4, []) > det(B); B[1, 1] B[2, 2] B[3, 3] B[4, 4] - B[1, 1] B[2, 2] B[3, 4] B[4, 3] - B[1, 1] B[3, 2] B[2, 3] B[4, 4] + B[1, 1] B[3, 2] B[2, 4] B[4, 3] + B[1, 1] B[4, 2] B[2, 3] B[3, 4] - B[1, 1] B[4, 2] B[2, 4] B[3, 3] - B[2, 1] B[1, 2] B[3, 3] B[4, 4] + B[2, 1] B[1, 2] B[3, 4] B[4, 3] + B[2, 1] B[3, 2] B[1, 3] B[4, 4] - B[2, 1] B[3, 2] B[1, 4] B[4, 3] - B[2, 1] B[4, 2] B[1, 3] B[3, 4] + B[2, 1] B[4, 2] B[1, 4] B[3, 3] + B[3, 1] B[1, 2] B[2, 3] B[4, 4] - B[3, 1] B[1, 2] B[2, 4] B[4, 3] - B[3, 1] B[2, 2] B[1, 3] B[4, 4] + B[3, 1] B[2, 2] B[1, 4] B[4, 3] + B[3, 1] B[4, 2] B[1, 3] B[2, 4] - B[3, 1] B[4, 2] B[1, 4] B[2, 3] - B[4, 1] B[1, 2] B[2, 3] B[3, 4] + B[4, 1] B[1, 2] B[2, 4] B[3, 3] + B[4, 1] B[2, 2] B[1, 3] B[3, 4] - B[4, 1] B[2, 2] B[1, 4] B[3, 3] - B[4, 1] B[3, 2] B[1, 3] B[2, 4] + B[4, 1] B[3, 2] B[1, 4] B[2, 3] > fibo := proc(n) > if (n > 2) then RETURN(fibo(n-1) + fibo(n-2)); > else RETURN(1); fi; > end; fibo := proc(n) if 2 < n then RETURN(fibo(n - 1) + fibo(n - 2)) else RETURN(1) fi end > for i from 1 to 10 do print(fibo(i)); od; 1 1 2 3 5 8 13 21 34 55 > fibo(25); 75025 > fibo(50); Warning, computation interrupted > combinat[fibonacci](50);