Методы решения задач математической физики. Рындин Е.А. - 70 стр.

UptoLike

Составители: 

70
% Задание матрицы-строки свободных членов СЛАУ
% размерности 1 x n,
% все элементы которой равны 0
b=zeros(1,n);
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих граничным условиям и проверка
% корректности значений параметров v1, v2
b(1)=G1(x(1));
if v1==1
a(1,1)=1;
elseif v1==2
a(1,1)=-1/dx;
a(1,2)=1/dx;
else
error('Parameter v1 have incorrect value');
end
b(n)=G2(x(n));
if v2==1
a(n,n)=1;
elseif v2==2
a(n,n)=1/dx;
a(n,n-1)=-1/dx;
else
error('Parameter v2 have incorrect value');
end
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих внутренним точкам области
for i=2:n-1
a(i,i)=-2/dx^2;
a(i,i+1)=1/dx^2;
a(i,i-1)=1/dx^2;
b(i)=F(x(i));
end
% Решение СЛАУ
u=b/a';