## Sets up the initial profile as a Gaussian 
## Usage : function y=initial_profile(x,x0,k)
## x is the array that contains values on the horizontal axis, x0 is the
## center of the Gaussian and k is a parameter that determines the width
## of Gaussian

function y=initial_profile(x,x0,k)
  
  y=exp(-k*(x-x0).^2);

endfunction
