
%
clear all
close all

 disp('Quetlet Data Set')
 lw = 2; 
 set(0, 'DefaultAxesFontSize', 16);
 fs = 15;
 msize = 10;
%

quetelet0=[
3,33;...
18,34;...
81,35;...
185,36;...
420,37;...
749,38;...
1073,39;...
1079,40;...
934,41;...
658,42;...
370,43;...
92,44;...
50,45;...
21,46;...
4,47;...
1,48];
figure(1)
bar(quetelet0(:,2), quetelet0(:,1))
colormap white

% 
quetelet1=[33 .* ones(1,3) ...
34 .* ones(1,18) ...
35 .* ones(1,81) ...
36 .* ones(1,185) ...
37 .* ones(1,420) ...
38 .* ones(1,749) ...
39 .* ones(1,1073) ...
40 .* ones(1,1079) ...
41 .* ones(1,934) ...
42 .* ones(1,658) ...
43 .* ones(1,370) ...
44 .* ones(1,92) ...
45 .* ones(1,50) ...
46 .* ones(1,21) ...
47 .* ones(1,4) ...
48 .* ones(1,1)];
figure(2)
histn(quetelet1,32,1,49)
colormap cool
mi=mean(quetelet1)
st=std(quetelet1)
hold on
x=32:0.01:49;
y = normpdf(x, mi, st);
plot(x, y, 'b-','linewidth',3)
% print -deps 'C:\Springer\Norm\Normeps\normquetelet.eps'




