%---------------------------------------------------
%    INITIAL.m
%
%    Set problem parameters when not default values.
%
%    NPBC   the boundary condition specification.
%    PHI    when specified as known.
%    Q      point sources.
%    QS     nodal values for flux.
%    HS     coeficient for convective heat transfer.
%    PHIa   nodal values of ambient temperature.
%
%-----------------------------------------------------

%--------------------------------------
%     This is written for the test
%     problem for which the exact
%     solution is:
%         PHI = 20 - 2x + 3y
%     Only node 9 is left unspecified.
%--------------------------------------
      for Ix=1:NUMNP
        if NPcode(Ix) == 1
	   NPBC(Ix)=1;
	   PHI(Ix) = 20-2*XORD(Ix)+3*YORD(Ix);
	else
	   NPBC(Ix)=0;
	end
      end

%---------------------------------------
%    Print out exact answer for node 9
%---------------------------------------
      exact = 20 - 2*XORD(9) + 3*YORD(9)
%---------------------------------------
