function y=flagr4(x)
y=(x.^2-6*x+5).^2;
