正常反向伽玛问题

Question

我尝试使用matalb实现Normal Inverse Gamma功能，但是我得到了错误。

function [ INGN ] = inversegamNormpdf(X,A,B,L,U,S )

INGN=(sqrt(L)/S*sqrt(2*pi))*(B^A/gamma(A))*((1/S^2)^A+1)*exp(-2*B+L*(X-U)^2/2*S^2);

end 

然后将函数调用为 其中x是由260 * 4组成的数据。

INGN  = inversegamNormpdf(x,A,B,L,U,S );

我不知道如何解决它？