#! /usr/bin/calc read 'SistemaLinearBib.calc' print 'Precisão? rode configura()'; apeteco(); A = mat[4,4] = { 3, 2, 0, 0, 1,-3, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1 }; SA = mat[4,9] = { 3, 2, 0, 0, 1, 1,0,0,0, 1,-3, 0, 0, 2, 0,1,0,0, 0, 0, 0, 0, 3, 0,0,1,0, 0, 0, 1, -1, 4, 0,0,0,1 }; Q = mat[4,4] = {1,0,0,0, 0,1,0,0, 0,0,1,0, 0,0,0,1 }; B = mat[4,4]; B = A; coQ = mat[4,4] = {1,0,0,0, 0,1,0,0, 0,0,1,0, 0,0,0,1 } ; Qinv = mat[4,4]; Qinv = Q; PrintMat(SA,4,9);print '--------------'; apeteco(); PrintMat(A,4,4);print '--------------';apeteco(); PrintMat(B,4,4);print '--------------';apeteco(); PrintMat(Q,4,4);print '--------------';apeteco(); PrintMat(Qinv,4,4);print '--------------';apeteco(); PrintMat(Q*Qinv,4,4);print '--------------';apeteco(); print det(A): print '--------------'; w =1; x_a = mat[4,1] = {7/11, -5/11, 0, -4 }; z = 1; x_ah = mat[4,1] = {0, 0, w, w }; print "verificando a homogênea "; PrintMat(A*x_ah ,4,1); apeteco(); print "verificando a não homogênea "; PrintMat(A*x_a ,4,1); apeteco(); print "verificando a solução geral "; w = -5;t=0; for(t=0; t<10; t++) { print t, 'w+t = ',w+t; x_a = {7/11, -5/11, w, -4+w }; PrintMat(A*x_a, 4,1); print '----------------------'; apeteco(); } apeteco();