function multiscale_eq(input, example)

  addpath('../../../timeLapseCode/TimeLapseVideo/helpers/image_pyramid');
  nLevel  = 5; 

  input_pyramid = LaplacianPyramids(input,nLevel); 
  example_pyramid = LaplacianPyramids(input,nLevel); 

  output_pyramid = {};

  for i = 1 : nLevel
    output_pyramid{i} = HistTransfer(input_pyramid{i}, example_pyramid{i});
  end

  output = reconstruct_laplacian_pyramid{output_pyramid}

  
function output = HistTransfer(input, example)
  output = imhistmatch(input+0.5, example+0.5)-0.5;

