%% prog003 promedio desviacion estandard

clearvars;    %borra variables en la memoria / workspace
close all;    %cierra ventanas de graficos
clc;          %borra todos los mensajes de conmmand window

%% entrada de datos
disp('programa que calcula promedio y desviacion estandard');
disp(' ');
disp('             N            ');
disp('           -----          ');
disp('           \              ');
disp('prom =  1   \     x(n)    ');    
disp('       ---  /             ');
disp('        N  /              ');
disp('           ------         ');
disp('            n = 1         ');
disp(' ');
disp(' ');
disp('         |        N                   | ^(1/2)');
disp('         |     ------                 |');
disp('         |     \                      |');
disp('dsvstd = |  1   \    ( x(n)-prom)^2   |');    
disp('         | ---  /                     |');
disp('         |  N  /                      |');
disp('         |     ------                 |');
disp('         |      n = 1                 |');

x1 = input('x1 = ');
x2 = input('x2 = ');
x3 = input('x3 = ');
x4 = input('x4 = ');
x5 = input('x5 = ');

%% calculos
promedio = (x1 + x2 + x3 + x4 + x5)/5;

%una forma es escribir el comando completo en una sola linea
%lo cual es incomodo y fuente de posibles errores

%dsvstd = sqrt( ((x1-promedio)^2 + (x2-promedio)^2 + (x3-promedio)^2 + (x4-promedio)^2 + (x5-promedio)^2)/5 ) 

%una linea de comando puede subdividirse en partes y puede ser concatenada 
%con el operador tres puntos (...)
%usaremos ademas una variables auxiliares suma y suma5

suma = (x1 - promedio)^2 + ...
       (x2 - promedio)^2 + ...
       (x3 - promedio)^2 + ...
       (x4 - promedio)^2 + ...
       (x5 - promedio)^2;
       
suma5 = suma/5;

dsvstd = sqrt(suma5);


%% escribir resultados
disp(' ');
disp('resultados tipo 1 usando disp()');
disp('promedio = ');
disp(promedio);
disp('desviacion estandard = ');
disp(dsvstd);

fprintf(' \n'); 
fprintf('resultados tipo 2 usando fprintf() \n');
fprintf('promedio             = %f\n',promedio);
fprintf('desviacion estandard = %f\n',dsvstd);