#!/usr/bin/octave -qf

if (size(argv) != 2)
  printf("Usage: prob.m <a> <b>\n");
  exit(1);
endif

function y = prob(x,a,b)
  y = 1/beta(a,b) * x**(a-1) * (1-x)**(b-1);
endfunction

x = [0:.01:1];
a = str2num(argv{1});
b = str2num(argv{2});

fd = fopen(sprintf("prob-a%g-b%g.dat", a, b), "wt");
for j = 1:size(x,2)
  y = prob(x(j), a, b);
  if (y != inf)
    fprintf(fd, "%g \t %g\n", x(j), y);
  endif
endfor