The use of Richardson extrapolation in PNs solutions of rocket nozzle flow