Метод Прогонки (Алгоритм Томаса) в Матлаб (MatLab)

[Саша Дмитриева]

Здравствуйте. У меня такая вот пролема: НЕ МОГУ НАЙТИ ОШИБКУ В ПРОГРАММЕ. ПОМОГИТЕ МНЕ, ПОЖАЛУЙСТА ЕЕ НАЙТИ И ИСПРАВИТЬ? Буду премного благодарна Программа должна осуществлять метод прогонки. Вот теория: http://mathalgo.blogspot.com/2007/11/blog-post.html
Вот алгоритм который прошу исправить:

function x=tomas(A,b);
A= [0.79, -0.12, 0.34, 0.16; -0.34, 1.18, -0.17, 0.18; -0.16, -0.34, 0.85, 0.31; -0.12, 0.26, 0.08, 0.75]
b= [-0.64; 1.42; -0.42; 0.83]
[n,n]=size(A);
z=zeros(n,1);
A=[z,A,z];
P(1)=1;
Q(1)=1;
for i=2:n
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
for i=n-1:-1:1
x(i+1)=1;
x(i)=P(i)*x(i+1)+Q(i);
end

ПОЖАЛУЙСТА!!! =)

[Саша Дмитриева]

Ну если никто не хочет искать ошибку (а я уверена что может на самом деле и именно не хочет), то может хоть кто нибудь напишет алгорит программы в матлаб на метод прогонки? м?

[Саша Дмитриева]

ВОТ ПРАВИЛЬНЫЙ КОД ПРОГРАММЫ!

ПРОВЕРЯЛА ВРУЧНУЮ - АЛГОРИТМ ВЕРНЫЙ.

function x=tomas(A,b);
A= [0.79, -0.12, 0.34, 0.16; -0.34, 1.18, -0.17, 0.18; -0.16, -0.34, 0.85, 0.31; -0.12, 0.26, 0.08, 0.75]
b= [-0.64; 1.42; -0.42; 0.83]
[n,n]=size(A);
%z=zeros(n,1);
%A=[z,A,z];
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

с любовью ваша Сашка!!!

[Саша Дмитриева]

ХОТЬ БЫ КТО ПОХВАЛИЛ МЕНЯ УМНИЦУ! )

[Sonny01]

молодец

[alisan]

А может кто-нибудь написать эту же программу, только в C#? Нигде не могу найти данный метод на нужном языке. SOS!

[adam999]

САша...МНЕ тож нужна та же программа..метод программы.....если есть поставь да всю прогу

[IRIKA]

Саша, у меня тоже проблемы с методом прогонки(((
можешь помочь?