Matematické Fórum

Zdravím,
mám takový problém. Měl bych vyřešit tento příklad , ale poté co napíši  nějaký program na GE tak mi to běží hodinu a to se strašně velkou chybou výpočtu takže to budu mít nejspíše celé blbě :D proto bych se rád zeptal zda někdo nevíte jak na to.
Matici mám nadefinovanou tímto způsobem:

function [ x ] = eliminace( v )
tic
%nadefinovani matice
A=spalloc(v,v,3*v-2)
for i=1:v
A(i,i)=4;
end
for i=1:v-1
A(i+1,i)=1;
A(i,v)=1;
end
disp('matice')
%nadefinovani vektoru
b=rand(v,1)
[m,n]=size(A);
disp('vektor')
%pro chybu

a pro výpočet chyby používám:

po nadefinování matic
AA=A;
bb=b;
a po skončení programu
Chyba = norm(AA*x-bb);
    Chyba

za jakoukoliv pomoc s programem bych byl moc vděčný.