Indeterminate Beam Analysis, how to create integrated piecewise function?

Hey guys, I'm having trouble creating shear and moment diagrams for the problem in the link. The following is my code : clc; clear;

E = 210; %GPA I = 210e-4; % MoI % KN/m x1= 0:1:5; V1 = -5*x1 + 1.5; %%KN ||| 0<x<5

x2 = 5:1:9; V2 = -5*x2 + 29.5; % 5<x<9

xt = [x1 x2]; zt = [V1 V2];

subplot(3,1,1) plot (xt,zt,'linewidth', 2),title('shear Diagram'),... xlabel('Position along the Beam, Ft'), ylabel('Shear Force, Lbs')

syms x y v1 = -5*x + 1.5; v2 = -5*y + 29.5;

f = int(v1); x = 0:1:5; F = subs(f);

f1 = int(v2); y = 5:1:9; F1 = subs(f1);

C = [x y]; F2 = [F F1];

subplot(3,1,2) plot (C,F2,'linewidth',2)

My thought process was as follows - i found the shear forces on the beam ( indeterminate beam theory method) and then by definition, use the integration operator in matlab to get the moment equation, symbolically. Now, I believe the shear force part of the code is right, but the discontinuity at x = 5 for the 2nd graph isn't supposed to be there. It's supposed to be a parabola that goes to zero. I think I'm either not defining that constraint. Any help would be appreciated, thank you!