function femsol
% finite element solution, updates global jmf
%
% input fields of global jmf:
%
% node(:,i)     coordinates of node i
% ien(a,e)      global node number of node a on element e
% id(i,ga)      global equation number for dof i on global node a
% ibdry(a,b)    global node number of node a on boundary element b
% bload(:,b)    load on boundary element i (assumed constant)
% vload(:,e)    load on element e (assumed constant)
% alternatively bload and vload may be functions of coordinates 
%
% output fields
%
% Kloc(:,:,e)   local stiffness matrix of element e

global jmf % updated