function f=testfunction(s,pts)
f=4^s*prod(pts.*(1-pts),2);
return