function MyReg()
close all
global P

P=[[-1,1];[1,-10];[2,1];[3,5]];

c0 = [0 0];

c = lsqnonlin(@Residuum,c0);

xx=min(P(:,1)):0.1:max(P(:,1));
yy=Ansatzfunktion(xx,c);
plot(xx,yy,'k-',P(:,1),P(:,2),'ro');
grid on

end

function y=Residuum(c)
global P

    y = P(:,2) - Ansatzfunktion(P(:,1),c);

end

function y=Ansatzfunktion(x,c)
    y = c(1)*exp(c(2)*x);
end