%------------------------------
%          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 = 2  % <- 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)          =  1;  % Left node
  NPBC(NUMNP)      =  1;  % Right node
  LHS(1)           = -1;  % deflection on left end
  LHS(2*NUMNP-1)   = 1;   % deflection on right end

elseif Prb == 2
   NPBC(1)          =  1;  % Left node
   NPBC(NUMNP)      =  1;  % Right node
  LHS(1)           = -1;  % deflection on left end
  LHS(2*NUMNP-1)   = 1;   % deflection on right end

end 
