This program compares the errors of different time-stepping schemes and
shows how they scale with dt, for linearized Navier-Stokes equations
and Orr-Sommerfeld eigenfunctions as initial data.

The DNS data is compared to the analytic expression for Orr-Sommerfeld 
growth in order to determine the integration error for t=0 to t=T. The
error is computed for a variety of values of dt and time-stepping schemes
to produce data for various plots. 

View results by running matlab script ploterrscaling.m

Terminology for the C++ code and Matlab scripts:
un is the DNS velocity field (u numerical) 
ul is the exact velocity field (u linear)