Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
93
Добавлен:
17.04.2013
Размер:
6.14 Кб
Скачать

%Задаем матрицы

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

Соседние файлы в папке other