Como calcular desvio padrão em C++ - C++ para Matemática e Estatística

#include <iostream>
#include <math.h>

using namespace std;

// função principal do programa
int main(int argc, char *argv[]){
  // conjunto dos dados
  double conjunto[] = {10, 30, 90, 30};
  double soma = 0.0; // soma dos elementos
  double desvio_padrao = 0.0; // desvio padrão
  int tam = 4; // tamanho dos dados
  double media;

  // vamos somar todos os elementos
  for(int i = 0; i < tam; i++){
    soma = soma + conjunto[i];
  }

  // agora obtemos a média do conjunto de dados    
  media = soma / tam;

  // e finalmente obtemos o desvio padrão
  for(int i = 0; i < tam; i++){