gaz(x,d)

function [y, dyda] = gaz(x,d)

d = d';
y0 =(1+x(1)*d(:,1));
y= y0.^x(2);
dyda(:,1) = x(2).*y0.^(x(2)-1).*d(:,1);
dyda(:,2) = y.*log(y0);
y = y-d(:,2);