function res = my_func(x)
res = x.^2 - 2*x + 1;
end % function my_func