x=linspace(0,10,1000);
uy=-5/6*(sing(x,4,0)-sing(x,4,5))+15/16*sing(x,3,8)+75*sing(x,2,7)+(57/6)*x.^3-238.25*x;
figure;
plot(x,uy);
title('Beam Displacement');
xlabel('x displacement (feet)');
ylabel('Uy (millifeet)');