function deltasixj = deltasixj(a,b,c)



deltasixj = sqrt( (factorial(a+b-c)*factorial(a-b+c)*factorial(-a+b+c))/factorial(a+b+c+1) );



