function [result]=rv_HW(x1,x2,x3)
n=x1+x2+x3;
parte1=fatorial(n+2)*fatorial(2*x1+x2)*fatorial(2*x3+x2)*2^(x2-1)/(fatorial(x1)*fatorial(x2)*fatorial(x3)*fatorial(2*n+1)*3);
parte2=5-12*(2*x1+x2+1)*(2*x3+x2+1)/((2*n+2)*(2*n+3));
result=parte1*parte2;
