## Copyright (C) 2017 Marina Campolo
## 
## This program is free software; you can redistribute it and/or modify it
## under the terms of the GNU General Public License as published by
## the Free Software Foundation; either version 3 of the License, or
## (at your option) any later version.
## 
## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
## GNU General Public License for more details.
## 
## You should have received a copy of the GNU General Public License
## along with this program.  If not, see <http://www.gnu.org/licenses/>.

## -*- texinfo -*- 
## @deftypefn {Function File} {@var{retval} =} conc_unbounded (@var{input1}, @var{input2})
##
## @seealso{}
## @end deftypefn

## Author: Marina Campolo <marina@boopsie>
## Created: 2017-03-27

function [c_u] = conc_unbounded (Nx,Ny,x,ys)
global M w ;
global u_x e_y;
for i=1:Nx;
xx=x(i);
for j=1:Ny;
yys=ys(j);
c_u(j,i)=M/w*1./(4.*pi*e_y*xx*u_x)^0.5*exp(-u_x*yys^2/(4*e_y*xx));
endfor
endfor
endfunction