just to pitch in, i think the discrepancy between the exact solution and the result from FEA is due mainly to the very idea that FEA is an approximate method of solving diferrential equations by discreteizing the problem domain and by using some functions to interpolate the results at gauss...