def gaussoptm(l,x): m= numpy.array([l[0],l[1]]) V= numpy.array([[l[2],l[4]],[l[4],l[3]]]) return -numpy.sum(logmultiskewnormal(x,m,V,[0.,0.]))
def optm(l,x): m= numpy.array([l[0],l[1]]) V= numpy.array([[l[2],l[4]],[l[4],l[3]]]) a= numpy.array([l[5],l[6]]) return -numpy.sum(logmultiskewnormal(x,m,V,a))