# funcao_dp
# Autor
# 12/03/2020
# Criar uma função para calcular a desvio padrao usando o source 
rm(list = ls())

source('C:/Users/Mie/Documents/ESALQ USP Disciplina/LCE 5783/scritp LCE 5783/script_media.r')
d_padrao <- function(vetor){
  # Tamanho do vetor
  n <- length(vetor)
  s1 <- 0
  for(i in 1:n){
    s1 <- s1+((vetor[i]-media(vetor))^2)
  }
  s2 <- s1/(n-1)
  sd <- sqrt(s2)
  return(sd)
}
leiteA <- c(38.2, 41.2, 41.1, 36.8, 43.6, 34.7, 36.6, 36.6, 36.3, 33.4)
d_padrao(leiteA)