function y=quartic(ggg)

y=((ggg-3).^2).*((ggg-1).^2);