model = studentFit(data);
    sT = @(X)exp(studentLogprob(model, X));
    model = laplaceFit(data);
    lap = @(X)exp(laplaceLogprob(model, X));   
    hold on;
    x = (-5:0.01:10)';
    h(1) = plot(x,gauss(x),'k:','LineWidth',3);
    h(2) = plot(x,sT(x),'r-','LineWidth',3);
    h(3) = plot(x,lap(x),'bREPLACE_WITH_DASH_DASH','LineWidth',3);
    axis([-5,10,0,0.5]);
    set(gca,'YTick',0:0.1:0.5);
    legendStr = {'gaussian', 'student T', 'laplace'};
    if isOctave(),
        legend(legendStr{:})
    else
