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.