function [ reconstructed_image, quantized_coefficients ] = mSimpleDCT8x8ImageAndCoeffs( image, step_size )

image_mean = round( mean(mean(image)) );

M = 8;

[ coefficients, image_size_vector ] = mDCT8x8( image - image_mean );

quantized_coefficients = quant( coefficients, step_size );

reconstruction_coefficients = quantized_coefficients;

DC_coefficients = quantized_coefficients( 1:M:image_size_vector( 1 ), 1:M:image_size_vector( 2 ) );
DPCM_DC_coefficients = mDPCMEncode( DC_coefficients );
quantized_coefficients( 1:M:image_size_vector( 1 ), 1:M:image_size_vector( 2 ) ) = DPCM_DC_coefficients;

reconstructed_image = round( mInverseDCT8x8( reconstruction_coefficients, image_size_vector ) ) + image_mean;