public static double GammaPDF(double t, GammaParams gammaParams)
{
//object gampdf_result = null;
//Matlab.Feval("gampdf", 1, out gampdf_result, t, gammaParams.a, gammaParams.b);
//object[] res = gampdf_result as object[];
//double result = ((double)res[0]);
//return result;
return(Gamma.PDF(gammaParams.a, 1 / gammaParams.b, t));
//MWNumericArray m_a = new MWNumericArray(gammaParams.a);
//MWNumericArray m_b = new MWNumericArray(gammaParams.b);
//MWNumericArray m_t = new MWNumericArray(t);
//MWArray res = matlabFunctions.GammaPDF(m_t, m_a, m_b);
//double result = ((double)(res.ToArray().GetValue(0,0)));
//return result;
}
public GammaParams Gamma; // если используем гамма-распределение
public PiParams(double _norm, GammaParams _gamma)
{
Norm = _norm;
Gamma = _gamma;
}
public static double GammaCDF(double t, GammaParams gammaParams)
{
return(Gamma.CDF(gammaParams.a, 1 / gammaParams.b, t));
}
