%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% boundary_conditions.m
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function u = boundary_conditions(x,y)

J = length(x) - 1;
K = length(y) - 1;

u = zeros(J+1,K+1);
  
for k = 1:K+1
    
    u(1,k) = sin(pi*y(k));
    u(J+1,k) = exp(pi)*sin(pi*y(k)) + 0.5*y(k)^2;
    
    for j = 1:J+1
        u(j,1) = 0.0;
        u(j,K+1) = 0.5*x(j)^2;
    end
end