function plot_beta(a,b)
  x = linspace(0, 1, 1000)';
  
  plot(x, betapdf(x, a, b));
  xlim([0 1]);
end