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