%------------------------------
%          INITIAL.m
%------------------------------
% Must specify non-default 
% values (all zeros) for node n:
%------------------------------
%                w    dw/dx
%------------------------------
% NPBC(n) = 0    U      U
%         = 1    K      U
%         = 2    U      K
%         = 3    K      K
%------------------------------
% RHS(2n-1) = nodal force 
% RHS(2n)   = nodal moment
%------------------------------
% LHS(2n-1) = nodal deflection
% LHS(2n)   = nodal dw/dx
%------------------------------


%----------------------------------------
%   Define test problem:
%----------------------------------------
    Prb = 1  % <- User's choice
%----------------------------------------
%   Prb  1    Cantilever beam with k = 0
%   Prb  2    Foundation support only
%----------------------------------------
%   Note, this selection will also be 
%   used in COEF.m
%----------------------------------------


%-----------------
% Note: NUMNP = 4
%-----------------

if Prb == 1
  NPBC(1)          =  3;  % Left node
  NPBC(4)          =  0;  % Right node
  RHS(2*4-1)       = -1;  % Force on right end

elseif Prb == 2
  NPBC(1)      =  0;  % Left node
  NPBC(4)      =  0;  % Right node

end 
