%Program to solve toy beam bending 
%problem

%Calculate reactions
A = [1 1; 0 12];
B = [400; 2400];
R = A\B; %Reaction vector

%Calculate sigma_x
M = -600;
ro = 2e-2;
ri = 1e-2;
I = pi*(ro^4 - ri^4)/4;
sigma_x = 1e-6*M*ro/I