Matlab ========================================
>> Plane_Sress_RS    % by numerical integration
Read 5 mesh coordinate pairs 
Read 4 elements with 3 nodes each 

X_disp    Y_disp    at 5 nodes 
   1.0e-04 *
         0    0.3403
         0         0
    0.5243    0.1701
    0.6667    0.1667
    0.6667    0.1736

Node, Component, Value, Eq. for 3 Reactions 
1 1 -5000 1 
2 1 -5000 3 
2 2 -9.37916e-13 4 

Elem, X_qp, Y_qp, Sig_x, Sig_y, Sig_xy 
1 0.333333 1 770833 -62500 1.01644e-10 
2 1 0.333333 500000 -5.82077e-11 62500 
3 1.66667 1 229167 62500 -2.03288e-11 
4 1 1.66667 500000 -2.91038e-11 -62500 

>> Plane_Sress_XY   % closed form arrays
Read 5 mesh coordinate pairs 
Read 4 elements with 3 nodes each 

X_disp    Y_disp    at 5 nodes 
   1.0e-04 *
         0    0.3403
         0         0
    0.5243    0.1701
    0.6667    0.1667
    0.6667    0.1736

Node, Component, Value, Eq. for 3 Reactions 
1 1 -5000 1 
2 1 -5000 3 
2 2 -9.37916e-13 4 

Elem, X_qp, Y_qp, Sig_x, Sig_y, Sig_xy 
1 0.333333 1 770833 -62500 1.01644e-10 
2 1 0.333333 500000 -5.82077e-11 62500 
3 1.66667 1 229167 62500 -2.03288e-11 
4 1 1.66667 500000 -2.91038e-11 -62500 
>> quit

UNIX ========================================
% more msh_bc_xyz.tmp
  10  0.00000E+00  0.00000E+00
  11  0.00000E+00  2.00000E+00
  00  1.00000E+00  1.00000E+00
  00  2.00000E+00  0.00000E+00
  00  2.00000E+00  2.00000E+00

% more msh_typ_nodes.tmp
    1    1    3    2
    1    1    4    3
    1    4    5    3
    1    3    5    2

% more node_results.tmp
0 3.40278e-05 
0 0 
5.24306e-05 1.70139e-05 
6.66667e-05 1.66667e-05 
6.66667e-05 1.73611e-05 

% more node_reaction.tmp
1 1 -5000 1 
2 1 -5000 3 
2 2 -9.37916e-13 4 

% more el_qp_xyz_fluxes.tmp
0.333333 1 770833 -62500 1.01644e-10 
1 0.333333 500000 -5.82077e-11 62500 
1.66667 1 229167 62500 -2.03288e-11 
1 1.66667 500000 -2.91038e-11 -62500 

% more MODEL_SOLUTION ===============================
*** REACTION RECOVERY ***
NODE, PARAMETER,     REACTION, EQUATION
   1,  DOF_1,      -5.0000E+03     1
   2,  DOF_1,      -5.0000E+03     3
   2,  DOF_2,      -8.5265E-14     4

*** EXTREME VALUES OF THE NODAL PARAMETERS ***
PARAMETER      MAXIMUM,   NODE      MINIMUM,   NODE
 DOF_1,     6.6667E-05,      4   0.0000E+00,      1
 DOF_2,     3.4028E-05,      1   0.0000E+00,      2

***  OUTPUT OF RESULTS IN NODAL ORDER  ***
NODE,  X-Coord,    Y-Coord,    DOF_1,      DOF_2,     
   1  0.0000E+00  0.0000E+00  0.0000E+00  3.4028E-05 
   2  0.0000E+00  2.0000E+00  0.0000E+00  0.0000E+00 
   3  1.0000E+00  1.0000E+00  5.2431E-05  1.7014E-05 
   4  2.0000E+00  0.0000E+00  6.6667E-05  1.6667E-05 
   5  2.0000E+00  2.0000E+00  6.6667E-05  1.7361E-05 
  
*** STRESSES AT INTEGRATION POINTS ***
           COORDINATES              STRESSES
POINT     X           Y         XX          YY
POINT                           XY          EFFECTIVE
         ELEMENT NUMBER  1
  1  3.333E-01  1.000E+00   7.70833E+05  -6.25000E+04
  1                         4.06576E-11   8.03908E+05
         ELEMENT NUMBER  2
  1  1.000E+00  3.333E-01   5.00000E+05  -5.82077E-11
  1                         6.25000E+04   5.11585E+05
         ELEMENT NUMBER  3
  1  1.667E+00  1.000E+00   2.29167E+05   6.25000E+04
  1                         0.00000E+00   2.05185E+05
         ELEMENT NUMBER  4
  1  1.000E+00  1.667E+00   5.00000E+05   2.91038E-11
  1                        -6.25000E+04   5.11585E+05
 NORMAL END OF MODEL_F90 (NO WARNINGS)