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