...
Now we will need to call the function bending_press
we just created to obtain the value for σx that is needed to create the plot. Do this by typing:
Run the file beam4.m and check the output You should get the same plot you got in the file beam3.m
...