Добавил:
korayakov
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:ЧМ (ЭКТ-3) / Лабы / Всякие сделанные варианты / other / laba3p1
.doc
%Задаем матрицы
A=rand(1,6)
B=rand(1,7)
C=rand(1,6)
F=rand(1,7)
al=0
be=0
%Задаем нулевые
al(1)=-C(1)/B(1)
be(1)=F(1)/B(1)
for i=2:size(F,2)-1
al(i)=-C(i)/((A(i-1)*al(i-1))+B(i))
be(i)=(F(i)-(A(i-1)*be(i-1)))/((A(i-1)*al(i-1))+B(i))
end
y=0
y(size(F,2),1)=(F(size(F,2))-(A(size(A,2))*be(size(be,2))))/((A(size(A,2))*al(size(al,2)))+B(size(B,2)))
for i=size(F,2)-1:-1:1
y(i,1)=al(i)*y(i+1,1)+be(i)
end
D=0
for i=1:size(B,2)
D(i,i)=B(i)
end
for i=1:size(A,2)
D(i+1,i)=A(i)
end
for i=1:size(C,2)
D(i,i+1)=C(i)
end
