clc;
close all;
clear all;

rath = load('rath.txt');
leuk = load('leuk.txt');

for k = 1:7
    
    figure;
    subplot(2, 1, 1);
    hist(rath(:, k), 50);
    subplot(2, 1, 2);
    hist(leuk(:, k), 50);
%     [N, Xr] = hist(rath(:, k)); Nr = N/trapz(Xr, N);            
%     [N, Xl] = hist(leuk(:, k)); Nl = N/trapz(Xl, N);
%     plot(Xr, Nr, 'b', Xl, Nl, 'r' );        
end        