function [zk,ck,r_inf,zk2,ck2,r2_inf]=cram_coeffs(n) % set the coefficients of the Chebyshev Rational Approximation Method (CRAM) for % exponentials % (c) 2012 Alain Hebert, Ecole Polytechnique de Montreal if n==2 zk=[0.585051560655138-1.185847251723686i, 0.585051560655138+1.185847251723686i]; ck=[0.169152633611546+0.809801111544530i, 0.169152633611546-0.809801111544530i]; r_inf=0.007390314885965; zk2=[1.555307840781615-1.514629903693033i, 1.555307840781615+1.514629903693033i]; ck2=[-0.429247080950924+1.110125209704511i, -0.429247080950924-1.110125209704511i]; r2_inf=0.001591934015977; elseif n == 4 zk=[-0.367838314399776-3.658133272063430i, 1.548400570539483-1.191825853927648i, ... 1.548400570539483+1.191825853927648i,-0.367838314399776+3.658133272063430i]; ck=[0.073392419234124-0.450004915853843i, -0.061683522554915+1.905059455980079i, ... -0.061683522554915-1.905059455980078i, 0.073392419234124+0.450004915853843i]; r_inf=8.647853972498787e-05+1.110223024625157e-16i; zk2=[0.648678605115787-4.229027021806990i, 2.497353528639151-1.379759162356832i, ... 2.497353528639151+1.379759162356832i, 0.648678605115787+4.229027021806990i]; ck2=[0.273989489697394-0.140735983610318i, -0.771863624189269+1.801958418903218i, ... -0.771863624189269-1.801958418903219i, 0.273989489697394+0.140735983610318i]; r2_inf=1.045285816358987e-05-8.326672684688674e-17i; elseif n == 6 zk=[-1.781988275920806-6.196512467345827i, 1.158552571719549-3.614772600819200i, ... 2.400602938933259-1.193129308402001i, 2.400602938933259+1.193129308402001i, ... 1.158552571719549+3.614772600819200i, -1.781988275920806+6.196512467345827i]; ck=[-0.083581617156413+0.106429260747821i, 0.663006870528762-1.451412919902496i, ... -0.579013004031147+4.286888564581573i, -0.579013004031152-4.286888564581576i, ... 0.663006870528764+1.451412919902496i, -0.083581617156413-0.106429260747821i]; r_inf=1.008500923549427e-06-2.844946500601964e-16i; zk2=[-0.737510927460956-6.867098018080055i, 2.136001046696699-4.009943800099992i, ... 3.342791232771612-1.324928206432743i, 3.342791232771612+1.324928206432743i, ... 2.136001046696699+4.009943800099992i, -0.737510927460956+6.867098018080055i]; ck2=[-0.056064744664036-0.009011749235991i, 0.911509924234118-0.379653018817892i, ... -1.355398796202293+3.121023452506399i, -1.355398796202295-3.121023452506398i, ... 0.911509924234118+0.379653018817893i, -0.056064744664036+0.009011749235991i]; r2_inf=8.452975430728316e-08+2.498001805406602e-16i; else error('invalid value of n') end