 = @(X) gaussProb(X, Xbar, sigma.^2);
    model = studentFit(data);
    sT = @(X)exp(studentLogprob(model, X));
    model = laplaceFit(data);
    lap = @(X)exp(laplaceLogprob(model, X));
    hold <span class="string">on</span>;
    x = (-5:0.01:10)';
    h(1) = plot(x,gauss(x),<span class="string">'k:'</span>,<span class="string">'LineWidth'</span>,3);
    h(2) = plot(x,sT(x),<span class="string">'r-'</span>,<span class="string">'LineWidth'</span>,3);
    h(3) = plot(x,lap(x),<span class="string">'b--'</spa