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

UptoLike

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

111
if exist('vt2')==0
vt2=1;
end
if exist('gt2')==0
gt2='sin(4*y)-sin(4*x)';
end
if exist('v1')==0
v1=2;
end
if exist('g1')==0
g1='0';
end
if exist('v2')==0
v2=2;
end
if exist('g2')==0
g2='0';
end
if exist('v3')==0
v3=2;
end
if exist('g3')==0
g3='0';
end
if exist('v4')==0
v4=2;
end
if exist('g4')==0
g4='0';
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);
t=t0:(ts-t0)/(s-1):ts; dt=t(2)-t(1);
% Вычисление значений функций, заданных символьно,
% в узлах координатной сетки
GT1=inline(gt1,'x','y');
GT2=inline(gt2,'x','y');
G1=inline(g1,'y');
G2=inline(g2,'y');
G3=inline(g3,'x');