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

UptoLike

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

71
% Построение графика искомой функции u(x,y)
plot(x,u)
xlabel('x')
ylabel('u')
grid on
Строки, начинающиеся символом %, являются комментариями. Приведен-
ное описание необходимо сохранить в виде текстового файла с именем
puass_1d.m и поместить его в каталог /WORK, находящийся в корневом катало-
ге системы MATLAB.
Вызов функции puass_1d может осуществляться следующими командами:
puass_1d;
x=puass_1d;
[x,u]=puass_1d;
[x,u]=puass_1d(x0,xn,n,f,v1,g1,v2,g2).
При использовании первой, второй или третьей команд функция будет вы-
водить графически решение задачи при входных параметрах, принятых по
умолчанию.
При вызове функции с помощью первой или второй команд без символа ;
после текста команды, на экран монитора помимо графического представления
решения будет выводиться вектор координат узлов сетки х. При использовании
символа ; вектор координат узлов сетки выводиться не будет.
При использовании третьей или четвертой команд без символа ; после тек-
ста команды, на экран монитора будет выводиться график решения, вектор ко-
ординат узлов сетки х и вектор значений искомой функции в узлах сетки u.
Например, при вызове функции командой
puass_1d;
на экране появится график решения задачи для функции f(х) в правой части
уравнения (235), определяемой выражением
)cos(10)sin(2
)(
22
xx
xf
+=
, (246)
с граничными условиями
0
)(
min
=
x
u
; (247)
5.0
max
=
dx
du
x
(248)
на границах x
min
= 0, x
max
= 5 на равномерной координатной сетке, содержащей
50 точек (рис. 21).