function Demo_Gaussian_Kernel


DIM = 2;


if DIM == 1
    
    sigma = 0.1;
    
    xi = 0;
    
    x = -1:0.01:1;
    
    for j = 1:length(x)
        k(j) = Gaussian_kernel(xi,x(j),sigma);
    end
    
    figure;
    plot(x,k,'r-');
    xlabel('x');
    ylabel('K(x_i,x)');
    title('Gaussian Kernel - 1D');
    
else
    % Now 2 D
    
    
    sigma = 0.5;
    
    xi = [0.0, 0.0];
    grid = -1:0.02:1;
    
    K = [];
    for i = 1:length(grid)
        for j = 1:length(grid)
            x = [grid(i),grid(j)];
            t1 =  Gaussian_kernel(xi,x,sigma);
            K(i,j) = t1;
        end
        i
    end
    
    
    figure;
    mesh(grid,grid,K);
    xlabel('x_1');
    ylabel('x_2');
    zlabel('K(x_i,x)');
    title('Gaussian Kernel - 2D');
    
    
end


tmp = 99;