xs = (1:10)/5;
ys = initmatrix(10)\simpleforces(10);
correct = simplecorrect(10);

format long
plot(xs, [ys, correct]);
axis([0.2 2 -1 0.5]);
title('Problem 2: Deflection, Simple Case');
xlabel('x (meters)');
ylabel('y (meters)');
disp('Error at x=2m: ');
disp(ys(10)-correct(10));