disc2class(normden1, normden2, prior1, prior2, range, n)

function disc2class(normden1, normden2, prior1, prior2, range, n)
%function disc2class(normden1, normden2, prior1, prior2, range, n)
% Create Decision Boundary between two classes from given density

x = linspace(range(1), range(2), n);
[row col] = size(normden1);
if (isequal(size(normden1),size(normden2))==1)&...
        (row == col)&(row~=1)
    gi = (log(normden1)+log(prior1))-(log(normden2)+log(prior2));
    [xax, yax] = find(partition(gi)~=0);
    plot(x(yax), x(xax), 'k.','MarkerSize',5);
    axis([range(1) range(2) range(1) range(2)])
elseif (isequal(size(normden1),size(normden2))==1)&(row==1)
    gi = (log(normden1)+log(prior1))-(log(normden2)+log(prior2));
    xax = find(partition(gi)~=0);
    for l = 1:length(xax)
        line([x(xax(l)) x(xax(l))], [0 max(max(normden1),max(normden2))], 'Color','g')
    end
else    
    errordlg('Matrix Sizes are not correct');
end