As mentioned in discussion, the only variables you need to consider in tsd.data are the following:

omega (this is the relaxation factor=1.8 at low speed and need to be reduced to 1.0 or less for tough cases)
avis=0.01 (remains at this value or close to it)
dm=0.0 (remains at zero)
em=0.1 (remains at 0.1)
ithick=3 (remains at 3 for Quasi-Joukowski thickness distribution)
alphad=10 (incidence in degrees remains at 10 deg.)
amu=2.4 (remains at (gamma+1)=2.4 for air)
amach0= (will vary with each question)

Edit and save tsd.data
Run tsd.exe. Iterate for some time (several 10,000 iterations in subsonic, much less in supersonic)
Save all the Cl, Cm,o, Cd values outputted on the screen.

Each plot will have a curve, which you will plot yourself using Excel and
the equations given in class. These are the linear results.

Run cpcont.m in Matlab, after setting the working directory to the
directory containing tsd.exe output files.

Good luck and be patient with the code!