n=1280;
L=2;
h=L/n;
ys = initclampedmatrix(n)\sinusoidalforces(n);
xs = h:h:L;
plot(xs,ys,'b');
axis([0 2 -.5 0.5]);
title('Problem 7: Deflection, sinusoidal clamped-clamped');
xlabel('x (meters)');
ylabel('y (meters)');