function loglik=loglik(lamda,dur,censor)

loglik=(lamda)*sum(1-censor)-exp(lamda)*sum(dur);