function fd=fder(beta,y,x,c);
[n,k]=size(x);
fd=-((1-c)*ones(1,k)).*x+(((y-80).*exp(-x*beta))*ones(1,k)).*x;
fd=-sum(fd)';