conv<- function(n){u<-matrix(nrow=n+1,ncol=1)

                   for (i in 0:n) { 

                     if(i==0) f<-function(i) {pnorm(i,m,dp)}

                     else     f<- function(i) {pnorm(i*1/n,m,dp)}

                     if(i==0) p<- function(i) {f(i)}

                     else     p<- function(i) {f(i)-f(i-1)}

                     I<- function (i) {-p(i)*log2(p(i))}

                     u[i,]<- I(i)}

                   y<-u[0:n,]; y}

Mas para o caso onde i==0 não está sendo feito corretamente, e preciso idendificar se não tenho problemas na condição f(i)-f(i-1), o que ele retorna no lugar do i=0 é o valor do i=1.