function ll_fder=loglik_fder(y,c,beta)
%ll=-beta*sum(1-c)-sum(y*exp(-beta));
ll_fder=-sum(1-c)+sum(y*exp(-beta));