format long
disp('Calculated solution at n=10:');
ys = initmatrix(10)\simpleforces(10)