function y = biwt (x) 
	y = ((1-x.^2).^2).*(abs(x)<1);



