function zz=Leg4(x)
zz=(1/sqrt(2/9))*1/8*(35*x.*x.*x.*x-30*x.*x+3*ones(size(x)));