Example is working
Change in plot() y-limits
Probit likelihood modified for plotting.