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