---

clear;
clc;
close all;

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


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


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


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

x = deval(sol,t);


figure;

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

grid on;

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

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