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

UptoLike

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

Рубрика:

% Задание равномерной сетки

x=x0:(xn-x0)/(n-1):xn; dx=x(2)-x(1);

y=y0:(ym-y0)/(m-1):ym; dy=y(2)-y(1);

t=t0:(ts-t0)/(s-1):ts; dt=t(2)-t(1);

% Вычисление значений функций, заданных символьно,

% в узлах координатной сетки

F=inline(f,'x','y');

R=inline(r,'x','y');

C=inline(c,'x','y');

K=inline(k,'x','y');

GT=inline(gt1,'x','y');

G1=inline(g1,'y');

G2=inline(g2,'y');

G3=inline(g3,'x');

G4=inline(g4,'x');

% Определение размерности СЛАУ

N=s*n*m;

% Задание матрицы коэффициентов СЛАУ размерности N x N,

% все элементы которой равны 0

a=zeros(N,N);

% Задание матрицы-строки свободных членов СЛАУ

% размерности 1 x N, все элементы которой равны 0

b=zeros(1,N);

% Определение коэффициентов и свободных членов СЛАУ,

% соответствующих граничным условиям и проверка

% корректности значений параметров vt, v1, v2, v3, v4

for i=1:n

for j=1:m

b(m*(i-1)+j)=GT(x(i),y(j));

if vt==1

a(m*(i-1)+j,m*(i-1)+j)=1;

elseif vt==2

a(m*(i-1)+j,m*(i-1)+j)=-1/dt;