Compare analytic calculation of d(beta)/d(delta) and d(eta)/d(delta) with "brute-force" method (change energy, recalculate 6x6 matrices, and find new Twiss parameters).

Program output: dbeta_ddelta.out