// questao4.sce

// importando as funcoes necessarias

getf("norma.sci");
getf("seidel_decomp.sci");
getf("inverte.sci");

// declarando as matrizes

A1 = [9,1,1,1; 1,8,1,1; 1,1,7,1; 1,1,1,6];
A2 = [2,1,0,0; 1,2,1,0; 0,1,2,1; 0,0,1,2];

// computando a decomposicao N - P de cada uma

[N1, P1] = seidel_decomp(A1);
[N2, P2] = seidel_decomp(A2);

// obtendo, em cada caso, a norma de N^{-1}P

m1 = norma(inverte(N1) * P1);
m2 = norma(inverte(N2) * P2);

// RESULTADO:
// m1 = 0.3333333
// m2 = 0.875
