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

UptoLike

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

Рубрика:

if exist('g4')==0

g4='10*sin(6*x)';

end

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

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

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

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

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

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

G1=inline(g1,'y');

G2=inline(g2,'y');

G3=inline(g3,'x');

G4=inline(g4,'x');

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

N=n*m;

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

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

a=zeros(N,N);

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

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

b=zeros(1,N);

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

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

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

for j=1:m

b(j)=G1(y(j));

if v1==1

a(j,j)=1;

elseif v1==2

a(j,j)=-1/dx;

a(j,m+j)=1/dx;

else

error('Parameter v1 have incorrect value');