equalizedImage = double(equalize(:));
    inputImage = double(image(:));
    A = ['mean: ', num2str(mean(inputImage)), ' standard deviation: ', num2str(std(inputImage))];
    B = ['mean: ', num2str(mean(equalizedImage)), ' standard deviation: ', num2str(std(equalizedImage))];
    disp('orignal image:');
    disp(A);
    disp('new image:');
    disp(B);