Hello.
2. 20Ne+p
3. a)
  so...
  RP 20 20 10 ;
  DL .1 ;
  DP 1. 45. .05 ;
  DL .1 ;
  
3. b)
   1.000000     0.1665335E-15 0.0000000E+00 0.0000000E+00-0.4144358     100000
  0.9071068      1.000000     0.0000000E+00 0.0000000E+00-0.1879688     010000
  0.0000000E+00 0.0000000E+00 0.6053100    -0.6936741     0.0000000E+00 001000
  0.0000000E+00 0.0000000E+00 0.9133969     0.6053100     0.0000000E+00 000100
  0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00  1.000000     000010
  0.1879688     0.4144358     0.0000000E+00 0.0000000E+00 0.2264914     000001
  
  We found it easier to change the line with 'PM' to 'WM', which gave us the below output.
  
       I  COEFFICIENT            ORDER EXPONENTS
     1   1.000000000000000       1   1 0  0 0  0 0 [magnification (x)]
     2  0.9071067811865474       1   0 1  0 0  0 0 [phase rotation]
     3  0.1879687538266329       1   0 0  0 0  0 1 [dispersion dx/de]
     ---------------------------------------------
     I  COEFFICIENT            ORDER EXPONENTS
     1  0.1665334536937735E-15   1   1 0  0 0  0 0 
     2   1.000000000000000       1   0 1  0 0  0 0 
     3  0.4144357813768277       1   0 0  0 0  0 1
     ---------------------------------------------
     I  COEFFICIENT            ORDER EXPONENTS
     1  0.6053100184945539       1   0 0  1 0  0 0
     2  0.9133969081326894       1   0 0  0 1  0 0
     ---------------------------------------------
     I  COEFFICIENT            ORDER EXPONENTS
     1  -.6936741036330284       1   0 0  1 0  0 0
     2  0.6053100184945541       1   0 0  0 1  0 0
     ---------------------------------------------
     I  COEFFICIENT            ORDER EXPONENTS
     1  -.4144357813768276       1   1 0  0 0  0 0
     2  -.1879687538266329       1   0 1  0 0  0 0
     3   1.000000000000000       1   0 0  0 0  1 0
     4  0.2264914342847914       1   0 0  0 0  0 1
     ---------------------------------------------
     I  COEFFICIENT            ORDER EXPONENTS
     1   1.000000000000000       1   0 0  0 0  0 1
     ---------------------------------------------
     ALL COMPONENTS ZERO 
     ------------------- 
     ALL COMPONENTS ZERO 
     ------------------- 
  
  
3. c)
   We changed the sixth term being passed to the 'SB' command. (PD below) to 0.01.
   SB <PX><PA><r12><PY><PB><r34><PT><PD><r56><PG><PZ>
   
   This pulled the three sets of rays in the image very close together.
   
  
4. a)
   Target                     1H
   Projectile                 20Ne
   Ejectile                   g
   Recoil                     21Na
   
   Projectile Kinetic Energy  40MeV
   
   We plotted theta4(recoil angle) vs theta3(gamma angle)
   It seems the larger the recoil cone, the larger the gamma cones (equally forward & backward)
   Ours did fit both separators. If it didn't, the larger the energy, the smaller the recoil cones.
   
4. b)
   We plotted theta4 (recoil angle) vs E4 (recoil E), and found dE/E = +/- 0.7%
   
5. 
   https://journals.aps.org/prc/abstract/10.1103/PhysRevC.97.065802
   
7. Q2 converged at, Mas. Res.Power at FP2= 751, Resol. by max. ray= 950. This was after the fit, minimizing the inverse of the first order resolving power. Which didnt differ from before changing length.
B before was=0.111727
B after was= 0.115052 -> 3% different
These were done with cosy 9.1 FRIB which is not the optimal version to use. Version 10 gave close numbers (747 and 946). did not have time to redo with other quads and other version. 
8.