---

clear;
clc;
close all;

% Вариант 14
% J = integral(0.002*t^3 + x' + x^2)dt
% Граничные условия:
% x(1)=5
% x(2)=1


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


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


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

x = deval(sol,t);


figure;
plot(t,x(1,:),'r','LineWidth',2);

grid on;

xlabel('Время t');
ylabel('Координата x(t)');

title('Оптимальная траектория движения (Вариант 14)');