---

clc;
clear;
close all;

% Вариант 14
% J = integral(0.002*t^3 + x' + x^2)dt
% x(1)=5
% x(2)=1


% Начальное приближение
solinit = bvpinit(linspace(1,2,20), [5 -1]);


% Решение краевой задачи
sol = bvp4c(@ode14, @bc14, solinit);


% Построение траектории
t = linspace(1,2,100);

x = deval(sol,t);


plot(t,x(1,:),'LineWidth',2)

grid on
xlabel('t')
ylabel('x(t)')
title('Оптимальная траектория x(t), вариант 14')


% =========================
% Система дифференциальных
% уравнений Эйлера
% =========================

function dx = ode14(t,x)

% x(1)=x
% x(2)=dx/dt

dx = [
    x(2)
    0
];

end


% =========================
% Граничные условия
% =========================

function res = bc14(x0,x1)

res = [
    x0(1)-5
    x1(1)-1
];

end