PROGON

n1=input('Введите номер группы: ');
n2=input('Bведите ваш номер в списке группы: ');
n=n1+n2;
rand('seed',n);
m=7;
m1=m-1;
c1=2;
c2=3;
e=ones(1,m);
a=rand(1,m1);
b=rand(1,m)+e*c1;
b=b';
c=rand(1,m1);
d=rand(1,m)+e;
R1=diag(a,-1);
R2=diag(b);
R3=diag(c,1);
disp('Матрица системы:');
A=(R1+R2+R3)*c2
disp('Правая часть системы:');
b
P=size(A);
n=P(1);
alfa(1)=0;
beta(1)=A(1,1);
gamma(1)=A(1,2);
alfa(n)=A(n,n-1);
beta(n)=A(n,n);
gamma(n)=0;
for i=2:n-1
alfa(i)=A(i,i-1);
beta(i)=A(i,i);
gamma(i)=A(i,i+1);
end
c(1)=0;
d(1)=0;
c(2)=-gamma(1)/beta(1);
d(2)=b(1)/beta(1);
for i=3:n
c(i)=-gamma(i-1)/(alfa(i-1)*c(i-1)+beta(i-1));
d(i)=(b(i-1)-alfa(i-1)*d(i-1))/(alfa(i-1)*c(i-1)+beta(i-1));
end
x(n)=(b(n)-alfa(n)*d(n))/(alfa(n)*c(n)+beta(n));
for i=n-1:(-1):1
x(i)=c(i+1)*x(i+1)+d(i+1);
end
disp('ОТВЕТ:');
x'
disp('Проверка, точное решение:');
A^(-1)*b