function new_images = gaussian_smooth(images, sigma)
% AVERAGE_SMOOTH: Given sequence of 1D data and the value sigma,
%                 Gaussian smooth will be applied.
%
% Code written by Katherine Smith, 2003
%
%    GENERAL
%
%      new_images = gaussian_smooth(images, sigma)
%
%    INPUT/S
%
%      -images:
%           Input images.
%
%      -sigma:
%           The sigma value (variance) of the Gaussian.
%           
%    OUTPUT/S
%
%      -new_images:
%           The images with Gaussian smooth applied.
%
%    PENDING WORK
%
%      -
%
%    KNOWN BUG/S
%
%      -
%
%    COMMENT/S
%
%      -This is an intermediate routine that uses Gaussian filter
%
%    RELATED FUNCTION/S
%
%      AVERAGE_SMOOTH
%
%    ABOUT
%
%      -Created:     November 23rd, 2003
%      -Last update: Novermber 25th, 2003
%      -Revision:    0.0.2
%      -Author:      R. S. Schestowitz, University of Manchester
% ==============================================================

filter = fspecial('gaussian',[9 1], sigma); %pass on to built-in function
new_images = imfilter(images,filter); %apply the filter set above