#algoritmo 1.1
f1<-function(sampl){u<-0;v<-0;for(i in 1:length(sampl)){u<-u+sampl[i];v<-v+sampl[i]^2};s2<-(v-u^2/length(sampl))/(length(sampl)-1)}
#
#instavel:
#exemplo:
#x<-rbeta(10000,1,1);
#var(x); #= 0.0831778417453782
#a<-f1(x); #=0.08317784174537904
#
#algoritmo 1.2
f2<-function(sampl){u<-0;for(i in 1:length(sampl)){u<-u+sampl[i]};uu<-u/length(sampl); v<-0; for(i in 1: length(sampl)){v<-v+(sampl[i]-uu)^2};s2<-v/(length(sampl)-1)}
#
#idem exemplo anterior
#b<-f2(x); #=0.0831778417453782
#
#algoritmo 1.3
f3<-function(sampl){u<-0;for(i in 1:length(sampl)){u<-u+sampl[i]};uu<-u/length(sampl); v<-0;w<-0;for(i in 1:length(sampl)){v<-v+(sampl[i]-uu)^2; w<-w+(sampl[i]-uu)};s2<-(v-w^2/length(sampl))/(length(sampl)-1) }
#