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

    D2X = subs(D2X,{t},{0:0.1:1})
    t = 0:0.1:1;
%-- 5.Визуализация U(t):
    U = exp(4.*t);
    subplot(4,2,1)
    plot(t,U,'r')
    xlabel('t')
    ylabel('U')
%-- 6.Визуализация X(t):
    subplot(4,2,3)
    plot(t,X,'r')
    xlabel('t')
    ylabel('X')
%-- 7.Визуализация DX(t):
    subplot(4,2,2)
    plot(t,DX,'r')
    xlabel('t')
    ylabel('DX')
%-- 8.Визуализация D2X(t):
    subplot(4,2,4)
    plot(t,D2X,'r')
    xlabel('t')
    ylabel('D2X')
%-- 9.Конец функции DSolveXDXD2X
function [X,DX,D2X] = EvalXDXD2X
%-- ФУНКЦИЯ ДЛЯ ВЫЧИСЕНИЯ ЗНАЧЕНИЙ U(t),X(t) =
%-- 1/5e(4t)+4/5e(-t) И ЕЕ ПРОИЗВОДНЫХ:
%-- 1.Задание вектора t:
    t = 0:0.1:1;
%-- 2.Вычисление U(t):
    U = exp(4.*t);
%-- 3.Вычисление X(t):
    X = 1./5.*exp(4.*t)+4./5.*exp(-t);
%-- 4.Вычисление DX(t):
    DX = 4./5.*exp(4.*t)-4./5.*exp(-t);
%-- 5.Вычисление D2X(t):
    D2X = 16./5.*exp(4.*t)+4./5.*exp(-t);
%-- 6.Визуализация U(t):
    subplot(4,2,1)
    plot(t,U,'r')
    xlabel('t')
    ylabel('U')
%-- 7.Визуализация X(t):
    subplot(4,2,3)
    plot(t,X,'r')