
function y=g1(x)

y=x.^2+4*x-1;