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

UptoLike

ВУЗ: 

ТТИ ЮФУ | Таганрог

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

Рубрика: 

112
G4=inline(g4,'x');
% Определение размерности СЛАУ
N=s*n*m;
% Задание матрицы коэффициентов СЛАУ размерности N x N,
% все элементы которой равны 0
a=zeros(N,N);
% Задание матрицы-строки свободных членов СЛАУ размерно-
сти 1 x N,
% все элементы которой равны 0
b=zeros(1,N);
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих начальным и граничным условиям
% и проверка корректности
% значений параметров vt1, vt2, v1, v2, v3, v4
for i=1:n
for j=1:m
b(m*(i-1)+j)=GT1(x(i),y(j));
if vt1==1
a(m*(i-1)+j,m*(i-1)+j)=1;
elseif vt1==2
a(m*(i-1)+j,m*(i-1)+j)=-1/dt;
a(m*(i-1)+j,n*m+m*(i-1)+j)=1/dt;
else
error('Parameter vt1 have incorrect value');
end
b(n*m*(s-1)+m*(i-1)+j)=GT2(x(i),y(j));
if vt2==1
a(n*m*(s-1)+m*(i-1)+j,n*m*(s-1)+m*(i-1)+j)=1;
elseif vt2==2
a(n*m*(s-1)+m*(i-1)+j,n*m*(s-1)+m*(i-1)+j)=1/dt;
a(n*m*(s-1)+m*(i-1)+j,n*m*(s-2)+m*(i-1)+j)= ...
-1/dt;
else
error('Parameter vt2 have incorrect value');
end
end
end

Страницы