function y=gn_residual(y_data, [a,c])

for i=1:size(y_data)
  y(i)=y_data(i)-compeq([a,c],i);
end