I have determined the characteristic equations, but am having a hard time implementing them into cpp.

the problem gives dy/dx(+/-) = f(x,y,u,v), and du/dy(+/-) = f(x,y,u,v) through a nozzle defined by y=1(+/-).04x^2 as boundaries. i have used method of characteristics to determine the characteristic equations, but i cant figure out how to write the code. If there is anyone who has a similar code, or an example that would help me, it would be apreciated. thanks

Found on the web (I think I wrote my 2D characteristic book based on the book of Zucrow): Gas Dynamics: by E. Rathakrishnan Prentice-Hall India Pvt Ltd, Pages: 408, Price: 195 2001. ISBN: 81-203-0952-9.

"Gas Dynamics" by Prof. E. Rathakrishnan is a simple and direct introduction to the study of compressible gas flows. The book is easily read and understood by an undergraduate student, as well as by a practicing engineer. It is written in a lucid style, with numerous solved numerical examples to illustrate physical principles. The author offers many practical insights into experimental gas dynamics at different places in the book, especially in Chapter 13, which is devoted to measurements of compressible flows. Normal shock tables, oblique shock tables (both weak and strong shocks), tables for one-dimensional viscous flows, and a FORTRAN program for designing axisymmetric supersonic nozzles by the method of characteristics (all included in an appendix) are both useful and instructive. The book is available at an affordable price and thus has the potential for becoming a popular textbook on gas dynamics. A logical arrangement of the essential topics, numerous figures, and the simple approach adopted in the book also make it an attractive and refreshing guide to a worker in the field of compressible fluid dynamics.

The contents of the book are similar to those found in the gas dynamics classics, such as the books by Shapiro (1953), Liepmann and Roshko (1957), and Zucrow (1976). It is the down-to-earth basic treatment which sets this book apart from the classics.