Exercicio Scilab

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

clc
A=[0.2,1,12;14.6,16.7,17;85.2,82.3,71]
b=[3;16;81]
M =[A b]
disp ("Matriz aumentada")
disp (M)
P=(rref(M))
disp("Transposta de X") 
disp (P')
//Letra B
inversa = inv(A)
disp("Matriz inversa")
disp(inversa)
disp ("Solução X = A-1 b")
sol=[inv(A) b]
disp (sol)
//Letra C
rref(M)
//Metodo de Gauss-Jacob
X0=[1;0.05;5]
eps1=0.0001
R=gaussjacob(A,b,X0,eps1)
//A matriz não é diagonalmente dominante