clc,clear all
A=[8 -2 0 0; -1 6 -2 0 ;0 2 10 -4 ; 0 0 -1 6]
b=[6 3 8 5]'
x = A\b
[n,n]=size(A);

P(1)=1;
Q(1)=1;
for i=2:n-1
zn=-A(i,i)-A(i,i-1)*P(i-1);
P(i)=A(i,i+1)/zn;
Q(i)=(A(i,i-1)*Q(i-1)-b(i))/zn;
end
Q(n)=(A(n,n-1)*Q(n-1)-b(n-1))/A(n,n)-A(n,n-1)*P(n-1);
x(n)=Q(n);
for i=n-1:-1:1
x(i)=P(i)*x(i+1)+Q(i);
end;





Чтобы не видеть здесь видео-рекламу достаточно стать зарегистрированным пользователем.
Чтобы не видеть никакую рекламу на сайте, нужно стать VIP-пользователем.
Это можно сделать совершенно бесплатно. Читайте подробности тут.