clear
flex_a = 1   //f
flex_b = 2    //f
//removing lower support and solving FBD
e = -2 -(2+1)//fp
//e_1 = (2+1+1)*R
//e_1 = -e Making the elongations zero since the both ends are fixed
R = -e/(2+1+1.0) //P
//since sum of forces are 0
printf("\n The reactions at bottom is %0.3f P",R)