Решение задач оптимального управления с использованием математической системы MATLAB и пакета имитационного моделирования SIMULINK. Сивохин А.В - 58 стр.

  % Вычисление нормы погрешности приближенного решения (u является
  % вектор-столбцом, а подфункция exsol вызывается от строк, поэтому
  % вектор с точным решением необходимо транспонировать)
  erhelp = norm(u-uex',Inf);
end
% Решение найдено с требуемой точностью
% Визуализация расчетной триангуляции и вывод контурного графика
% решения, залитого цветом
figure
subplot(2,1,1)
pdemesh(p,e,t)
subplot(2,1,2)
pdeplot(p, e, t, 'xydata', u, 'colormap', 'gray', 'colorbar', 'off')

function z = exsol(x,y)
% Подфункция для вычисления точного решения
z = sin(pi*x).*sin(2*pi*y);

    Файл-программа для решения задачи о составлении рациона:

% Задание матрицы и вектора правой части системы неравенств
A = [4 6 15
   2 2 0
   5 3 4
   7 3 12];
A = -A;
b = [250; 60; 100; 220];
b = -b;
% Определение коэффициентов целевой функции
f = [44; 35; 100];
% Задание ограничений снизу на переменные
lb =[0; 0; 0];
% Решение и вывод результата в командное окно
x = linprog(f, A, b, [], [], lb, [])


            Файл-функция, вычисляющая минимизируемую функцию:

function f = myfun(x)
f = 3*x(1)^2+2*x(2)^2;

    Файл-функция с ограничениями:

function [c, ceq] = mycon(x)