clear all
clc

%% Zadání
F=1;
q=4;
a=0.5;

%% Proměnné x ve třech intervalech:
x1=0:0.01:2*a;
x2=0:0.01:a;
x3=a:0.01:2*a;

%% Na cviceni vypoctene prubehy VSU:
T1=(1/4)*F+(3/2)*q*a-q.*x1;
T2(1:length(x2))=-(3/4)*F-0.5*q*a;
T3(1:length(x3))=(1/4)*F-0.5*q*a;

M1=((1/4)*F+(3/2)*q*a).*x1-0.5*q.*x1.^2;
M2=((3/4)*F+0.5*q*a).*x2;
M3=((3/4)*F+0.5*q*a).*x3-F.*(x3-a);



%% Vykresleni
subplot(2,1,1)
plot(x1,T1,'Color','b','LineWidth',2)
hold on
plot(3*a+flip(x2),T2,'Color','b','LineWidth',2)
plot(2*a+flip(x2),T3,'Color','b','LineWidth',2)
xline(0)
yline(0)
ylabel('T(x) (N)')

subplot(2,1,2)
plot(x1,M1,'Color','b','LineWidth',2)
hold on
plot(3*a+flip(x2),M2,'Color','b','LineWidth',2)
plot(2*a+flip(x2),M3,'Color','b','LineWidth',2)
xline(0)
yline(0)
xlabel('x (m)')
ylabel('M(x) (N.m)')
movegui('center')