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

UptoLike

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

79
end
b(m*(n-1)+j)=G2(y(j));
if v2==1
a(m*(n-1)+j,m*(n-1)+j)=1;
elseif v2==2
a(m*(n-1)+j,m*(n-1)+j)=1/dx;
a(m*(n-1)+j,m*(n-2)+j)=-1/dx;
else
error('Parameter v2 have incorrect value');
end
end
for i=2:n-1
b(m*(i-1)+1)=G3(x(i));
if v3==1
a(m*(i-1)+1,m*(i-1)+1)=1;
elseif v3==2
a(m*(i-1)+1,m*(i-1)+1)=-1/dy;
a(m*(i-1)+1,m*(i-1)+2)=1/dy;
else
error('Parameter v3 have incorrect value');
end
b(m*(i-1)+m)=G4(x(i));
if v4==1
a(m*(i-1)+m,m*(i-1)+m)=1;
elseif v4==2
a(m*(i-1)+m,m*(i-1)+m)=1/dy;
a(m*(i-1)+m,m*(i-1)+m-1)=-1/dy;
else
error('Parameter v4 have incorrect value');
end
end
% Определение коэффициентов и свободных членов СЛАУ,
% соответствующих внутренним точкам области
for i=2:n-1
for j=2:m-1
a(m*(i-1)+j,m*(i-1)+j)=-2/dx^2-2/dy^2;
a(m*(i-1)+j,m*(i)+j)=1/dx^2;
a(m*(i-1)+j,m*(i-2)+j)=1/dx^2;
a(m*(i-1)+j,m*(i-1)+j+1)=1/dy^2;
a(m*(i-1)+j,m*(i-1)+j-1)=1/dy^2;
b(m*(i-1)+j)=F(x(i),y(j));
end
end