function su = vonmises_likelihood(theta, kappa_max, samples)
t = reshape(theta, [1 1 length(theta)]);
[mu,kappa] = meshgrid(linspace(-pi,pi,samples), linspace(0,kappa_max,samples));
S = prod(vonmises(t, mu, kappa), 3);
su = surf(mu,kappa,S, 'EdgeColor', 'none');