Xlabel('X')

ylabel('y')

zlabel('z(x,y)')

v = axis;

da = daspect;

da(1:2) = min(da(1:2));

daspect(da);

axis(v);

colormap(gray) % палитра

Рис. 12.2. Решение примера 12b

Ответ. Для сетки 2040 система уравнений МКР вида (12.3) имеет вид

(12.0)

при этом входящие в (12.8) величины z при i=1, i=20, k=1, k=40 находятся из граничного условия в соответствующей точке. График экстремальной поверхности показан на рисунке в зоне вывода. Видно, что он очень похож на поверхность, построенной по МКЭ в главе 4.

3. Задание

Решить примеры 1a и 4 методом конечных разностей. Сравнить решения с аналитическими.

8