plane_truss: 1 based indexing at user interface
fem2d: revamp of linear solver.  call update_state now
       to update stiffness.  stress is computed at end
       of step.