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