Solution of compressible flow inside a CD nozzle with shock formation. Implementation with Mac-Cormack scheme