beam2.f90 running E= 10.0000000000000000 I= 9.0000000000000000 W= 7.0000000000000000 L= 8.0000000000000000 h= 1.0000000000000000 n= 9 Analytic Solution. beam x= 0.0000000000000000E+000 y= -0.0000000000000000E+000 slope= 0.0000000000000000E+000 beam x= 1.0000000000000000 y= -1.9849537037037037E-02 slope= -3.4027777777777775E-02 beam x= 2.0000000000000000 y= -5.8333333333333334E-02 slope= -3.8888888888888896E-02 beam x= 3.0000000000000000 y= -9.1145833333333329E-02 slope= -2.4305555555555556E-02 beam x= 4.0000000000000000 y= -0.1037037037037037 slope= -3.9474596431116678E-17 beam x= 5.0000000000000000 y= -9.1145833333333329E-02 slope= 2.4305555555555556E-02 beam x= 6.0000000000000000 y= -5.8333333333333334E-02 slope= 3.8888888888888890E-02 beam x= 7.0000000000000000 y= -1.9849537037037037E-02 slope= 3.4027777777777775E-02 beam x= 8.0000000000000000 y= -0.0000000000000000E+000 slope= -7.8949192862233357E-17 Numerical Solution. Fourth order difference equation solution boundary dy/dx=0 fourth order n= 7, h=.10000E+01, C= 1.00 -4.00 6.00 -4.00 1.00 solve A y = F for y A F 48.000 -36.000 16.000 -3.000 0.000 0.000 0.000 0.0000000000000000E+000 -4.000 6.000 -4.000 1.000 0.000 0.000 0.000 -9.7222222222222224E-03 1.000 -4.000 6.000 -4.000 1.000 0.000 0.000 -9.7222222222222224E-03 0.000 1.000 -4.000 6.000 -4.000 1.000 0.000 -9.7222222222222224E-03 0.000 0.000 1.000 -4.000 6.000 -4.000 1.000 -9.7222222222222224E-03 0.000 0.000 0.000 1.000 -4.000 6.000 -4.000 -9.7222222222222224E-03 0.000 0.000 0.000 -3.000 16.000 -36.000 48.000 0.0000000000000000E+000 fourth order difference equations are exact within numerical accuracy when solution is fourth order. x=0.0 y=0.0 boundary x= 1.0000000000000000 y= -1.9849537037036919E-02 x= 2.0000000000000000 y= -5.8333333333332890E-02 x= 3.0000000000000000 y= -9.1145833333332621E-02 x= 4.0000000000000000 y= -0.1037037037037029 x= 5.0000000000000000 y= -9.1145833333332649E-02 x= 6.0000000000000000 y= -5.8333333333332911E-02 x= 7.0000000000000000 y= -1.9849537037036909E-02 x= 8.0000000000000000 y=0.0 boundary