eCulpritIntegrals.m

function [ I,dIdX, Idldg1, Idldg2, I1xx, I2xx,I1lgx, I2lgx ] = eCulpritIntegrals( delta,ddeltada,d2deltada2,e,~,l,n,m,p , Jacobian2On)
%eCulpritIntegrals Culprit Integrals and its partial derivatives using
%eccentricity expansions generated using Maple symbolic software.
% Author: Bharat Mahajan (https://github.com/princemahajan)


t1 = (n - p);
t2 = (p * t1);
t3 = (t2 * n);
t4 = (p ^ 2);
t5 = (t4 ^ 2);
t6 = (p * t4);
t7 = 27 * n;
t8 = -432 * p;
t9 = (((((t7 + t8 + 261) * n + (2736 * p - 2688) * p + 274) * n + ((-8704 * p + 9600) * p + 524) * p - 1476) * n + t4 * ((14592 * p - 13824) * p - 9804) - 90) * n + t6 * ((-12288 * p + 6912) * p + 18560) + 1764) * n + t5 * (4096 * t4 - 9280) - 6136 * t2 + 8844 * t3 - 720;
t10 = (delta * m);
t11 = n - t10;
t12 = 2;
t13 = t12 * p;
t14 = t13 - t11;
t15 = -l * t11;
t16 = t13 * l;
t17 = t16 + t15;
t18 = n ^ 2;
t19 = t18 ^ 2;
t20 = n * t18;
t1 = t1 * t6;
t21 = 4096 * p;
t22 = 468952 * p;
t23 = t13 - 6 - t11;
t24 = (t13 - 6) * l + t15;
t25 = 243 * n;
t26 = 59192 * p;
t27 = t13 - 4 - t11;
t28 = (t13 - 4) * l + t15;
t29 = sin(t28);
t30 = 81 * n;
t31 = -3496 * p;
t32 = 1 - p;
t33 = -t12 * t32 - t11;
t34 = -t12 * l * t32 + t15;
t35 = sin(t34);
t36 = 9 * n;
t37 = 1 + p;
t38 = t12 * t37 - t11;
t39 = t12 * l * t37 + t15;
t40 = sin(t39);
t41 = 1792;
t42 = 3 * n;
t43 = -64 * p;
t44 = t13 + 4 - t11;
t45 = (t13 + 4) * l + t15;
t46 = sin(t45);
t47 = 1 - n;
t48 = 1280;
t49 = 24;
t50 = t13 + 6 - t11;
t51 = (t13 + 6) * l + t15;
t50 = 1 / t50;
t23 = 1 / t23;
t23 = (-34560 * t1 * t18 + p * ((((-73710 * n + 194400 * p) * n - 256320 * t4 - 760620) * n + (168960 * t4 + 975480) * p - 745376) * n + t22) + (((((729 * n - 5832 * p + 11178) * n + 19440 * t4 + 66960) * n - 348840 * p + 197580) * n + 681300 * t4 + 295864) * n + t6 * (-18432 * t4 - 591200) + 204576) * n + (t4 * (((t21 - 44544) * p + 192320) * p - 416760) - 251232) * p + 46656) * t23;
t22 = (p * ((((750 * n - 6240 * p) * n + 25920 * t4 + 59940) * n + (-53760 * t4 - 274800) * p - 192528) * n + t22) + 46656 * t47 + ((((-t18 * t48 - 178080) * n + ((-6144 * n + t21 + 44544) * p + 3840 * t18 + 192320) * p + 416760) * p + t18 * (240 * t18 + 61620)) * p + t20 * (-t18 * t49 - 9440) + 251232) * p + t18 * ((((-36 + n) * n + 540) * n - 4320) * n + 19440)) * t50;
t27 = 1 / t27;
t44 = 1 / t44;
t50 = (p * ((((940 * n - 7040 * p - 5340) * n + 25600 * t4 + 14900) * n + (-44800 * t4 - 69900) * p - 21644) * n + t26) + (((((t43 + t42 - 49) * n + 560 * t4 + 310) * n - 2560 * t6 - 960) * n + t4 * (6400 * t4 + 32980) + 1600) * n + t6 * (-8192 * t4 - 86400) - t41) * n + (t4 * (((t21 + 29696) * p + 80320) * p + 100800) + 14808) * p + 1536) * t44;
t26 = (p * ((((-32940 * n + 97920 * p - 122700) * n - 143360 * t4 - 177500) * n + (103680 * t4 + 232500) * p - 96740) * n + t26) + (((((-2592 * p + 4347 + t25) * n + 10800 * t4 + 21870) * n - 23040 * t6 + 44840) * n + t4 * (26880 * t4 + 255700) + 39148) * n + t6 * (-16384 * t4 - 234880) + 13016) * n + (t4 * (((t21 - 29696) * p + 80320) * p - 100800) - 14808) * p + 1536) * t27;
t33 = 1 / t33;
t38 = 1 / t38;
t32 = (-784 * n * t32 + p * ((((378 * n - 3616 * p + 1032) * n + 14272 * t4 - 2924) * n + p * (-24576 * t4 + 7336)) * n + t31) + (-168 * n * t19 + (((-4864 * t18 - 8160) * n + ((-10240 * n + t21 + 14848) * p + 9984 * t18 + 13120) * p - 2712) * p + t18 * (1264 * t18 - 684)) * p + 1248) * p + t18 * ((((t36 - 10) * n - 136) * n + 276) * n + 224) + 384) * t38;
t21 = (p * ((((-11898 * n + 40224 * p - 27664) * n - 64448 * t4 - 3612) * n + (49664 * t4 - 800) * p + 6208) * n + t31) + (((((-1080 * p + 1296 + t30) * n + 5616 * t4 + 5172) * n - 14592 * t6 + 1976) * n + t4 * (20224 * t4 + 53556) - 2488) * n + t6 * (-14336 * t4 - 44320) + 464) * n + (t4 * (((t21 - 14848) * p + 13120) * p + 2712) - 1248) * p + 384) * t33;
t14 = 1 / t14;
t31 = -(0.1e1 / 0.7680e4);
t52 = 0.1e1 / 0.3072e4;
t53 = 0.1e1 / 0.46080e5;
t54 = t14 * sin(t17);
t55 = -1024 * p;
t56 = t13 - 5 - t11;
t57 = (t13 - 5) * l + t15;
t58 = t13 - 3 - t11;
t59 = (t13 - 3) * l + t15;
t60 = sin(t59);
t61 = 1 + n - t10;
t62 = t13 - t61;
t61 = -t61 * l + t16;
t63 = sin(t61);
t64 = 80;
t10 = 1 - n + t10;
t65 = t13 + t10;
t10 = t10 * l + t16;
t16 = sin(t10);
t66 = t13 + 3 - t11;
t67 = (t13 + 3) * l + t15;
t68 = sin(t67);
t69 = -640 * n;
t11 = t13 + 5 - t11;
t13 = (t13 + 5) * l + t15;
t15 = 1 / t58;
t58 = 1 / t66;
t66 = (((((t30 + 999) * n + 2814) * n + 2598) * n + 649) * n + ((((-756 * n - 6012) * n - 11500) * n - 6432) * n + (((2592 * n + 13296) * n + 15500) * n + ((-4224 * n - 12864) * n + (3328 * n + t55 + 4608) * p - 6896) * p + 3924) * p - 676) * p + 27) * t15;
t41 = (((((-52 * n + 420) * n - 1188) * n + 1416) * n + (((352 * n - 2352) * n + 5188) * n + ((-1152 * n + 5568) * n + (n * t41 + t55 - 4608) * p - 6896) * p - 3924) * p - 676) * p + t18 * ((3 * t18 + 82) * n - 90) - 27 * (-1 + t20) * n - 27) * t58;
t56 = 1 / t56;
t11 = 1 / t11;
t11 = (((((-20 * n + 420) * n - 3340) * n + 11960) * n + (((160 * n - 2640) * n + 14780) * n + ((t69 + 7360) * n + (n * t48 + t55 - 7680) * p - 21680) * p - 28220) * p - 16348) * p + t18 * (((-25 + n) * n + 250) * n - 1250) - 3125 * t47) * t11;
t25 = (((((t25 + 2565) * n + 9990) * n + 17510) * n + 13223) * n + ((((-1620 * n - 13500) * n - 38820) * n - 44480) * n + (((4320 * n + 26640) * n + 50260) * n + ((-5760 * n - 23360) * n + (3840 * n + 7680 + t55) * p - 21680) * p + 28220) * p - 16348) * p + 3125) * t56;
t48 = 1 / t62;
t56 = 1 / t65;
t62 = (((((t36 + 23) * n - t64) * n - 30) * n + 139) * n + ((((-132 * n - 108) * n + 708) * n - 284) * n + (((736 * n - 240) * n - 1548) * n + ((-1920 * n + 1600) * n + (2304 * n + t55 - 1536) * p + 496) * p + 860) * p - 184) * p - 61) * t56;
t55 = (((((t7 + 261) * n + 424) * n - 546) * n + 45) * n + ((((-324 * n - 1932) * n - 900) * n + 1436) * n + (((1440 * n + 4656) * n + 60) * n + ((-2944 * n - 4544) * n + (2816 * n + t55 + 1536) * p + 496) * p - 860) * p - 184) * p + 61) * t48;
t65 = (0.1e1 / 0.768e3);
t70 = 0.1e1 / 0.3840e4;
t71 = 0.1e1 / 0.384e3;
t3 = (((-96 * p + t36 + 38) * n + 352 * t4 - 14) * n - 512 * t6 - 50) * n + 196 * t2 - 192 * t3 + t49 + 256 * t5;
t6 = 1216 * p;
t72 = 16;
t73 = 256 * p;
t1 = (-t72 * t20 * t37 + ((204 * n - 880) * n + (-864 * n + 1996 + t6) * p + 1292) * p - 256 * n - 256 * t1 + t19 + 96 * t18 * (1 + t4) + 256) * t44;
t4 = ((((t30 + t8 + 540) * n + (864 * p - 2124) * p + 1212) * n + ((-768 * p + 2784) * p - 3112) * p + 1036) * n + ((-t6 + 1996) * p - 1292) * p + 256 * t5 + 256) * t27;
t5 = ((((t42 - 11) * n + 6) * n + 12) * n + (((-40 * n + 138) * n - 126) * n + ((192 * n - 528) * n + (-384 * n + 608 + t73) * p + 352) * p + t72) * p - 8) * t38;
t6 = ((((t7 + 207) * n + 232) * n + 28) * n + (((-216 * n - 906) * n - 578) * n + ((576 * n + 1296) * n + (t73 + t69 - 608) * p + 352) * p - t72) * p - 8) * t33;
t7 = -(0.1e1 / 0.96e2);
t8 = 1 + t20;
t15 = (n * (99 * n + 103) + ((-108 * n - 258) * n + (144 * n + t43 + 168) * p - 130) * p + 27 * t8) * t15;
t19 = (((-12 * n + 78) * n + (48 * n + t43 - 168) * p - 130) * p + t18 * (n - 9) - 27 * t47) * t58;
t8 = (-n * (n + 5) + ((-28 * n + 22) * n + (n * t64 + t43 - 56) * p + 6) * p + 3 * t8) * t56;
t20 = (((t36 + 35) * n - 1) * n + ((-60 * n - 90) * n + (112 * n + t43 + 56) * p + 6) * p - 3) * t48;
t27 = -(0.1e1 / 0.16e2);
t30 = (0.1e1 / 0.48e2);
t2 = -t2 * t72 - t12 + 3 * n * (n + 1);
t12 = t72 * p;
t18 = ((-8 * n + 18 + t12) * p + t18 + 4 * t47) * t38;
t12 = (n * (t36 + 14) + (-n * t49 + t12 - 18) * p + 4) * t33;
t33 = (0.1e1 / 0.8e1);
t36 = 4 * p;
t37 = (-t36 - t47) * t56;
t36 = (1 + t42 - t36) * t48;
t38 = e ^ 2;
t42 = (0.1e1 / 0.2e1);
t28 = cos(t28);
t34 = cos(t34);
t39 = cos(t39);
t43 = cos(t45);
t14 = (t14 * cos(t17));
t17 = cos(t59);
t44 = cos(t61);
t10 = cos(t10);
t45 = cos(t67);
t10 = ((t42 * (-t37 * t10 + t36 * t44) * e) + t38 * ((((((t31 * (t26 * t28 + t50 * t43)) + t52 * (t21 * t34 + t32 * t39) + t53 * (t23 * cos(t24) + t22 * cos(t51)) - (t14 * t9) / 0.2304e4) * e - (t65 * (t66 * t17 - t41 * t45)) + t70 * (t25 * cos(t57) - t11 * cos(t13)) + t71 * (-t62 * t10 + t55 * t44)) * e + (t7 * (t6 * t34 + t5 * t39)) + t71 * (t1 * t43 + t4 * t28) + (t14 * t3) / 0.64e2) * e + (t27 * (-t8 * t10 + t20 * t44)) + (t30 * (t15 * t17 - t19 * t45))) * e + (t33 * (t12 * t34 + t18 * t39)) - (t14 * t2) / 0.4e1) + t14);
I12 = [t42 * (-(t37 * t16) + t36 * t63) * e + t38 * (((((t53 * (t22 * sin(t51) + t23 * sin(t24)) + t31 * (t26 * t29 + t50 * t46) + t52 * (t21 * t35 + t32 * t40) - t54 * t9 / 0.2304e4) * e + t65 * (t41 * t68 - t66 * t60) - t70 * (t11 * sin(t13) - t25 * sin(t57)) + t71 * (-(t62 * t16) + t55 * t63)) * e + t71 * (t1 * t46 + t4 * t29) + t7 * (t6 * t35 + t5 * t40) + t54 * t3 / 0.64e2) * e + t30 * (t15 * t60 - t19 * t68) + t27 * (-(t8 * t16) + t20 * t63)) * e + t33 * (t12 * t35 + t18 * t40) - t54 * t2 / 0.4e1) + t54 t10];

t1 = (n - p);
t2 = (p * t1);
t3 = (t2 * n);
t4 = (p ^ 2);
t5 = (t4 ^ 2);
t6 = (p * t4);
t7 = 27 * n;
t8 = -432 * p;
t9 = (delta * m);
t10 = n - t9;
t11 = 2;
t12 = t11 * p;
t13 = -t10 + t12;
t14 = -l * t10;
t15 = t12 * l;
t16 = t14 + t15;
t17 = sin(t16);
t16 = cos(t16);
t18 = n ^ 2;
t19 = t18 ^ 2;
t20 = n * t18;
t1 = (t1 * t6);
t21 = 4096 * p;
t22 = 468952 * p;
t23 = -34560 * t1 * t18 + p * ((((-73710 * n + 194400 * p) * n - 256320 * t4 - 760620) * n + (168960 * t4 + 975480) * p - 745376) * n + t22) + (((((729 * n - 5832 * p + 11178) * n + 19440 * t4 + 66960) * n - 348840 * p + 197580) * n + 681300 * t4 + 295864) * n + t6 * (-18432 * t4 - 591200) + 204576) * n + (t4 * (((-44544 + t21) * p + 192320) * p - 416760) - 251232) * p + 46656;
t24 = -t10 - 6 + t12;
t25 = (-6 + t12) * l + t14;
t26 = sin(t25);
t25 = cos(t25);
t27 = 243 * n;
t28 = 59192 * p;
t29 = p * ((((-32940 * n + 97920 * p - 122700) * n - 143360 * t4 - 177500) * n + (103680 * t4 + 232500) * p - 96740) * n + t28) + (((((-2592 * p + 4347 + t27) * n + 10800 * t4 + 21870) * n - 23040 * t6 + 44840) * n + t4 * (26880 * t4 + 255700) + 39148) * n + t6 * (-16384 * t4 - 234880) + 13016) * n + (t4 * (((-29696 + t21) * p + 80320) * p - 100800) - 14808) * p + 1536;
t30 = -t10 - 4 + t12;
t31 = (-4 + t12) * l + t14;
t32 = sin(t31);
t31 = cos(t31);
t33 = 81 * n;
t34 = -3496 * p;
t35 = p * ((((-11898 * n + 40224 * p - 27664) * n - 64448 * t4 - 3612) * n + (49664 * t4 - 800) * p + 6208) * n + t34) + (((((-1080 * p + 1296 + t33) * n + 5616 * t4 + 5172) * n - 14592 * t6 + 1976) * n + t4 * (20224 * t4 + 53556) - 2488) * n + t6 * (-14336 * t4 - 44320) + 464) * n + (t4 * (((-14848 + t21) * p + 13120) * p + 2712) - 1248) * p + 384;
t36 = 1 - p;
t37 = -t11 * t36 - t10;
t38 = -t11 * l * t36 + t14;
t39 = sin(t38);
t38 = cos(t38);
t40 = 9 * n;
t41 = 1 + p;
t42 = t11 * t41 - t10;
t43 = t11 * l * t41 + t14;
t44 = sin(t43);
t43 = cos(t43);
t45 = 1792;
t46 = 3 * n;
t47 = -64 * p;
t48 = -t10 + 4 + t12;
t49 = (4 + t12) * l + t14;
t50 = sin(t49);
t49 = cos(t49);
t51 = 1 - n;
t52 = 1280;
t53 = 24;
t54 = -t10 + 6 + t12;
t55 = (6 + t12) * l + t14;
t56 = sin(t55);
t55 = cos(t55);
t13 = 1 / t13;
t57 = (l * t16) - t13 * t17;
t58 = ((((((261 + t8 + t7) * n + (2736 * p - 2688) * p + 274) * n + ((-8704 * p + 9600) * p + 524) * p - 1476) * n + t4 * ((14592 * p - 13824) * p - 9804) - 90) * n + t6 * ((-12288 * p + 6912) * p + 18560) + 1764) * n + t5 * (4096 * t4 - 9280) - 6136 * t2 + 8844 * t3 - 720) * t13 * m;
t54 = 1 / t54;
t24 = 1 / t24;
t59 = t24 ^ 2;
t22 = (p * ((((750 * n - 6240 * p) * n + 25920 * t4 + 59940) * n + (-53760 * t4 - 274800) * p - 192528) * n + t22) + 46656 * t51 + ((((-t18 * t52 - 178080) * n + ((-6144 * n + t21 + 44544) * p + 3840 * t18 + 192320) * p + 416760) * p + t18 * (240 * t18 + 61620)) * p + t20 * (-t18 * t53 - 9440) + 251232) * p + t18 * ((((-36 + n) * n + 540) * n - 4320) * n + 19440)) * t54;
t48 = 1 / t48;
t30 = 1 / t30;
t60 = t30 ^ 2;
t61 = (l * t31 * t30) - t32 * t60;
t62 = (l * t49) - t48 * t50;
t28 = (p * ((((940 * n - 7040 * p - 5340) * n + 25600 * t4 + 14900) * n + (-44800 * t4 - 69900) * p - 21644) * n + t28) + (((((-49 + t47 + t46) * n + 560 * t4 + 310) * n - 2560 * t6 - 960) * n + t4 * (6400 * t4 + 32980) + 1600) * n + t6 * (-8192 * t4 - 86400) - t45) * n + (t4 * (((29696 + t21) * p + 80320) * p + 100800) + 14808) * p + 1536) * t48;
t42 = 1 / t42;
t37 = 1 / t37;
t63 = t37 ^ 2;
t64 = (l * t38 * t37) - t39 * t63;
t65 = (l * t43) - t42 * t44;
t21 = (-784 * n * t36 + p * ((((378 * n - 3616 * p + 1032) * n + 14272 * t4 - 2924) * n + p * (-24576 * t4 + 7336)) * n + t34) + (-168 * n * t19 + (((-4864 * t18 - 8160) * n + ((-10240 * n + t21 + 14848) * p + 9984 * t18 + 13120) * p - 2712) * p + t18 * (1264 * t18 - 684)) * p + 1248) * p + t18 * ((((-10 + t40) * n - 136) * n + 276) * n + 224) + 384) * t42;
t34 = (0.1e1 / 0.2304e4);
t36 = (0.1e1 / 0.3072e4);
t66 = 0.1e1 / 0.7680e4;
t67 = 0.1e1 / 0.46080e5;
t68 = -1024 * p;
t27 = ((((2565 + t27) * n + 9990) * n + 17510) * n + 13223) * n + ((((-1620 * n - 13500) * n - 38820) * n - 44480) * n + (((4320 * n + 26640) * n + 50260) * n + ((-5760 * n - 23360) * n + (3840 * n + 7680 + t68) * p - 21680) * p + 28220) * p - 16348) * p + 3125;
t69 = -t10 - 5 + t12;
t70 = (-5 + t12) * l + t14;
t71 = cos(t70);
t72 = ((((999 + t33) * n + 2814) * n + 2598) * n + 649) * n + ((((-756 * n - 6012) * n - 11500) * n - 6432) * n + (((2592 * n + 13296) * n + 15500) * n + ((-4224 * n - 12864) * n + (3328 * n + t68 + 4608) * p - 6896) * p + 3924) * p - 676) * p + 27;
t73 = -t10 - 3 + t12;
t74 = (-3 + t12) * l + t14;
t75 = cos(t74);
t76 = ((((261 + t7) * n + 424) * n - 546) * n + 45) * n + ((((-324 * n - 1932) * n - 900) * n + 1436) * n + (((1440 * n + 4656) * n + 60) * n + ((-2944 * n - 4544) * n + (2816 * n + t68 + 1536) * p + 496) * p - 860) * p - 184) * p + 61;
t77 = 1 + n - t9;
t78 = -t77 + t12;
t77 = -t77 * l + t15;
t79 = cos(t77);
t80 = 80;
t81 = ((((23 + t40) * n - t80) * n - 30) * n + 139) * n + ((((-132 * n - 108) * n + 708) * n - 284) * n + (((736 * n - 240) * n - 1548) * n + ((-1920 * n + 1600) * n + (2304 * n + t68 - 1536) * p + 496) * p + 860) * p - 184) * p - 61;
t9 = 1 - n + t9;
t82 = t9 + t12;
t9 = t9 * l + t15;
t15 = sin(t9);
t45 = ((((-52 * n + 420) * n - 1188) * n + 1416) * n + (((352 * n - 2352) * n + 5188) * n + ((-1152 * n + 5568) * n + (n * t45 + t68 - 4608) * p - 6896) * p - 3924) * p - 676) * p + t18 * ((3 * t18 + 82) * n - 90) - 27 * (-1 + t20) * n - 27;
t83 = -t10 + 3 + t12;
t84 = (3 + t12) * l + t14;
t85 = sin(t84);
t86 = -640 * n;
t52 = ((((-20 * n + 420) * n - 3340) * n + 11960) * n + (((160 * n - 2640) * n + 14780) * n + ((7360 + t86) * n + (n * t52 + t68 - 7680) * p - 21680) * p - 28220) * p - 16348) * p + t18 * (((-25 + n) * n + 250) * n - 1250) - 3125 * t51;
t10 = -t10 + 5 + t12;
t12 = (5 + t12) * l + t14;
t14 = sin(t12);
t10 = 1 / t10;
t68 = 1 / t69;
t12 = (t10 * cos(t12));
t69 = (t68 * sin(t70));
t70 = 1 / t73;
t73 = 1 / t83;
t83 = (t73 * cos(t84));
t84 = (t85 * (t73 ^ 2) - (t83 * l));
t74 = (t70 * sin(t74));
t87 = l * t75 - t74;
t82 = 1 / t82;
t78 = 1 / t78;
t9 = (t82 * cos(t9));
t88 = t15 * t82 ^ 2 - t9 * l;
t77 = (t78 * sin(t77));
t89 = l * t79 - t77;
t90 = 0.1e1 / 0.384e3;
t91 = 0.1e1 / 0.768e3;
t92 = 0.1e1 / 0.3840e4;
t93 = 1216 * p;
t8 = (((540 + t8 + t33) * n + (864 * p - 2124) * p + 1212) * n + ((-768 * p + 2784) * p - 3112) * p + 1036) * n + ((1996 - t93) * p - 1292) * p + 256 * t5 + 256;
t33 = 16;
t94 = 256 * p;
t7 = (((207 + t7) * n + 232) * n + 28) * n + (((-216 * n - 906) * n - 578) * n + ((576 * n + 1296) * n + (-608 + t94 + t86) * p + 352) * p - t33) * p - 8;
t3 = ((((-96 * p + t40 + 38) * n + 352 * t4 - 14) * n - 512 * t6 - 50) * n + 196 * t2 - 192 * t3 + 256 * t5 + t53) * t13 * m;
t1 = (-t33 * t20 * t41 + ((204 * n - 880) * n + (-864 * n + 1996 + t93) * p + 1292) * p - 256 * n - 256 * t1 + t19 + 96 * t18 * (1 + t4) + 256) * t48;
t4 = ((((-11 + t46) * n + 6) * n + 12) * n + (((-40 * n + 138) * n - 126) * n + ((192 * n - 528) * n + (-384 * n + 608 + t94) * p + 352) * p + t33) * p - 8) * t42;
t5 = (0.1e1 / 0.64e2);
t6 = (0.1e1 / 0.96e2);
t19 = 1 + t20;
t20 = n * (99 * n + 103) + ((-108 * n - 258) * n + (144 * n + t47 + 168) * p - 130) * p + 27 * t19;
t41 = ((35 + t40) * n - 1) * n + ((-60 * n - 90) * n + (112 * n + t47 + 56) * p + 6) * p - 3;
t19 = -n * (n + 5) + ((-28 * n + 22) * n + (n * t80 + t47 - 56) * p + 6) * p + 3 * t19;
t47 = ((-12 * n + 78) * n + (48 * n + t47 - 168) * p - 130) * p + t18 * (n - 9) - 27 * t51;
t80 = t33 * p;
t40 = n * (14 + t40) + (-n * t53 + t80 - 18) * p + 4;
t2 = (-t2 * t33 - t11 + 3 * n * (n + 1)) * t13 * m;
t11 = ((-8 * n + 18 + t80) * p + t18 + 4 * t51) * t42;
t18 = 4 * p;
t33 = 1 + t46 - t18;
t18 = -t51 - t18;
t46 = e ^ 2;
t51 = t13 * m;
t13 = (l * t17 + (t13 * t16));
t16 = (l * t32 * t30 + (t31 * t60));
t17 = l * t50 + (t48 * t49);
t30 = (l * t39 * t37 + (t38 * t63));
t31 = (l * t44 + (t42 * t43));
t32 = t75 * (t70 ^ 2) + (t74 * l);
t37 = (l * t85 + t83);
t38 = (t79 * (t78 ^ 2) + (t77 * l));
t9 = l * t15 + t9;
t9 = (-t46 * ((((-((t34 * t58 * t13) - (t36 * m * (t21 * t31 + t30 * t35)) + t66 * m * ((t16 * t29) + t28 * t17) - t67 * m * (t23 * (l * t26 * t24 + (t25 * t59)) + t22 * (l * t56 + (t54 * t55)))) * e + t90 * m * (-t81 * t82 * t9 + t38 * t76) - t91 * m * (-(t45 * t73 * t37) + t32 * t72) + t92 * m * (t27 * (t71 * (t68 ^ 2) + (t69 * l)) - (t52 * t10 * (l * t14 + t12)))) * e + (t5 * t3 * t13) - (t6 * m * (t30 * t7 + t4 * t31)) + t90 * m * (t1 * t17 + (t16 * t8))) * e + m * (-(t47 * t73 * t37) + t20 * t32) / 0.48e2 - (m * (-t19 * t82 * t9 + t38 * t41)) / 0.16e2) * e + (m * (t11 * t31 + t30 * t40)) / 0.8e1 - (t2 * t13) / 0.4e1) - (t51 * t13) - (m * (-t18 * t82 * t9 + t33 * t38) * e) / 0.2e1);
Idel = [t46 * (((((-t34 * t58 * t57 + t67 * m * (t23 * ((l * t25 * t24) - t26 * t59) + t22 * ((l * t55) - t54 * t56)) - t66 * m * (t28 * t62 + t29 * t61) + t36 * m * (t21 * t65 + t35 * t64)) * e + t92 * m * ((t52 * (t14 * t10 ^ 2 - t12 * l)) + t27 * t68 * (l * t71 - t69)) - t91 * m * (t72 * t70 * t87 + (t45 * t84)) + t90 * m * (t76 * t78 * t89 + (t81 * t88))) * e + t5 * t3 * t57 + t90 * m * (t1 * t62 + t61 * t8) - t6 * m * (t4 * t65 + t64 * t7)) * e + m * (t20 * t70 * t87 + (t47 * t84)) / 0.48e2 - m * (t41 * t78 * t89 + (t19 * t88)) / 0.16e2) * e - t2 * t57 / 0.4e1 + m * (t11 * t65 + t40 * t64) / 0.8e1) + t51 * t57 + m * (t33 * t78 * t89 + (t18 * t88)) * e / 0.2e1 t9];

t1 = (n - p);
t2 = (p * t1);
t3 = (t2 * n);
t4 = (p ^ 2);
t5 = (t4 ^ 2);
t6 = (p * t4);
t7 = -432 * p;
t8 = 27 * n;
t9 = (delta * m);
t10 = n - t9;
t11 = 2;
t12 = t11 * p;
t13 = t12 - t10;
t14 = -l * t10;
t15 = t12 * l;
t16 = t15 + t14;
t17 = sin(t16);
t18 = n ^ 2;
t19 = t18 ^ 2;
t20 = n * t18;
t1 = (t1 * t6);
t21 = 4096 * p;
t22 = 468952 * p;
t23 = -6 + t12 - t10;
t24 = (-6 + t12) * l + t14;
t25 = 243 * n;
t26 = 59192 * p;
t27 = -4 + t12 - t10;
t28 = (-4 + t12) * l + t14;
t29 = sin(t28);
t30 = 81 * n;
t31 = -3496 * p;
t32 = 1 - p;
t33 = -t11 * t32 - t10;
t34 = -t11 * l * t32 + t14;
t35 = sin(t34);
t36 = 9 * n;
t37 = 1 + p;
t38 = t11 * t37 - t10;
t39 = t11 * l * t37 + t14;
t40 = sin(t39);
t41 = 1792;
t42 = -64 * p;
t43 = 3 * n;
t44 = 4 + t12 - t10;
t45 = (4 + t12) * l + t14;
t46 = sin(t45);
t47 = 1 - n;
t48 = 1280;
t49 = 24;
t50 = 6 + t12 - t10;
t51 = (6 + t12) * l + t14;
t44 = 1 / t44;
t27 = 1 / t27;
t52 = (p * ((((940 * n - 7040 * p - 5340) * n + 25600 * t4 + 14900) * n + (-44800 * t4 - 69900) * p - 21644) * n + t26) + (((((-49 + t43 + t42) * n + 560 * t4 + 310) * n - 2560 * t6 - 960) * n + t4 * (6400 * t4 + 32980) + 1600) * n + t6 * (-8192 * t4 - 86400) - t41) * n + (t4 * (((29696 + t21) * p + 80320) * p + 100800) + 14808) * p + 1536) * t44;
t26 = (p * ((((-32940 * n + 97920 * p - 122700) * n - 143360 * t4 - 177500) * n + (103680 * t4 + 232500) * p - 96740) * n + t26) + (((((-2592 * p + 4347 + t25) * n + 10800 * t4 + 21870) * n - 23040 * t6 + 44840) * n + t4 * (26880 * t4 + 255700) + 39148) * n + t6 * (-16384 * t4 - 234880) + 13016) * n + (t4 * (((-29696 + t21) * p + 80320) * p - 100800) - 14808) * p + 1536) * t27;
t33 = 1 / t33;
t38 = 1 / t38;
t53 = (p * ((((-11898 * n + 40224 * p - 27664) * n - 64448 * t4 - 3612) * n + (49664 * t4 - 800) * p + 6208) * n + t31) + (((((-1080 * p + 1296 + t30) * n + 5616 * t4 + 5172) * n - 14592 * t6 + 1976) * n + t4 * (20224 * t4 + 53556) - 2488) * n + t6 * (-14336 * t4 - 44320) + 464) * n + (t4 * (((-14848 + t21) * p + 13120) * p + 2712) - 1248) * p + 384) * t33;
t31 = (-784 * n * t32 + p * ((((378 * n - 3616 * p + 1032) * n + 14272 * t4 - 2924) * n + p * (-24576 * t4 + 7336)) * n + t31) + (-168 * n * t19 + (((-4864 * t18 - 8160) * n + ((-10240 * n + t21 + 14848) * p + 9984 * t18 + 13120) * p - 2712) * p + t18 * (1264 * t18 - 684)) * p + 1248) * p + t18 * ((((-10 + t36) * n - 136) * n + 276) * n + 224) + 384) * t38;
t32 = 1 / t50;
t23 = 1 / t23;
t32 = (p * ((((750 * n - 6240 * p) * n + 25920 * t4 + 59940) * n + (-53760 * t4 - 274800) * p - 192528) * n + t22) + 46656 * t47 + ((((-t18 * t48 - 178080) * n + ((-6144 * n + t21 + 44544) * p + 3840 * t18 + 192320) * p + 416760) * p + t18 * (240 * t18 + 61620)) * p + t20 * (-t18 * t49 - 9440) + 251232) * p + t18 * ((((-36 + n) * n + 540) * n - 4320) * n + 19440)) * t32;
t21 = (-34560 * t1 * t18 + p * ((((-73710 * n + 194400 * p) * n - 256320 * t4 - 760620) * n + (168960 * t4 + 975480) * p - 745376) * n + t22) + (((((729 * n - 5832 * p + 11178) * n + 19440 * t4 + 66960) * n - 348840 * p + 197580) * n + 681300 * t4 + 295864) * n + t6 * (-18432 * t4 - 591200) + 204576) * n + (t4 * (((-44544 + t21) * p + 192320) * p - 416760) - 251232) * p + 46656) * t23;
t13 = 1 / t13;
t22 = -(0.1e1 / 0.384e3);
t23 = -(0.1e1 / 0.1280e4);
t50 = (0.1e1 / 0.512e3);
t54 = 0.1e1 / 0.7680e4;
t22 = t22 * ((((((261 + t8 + t7) * n + (2736 * p - 2688) * p + 274) * n + ((-8704 * p + 9600) * p + 524) * p - 1476) * n + t4 * ((14592 * p - 13824) * p - 9804) - 90) * n + t6 * ((-12288 * p + 6912) * p + 18560) + 1764) * n + t5 * (4096 * t4 - 9280) - 6136 * t2 + 8844 * t3 - 720) * t13;
t55 = -1024 * p;
t56 = -5 + t12 - t10;
t57 = (-5 + t12) * l + t14;
t58 = -3 + t12 - t10;
t59 = (-3 + t12) * l + t14;
t60 = sin(t59);
t61 = 1 + n - t9;
t62 = t12 - t61;
t61 = -t61 * l + t15;
t63 = sin(t61);
t64 = 80;
t9 = 1 - n + t9;
t65 = t12 + t9;
t9 = t9 * l + t15;
t15 = sin(t9);
t66 = 3 + t12 - t10;
t67 = (3 + t12) * l + t14;
t68 = sin(t67);
t69 = -640 * n;
t10 = 5 + t12 - t10;
t12 = (5 + t12) * l + t14;
t10 = 1 / t10;
t14 = 1 / t56;
t10 = (((((-20 * n + 420) * n - 3340) * n + 11960) * n + (((160 * n - 2640) * n + 14780) * n + ((7360 + t69) * n + (n * t48 + t55 - 7680) * p - 21680) * p - 28220) * p - 16348) * p + t18 * (((-25 + n) * n + 250) * n - 1250) - 3125 * t47) * t10;
t14 = (((((2565 + t25) * n + 9990) * n + 17510) * n + 13223) * n + ((((-1620 * n - 13500) * n - 38820) * n - 44480) * n + (((4320 * n + 26640) * n + 50260) * n + ((-5760 * n - 23360) * n + (3840 * n + 7680 + t55) * p - 21680) * p + 28220) * p - 16348) * p + 3125) * t14;
t25 = 1 / t58;
t48 = 1 / t66;
t41 = (((((-52 * n + 420) * n - 1188) * n + 1416) * n + (((352 * n - 2352) * n + 5188) * n + ((-1152 * n + 5568) * n + (n * t41 + t55 - 4608) * p - 6896) * p - 3924) * p - 676) * p + t18 * ((3 * t18 + 82) * n - 90) - 27 * (-1 + t20) * n - 27) * t48;
t56 = (((((999 + t30) * n + 2814) * n + 2598) * n + 649) * n + ((((-756 * n - 6012) * n - 11500) * n - 6432) * n + (((2592 * n + 13296) * n + 15500) * n + ((-4224 * n - 12864) * n + (3328 * n + t55 + 4608) * p - 6896) * p + 3924) * p - 676) * p + 27) * t25;
t58 = 1 / t62;
t62 = 1 / t65;
t65 = (((((23 + t36) * n - t64) * n - 30) * n + 139) * n + ((((-132 * n - 108) * n + 708) * n - 284) * n + (((736 * n - 240) * n - 1548) * n + ((-1920 * n + 1600) * n + (2304 * n + t55 - 1536) * p + 496) * p + 860) * p - 184) * p - 61) * t62;
t55 = (((((261 + t8) * n + 424) * n - 546) * n + 45) * n + ((((-324 * n - 1932) * n - 900) * n + 1436) * n + (((1440 * n + 4656) * n + 60) * n + ((-2944 * n - 4544) * n + (2816 * n + t55 + 1536) * p + 496) * p - 860) * p - 184) * p + 61) * t58;
t66 = -(0.5e1 / 0.768e3);
t70 = 0.1e1 / 0.768e3;
t71 = 0.5e1 / 0.384e3;
t72 = 1216 * p;
t73 = 16;
t74 = 256 * p;
t1 = (-t73 * t20 * t37 + ((204 * n - 880) * n + (-864 * n + 1996 + t72) * p + 1292) * p - 256 * n - 256 * t1 + t19 + 96 * t18 * (1 + t4) + 256) * t44;
t7 = ((((540 + t30 + t7) * n + (864 * p - 2124) * p + 1212) * n + ((-768 * p + 2784) * p - 3112) * p + 1036) * n + ((1996 - t72) * p - 1292) * p + 256 * t5 + 256) * t27;
t19 = ((((-11 + t43) * n + 6) * n + 12) * n + (((-40 * n + 138) * n - 126) * n + ((192 * n - 528) * n + (-384 * n + 608 + t74) * p + 352) * p + t73) * p - 8) * t38;
t8 = ((((207 + t8) * n + 232) * n + 28) * n + (((-216 * n - 906) * n - 578) * n + ((576 * n + 1296) * n + (-608 + t74 + t69) * p + 352) * p - t73) * p - 8) * t33;
t27 = -(0.1e1 / 0.24e2);
t30 = (0.1e1 / 0.16e2);
t37 = (0.1e1 / 0.96e2);
t3 = t30 * ((((-96 * p + t36 + 38) * n + 352 * t4 - 14) * n - 512 * t6 - 50) * n + 196 * t2 - 192 * t3 + t49 + 256 * t5) * t13;
t4 = 1 + t20;
t5 = (((-12 * n + 78) * n + (48 * n + t42 - 168) * p - 130) * p + t18 * (n - 9) - 27 * t47) * t48;
t6 = (n * (99 * n + 103) + ((-108 * n - 258) * n + (144 * n + t42 + 168) * p - 130) * p + 27 * t4) * t25;
t4 = (-n * (n + 5) + ((-28 * n + 22) * n + (n * t64 + t42 - 56) * p + 6) * p + 3 * t4) * t62;
t20 = (((35 + t36) * n - 1) * n + ((-60 * n - 90) * n + (112 * n + t42 + 56) * p + 6) * p - 3) * t58;
t25 = -(0.3e1 / 0.16e2);
t42 = t73 * p;
t18 = ((-8 * n + 18 + t42) * p + t18 + 4 * t47) * t38;
t33 = (n * (14 + t36) + (-n * t49 + t42 - 18) * p + 4) * t33;
t36 = (0.1e1 / 0.2e1);
t38 = (0.1e1 / 0.4e1);
t2 = t36 * (-t2 * t73 - t11 + 3 * n * (n + 1)) * t13;
t11 = 4 * p;
t13 = (-t11 - t47) * t62;
t11 = (1 - t11 + t43) * t58;
t16 = cos(t16);
t28 = cos(t28);
t34 = cos(t34);
t39 = cos(t39);
t42 = cos(t45);
t43 = cos(t59);
t44 = cos(t61);
t9 = cos(t9);
t45 = cos(t67);
t9 = ((t36 * (t11 * t44 - t13 * t9)) + ((((((t23 * (t26 * t28 + t52 * t42)) + (t50 * (t31 * t39 + t53 * t34)) + t54 * (t21 * cos(t24) + t32 * cos(t51)) + (t22 * t16)) * e + (t66 * (-t41 * t45 + t56 * t43)) + t70 * (t14 * cos(t57) - t10 * cos(t12)) + t71 * (t55 * t44 - t65 * t9)) * e + (t27 * (t19 * t39 + t8 * t34)) + (t37 * (t1 * t42 + t7 * t28)) + (t3 * t16)) * e + (t25 * (t20 * t44 - t4 * t9)) + (t30 * (t6 * t43 - t5 * t45))) * e + (t38 * (t18 * t39 + t33 * t34)) - (t2 * t16)) * e);
Ie = [t36 * (t11 * t63 - (t13 * t15)) + (((((t23 * (t26 * t29 + t52 * t46) + t50 * (t31 * t40 + t53 * t35) + t54 * (t21 * sin(t24) + t32 * sin(t51)) + t22 * t17) * e - t70 * (-t14 * sin(t57) + t10 * sin(t12)) + t66 * (-t41 * t68 + t56 * t60) + t71 * (-(t65 * t15) + t55 * t63)) * e + t37 * (t1 * t46 + t7 * t29) + t27 * (t19 * t40 + t8 * t35) + t3 * t17) * e + t30 * (-t5 * t68 + t6 * t60) + t25 * (-(t4 * t15) + t20 * t63)) * e + t38 * (t18 * t40 + t33 * t35) - t2 * t17) * e t9];

t1 = n - p;
t2 = p * t1;
t3 = t2 * n;
t4 = -0.29e2 / 0.256e3;
t5 = 0.19e2 / 0.3e1;
t6 = 0.3e1 / 0.16e2;
t7 = p ^ 2;
t8 = t7 ^ 2;
t9 = p * t7;
t4 = (((((t6 * p - 0.3e1 / 0.256e3 * n + t4) * n + (-0.19e2 / 0.16e2 * p + 0.7e1 / 0.6e1) * p - 0.137e3 / 0.1152e4) * n + ((0.34e2 / 0.9e1 * p - 0.25e2 / 0.6e1) * p - 0.131e3 / 0.576e3) * p + 0.41e2 / 0.64e2) * n + ((-p * t5 + 0.6e1) * p + 0.817e3 / 0.192e3) * t7 + 0.5e1 / 0.128e3) * n + ((0.16e2 / 0.3e1 * p - 0.3e1) * p - 0.145e3 / 0.18e2) * t9 - 0.49e2 / 0.64e2) * n + 0.767e3 / 0.288e3 * t2 + t8 * (0.145e3 / 0.36e2 - 0.16e2 / 0.9e1 * t7) - 0.737e3 / 0.192e3 * t3 + 0.5e1 / 0.16e2;
t10 = delta * m;
t11 = l * (-n + t10);
t12 = 0.2e1 * p;
t13 = t12 * l;
t14 = t13 + t11;
t15 = cos(t14);
t16 = n ^ 2;
t17 = t16 ^ 2;
t18 = n * t16;
t1 = t1 * t9;
t19 = 0.4e1 / 0.45e2 * p;
t20 = (((((0.81e2 / 0.5120e4 * n - 0.81e2 / 0.640e3 * p + 0.621e3 / 0.2560e4) * n + 0.27e2 / 0.64e2 * t7 + 0.93e2 / 0.64e2) * n - 0.969e3 / 0.128e3 * p + 0.3293e4 / 0.768e3) * n + 0.3785e4 / 0.256e3 * t7 + 0.36983e5 / 0.5760e4) * n + t9 * (-0.2e1 / 0.5e1 * t7 - 0.3695e4 / 0.288e3) + 0.2131e4 / 0.480e3) * n + (t7 * (((-0.29e2 / 0.30e2 + t19) * p + 0.601e3 / 0.144e3) * p - 0.3473e4 / 0.384e3) - 0.2617e4 / 0.480e3) * p + 0.81e2 / 0.80e2;
t21 = 0.58619e5 / 0.5760e4 * p;
t22 = p * ((((0.135e3 / 0.32e2 * p - 0.819e3 / 0.512e3 * n) * n - 0.89e2 / 0.16e2 * t7 - 0.12677e5 / 0.768e3) * n + (0.11e2 / 0.3e1 * t7 + 0.8129e4 / 0.384e3) * p - 0.23293e5 / 0.1440e4) * n + t21) - 0.3e1 / 0.4e1 * t1 * t16;
t23 = (t12 - 0.6e1) * l + t11;
t24 = -0.12785e5 / 0.384e3;
t25 = -0.729e3 / 0.256e3;
t26 = 0.45e2 / 0.32e2;
t27 = -0.8e1 / 0.15e2 * p;
t24 = (((((-0.81e2 / 0.2560e4 * n + 0.27e2 / 0.80e2 * p - 0.1449e4 / 0.2560e4) * n - t26 * t7 + t25) * n + 0.3e1 * t9 - 0.1121e4 / 0.192e3) * n + t7 * (t24 - 0.7e1 / 0.2e1 * t7) - 0.9787e4 / 0.1920e4) * n + t9 * (0.32e2 / 0.15e2 * t7 + 0.367e3 / 0.12e2) - 0.1627e4 / 0.960e3) * n + (t7 * (((t27 + 0.58e2 / 0.15e2) * p - 0.251e3 / 0.24e2) * p + 0.105e3 / 0.8e1) + 0.617e3 / 0.320e3) * p - 0.1e1 / 0.5e1;
t25 = -0.7399e4 / 0.960e3 * p;
t28 = p * ((((0.549e3 / 0.128e3 * n - 0.51e2 / 0.4e1 * p + 0.2045e4 / 0.128e3) * n + 0.56e2 / 0.3e1 * t7 + 0.8875e4 / 0.384e3) * n + (-0.27e2 / 0.2e1 * t7 - 0.3875e4 / 0.128e3) * p + 0.4837e4 / 0.384e3) * n + t25);
t29 = 4;
t30 = (t12 - t29) * l + t11;
t31 = cos(t30);
t32 = 0.4e1 / 0.3e1 * p;
t33 = (((((0.27e2 / 0.1024e4 * n - 0.45e2 / 0.128e3 * p + 0.27e2 / 0.64e2) * n + 0.117e3 / 0.64e2 * t7 + 0.431e3 / 0.256e3) * n - 0.19e2 / 0.4e1 * t9 + 0.247e3 / 0.384e3) * n + t7 * (0.79e2 / 0.12e2 * t7 + 0.4463e4 / 0.256e3) - 0.311e3 / 0.384e3) * n + t9 * (-0.14e2 / 0.3e1 * t7 - 0.1385e4 / 0.96e2) + 0.29e2 / 0.192e3) * n + (t7 * (((t32 - 0.29e2 / 0.6e1) * p + 0.205e3 / 0.48e2) * p + 0.113e3 / 0.128e3) - 0.13e2 / 0.32e2) * p + 0.1e1 / 0.8e1;
t34 = -0.437e3 / 0.384e3 * p;
t35 = p * ((((0.419e3 / 0.32e2 * p - 0.1983e4 / 0.512e3 * n - 0.1729e4 / 0.192e3) * n - 0.1007e4 / 0.48e2 * t7 - 0.301e3 / 0.256e3) * n + (0.97e2 / 0.6e1 * t7 - 0.25e2 / 0.96e2) * p + 0.97e2 / 0.48e2) * n + t34);
t36 = 0.1e1 - p;
t37 = -0.2e1 * l * t36 + t11;
t38 = cos(t37);
t39 = (-0.7e1 / 0.128e3 * n * t17 + (((-0.19e2 / 0.12e2 * t16 - 0.85e2 / 0.32e2) * n + ((-0.10e2 / 0.3e1 * n + t32 + 0.29e2 / 0.6e1) * p + 0.13e2 / 0.4e1 * t16 + 0.205e3 / 0.48e2) * p - 0.113e3 / 0.128e3) * p + t16 * (0.79e2 / 0.192e3 * t16 - 0.57e2 / 0.256e3)) * p + 0.13e2 / 0.32e2) * p + t16 * ((((0.3e1 / 0.1024e4 * n - 0.5e1 / 0.1536e4) * n - 0.17e2 / 0.384e3) * n + 0.23e2 / 0.256e3) * n + 0.7e1 / 0.96e2) + 0.1e1 / 0.8e1;
t34 = -0.49e2 / 0.192e3 * n * t36 + p * ((((0.63e2 / 0.512e3 * n - 0.113e3 / 0.96e2 * p + 0.43e2 / 0.128e3) * n + 0.223e3 / 0.48e2 * t7 - 0.731e3 / 0.768e3) * n + p * (-0.8e1 * t7 + 0.917e3 / 0.384e3)) * n + t34);
t36 = 0.1e1 + p;
t40 = 0.2e1 * l * t36 + t11;
t41 = cos(t40);
t27 = (((((p / 0.120e3 - n / 0.2560e4 + 0.49e2 / 0.7680e4) * n - 0.7e1 / 0.96e2 * t7 - 0.31e2 / 0.768e3) * n + t9 / 0.3e1 + 0.1e1 / 0.8e1) * n + t7 * (-0.5e1 / 0.6e1 * t7 - 0.1649e4 / 0.384e3) - 0.5e1 / 0.24e2) * n + t9 * (0.16e2 / 0.15e2 * t7 + 0.45e2 / 0.4e1) + 0.7e1 / 0.30e2) * n + (t7 * (((t27 - 0.58e2 / 0.15e2) * p - 0.251e3 / 0.24e2) * p - 0.105e3 / 0.8e1) - 0.617e3 / 0.320e3) * p - 0.1e1 / 0.5e1;
t25 = p * ((((0.11e2 / 0.12e2 * p - 0.47e2 / 0.384e3 * n + 0.89e2 / 0.128e3) * n - 0.10e2 / 0.3e1 * t7 - 0.745e3 / 0.384e3) * n + (0.35e2 / 0.6e1 * t7 + 0.1165e4 / 0.128e3) * p + 0.5411e4 / 0.1920e4) * n + t25);
t42 = (t12 + t29) * l + t11;
t43 = cos(t42);
t44 = 0.1e1 - n;
t19 = ((((-t16 / 0.36e2 - 0.371e3 / 0.96e2) * n + ((-0.2e1 / 0.15e2 * n + t19 + 0.29e2 / 0.30e2) * p + t16 / 0.12e2 + 0.601e3 / 0.144e3) * p + 0.3473e4 / 0.384e3) * p + t16 * (t16 / 0.192e3 + 0.1027e4 / 0.768e3)) * p + t18 * (-t16 / 0.1920e4 - 0.59e2 / 0.288e3) + 0.2617e4 / 0.480e3) * p + t16 * ((((n / 0.46080e5 - 0.1e1 / 0.1280e4) * n + 0.3e1 / 0.256e3) * n - 0.3e1 / 0.32e2) * n + 0.27e2 / 0.64e2);
t21 = p * ((((0.25e2 / 0.1536e4 * n - 0.13e2 / 0.96e2 * p) * n + 0.9e1 / 0.16e2 * t7 + 0.333e3 / 0.256e3) * n + (-0.1145e4 / 0.192e3 - 0.7e1 / 0.6e1 * t7) * p - 0.1337e4 / 0.320e3) * n + t21) + 0.81e2 / 0.80e2 * t44;
t45 = (t12 + 0.6e1) * l + t11;
t46 = 0.4e1 / 0.15e2 * p;
t47 = 0.3e1 / 0.2e1 * n;
t48 = 0.9e1 / 0.8e1 * n;
t49 = ((((0.81e2 / 0.1280e4 * n + 0.171e3 / 0.256e3) * n + 0.333e3 / 0.128e3) * n + 0.1751e4 / 0.384e3) * n + 0.13223e5 / 0.3840e4) * n + ((((-0.27e2 / 0.64e2 * n - 0.225e3 / 0.64e2) * n - 0.647e3 / 0.64e2) * n - 0.139e3 / 0.12e2) * n + (((0.111e3 / 0.16e2 + t48) * n + 0.2513e4 / 0.192e3) * n + ((-t47 - 0.73e2 / 0.12e2) * n + (n - t46 + 0.2e1) * p - 0.271e3 / 0.48e2) * p + 0.1411e4 / 0.192e3) * p - 0.4087e4 / 0.960e3) * p + 0.625e3 / 0.768e3;
t50 = (t12 - 0.5e1) * l + t11;
t51 = ((((-0.27e2 / 0.256e3 * n - 0.333e3 / 0.256e3) * n - 0.469e3 / 0.128e3) * n - 0.433e3 / 0.128e3) * n - 0.649e3 / 0.768e3) * n + ((((0.63e2 / 0.64e2 * n + 0.501e3 / 0.64e2) * n + 0.2875e4 / 0.192e3) * n + 0.67e2 / 0.8e1) * n + (((-0.27e2 / 0.8e1 * n - 0.277e3 / 0.16e2) * n - 0.3875e4 / 0.192e3) * n + ((0.11e2 / 0.2e1 * n + 0.67e2 / 0.4e1) * n + (-0.13e2 / 0.3e1 * n + t32 - 0.6e1) * p + 0.431e3 / 0.48e2) * p - 0.327e3 / 0.64e2) * p + 0.169e3 / 0.192e3) * p - 0.9e1 / 0.256e3;
t52 = (t12 - 0.3e1) * l + t11;
t53 = cos(t52);
t54 = 0.22e2 / 0.3e1;
t55 = 0.8e1 / 0.3e1 * p;
t56 = 0.15e2 / 0.4e1 * n;
t54 = ((((0.9e1 / 0.128e3 * n + 0.87e2 / 0.128e3) * n + 0.53e2 / 0.48e2) * n - 0.91e2 / 0.64e2) * n + 0.15e2 / 0.128e3) * n + ((((-0.27e2 / 0.32e2 * n - 0.161e3 / 0.32e2) * n - 0.75e2 / 0.32e2) * n + 0.359e3 / 0.96e2) * n + (((0.97e2 / 0.8e1 + t56) * n + 0.5e1 / 0.32e2) * n + ((-0.23e2 / 0.3e1 * n - 0.71e2 / 0.6e1) * n + (n * t54 + t29 - t55) * p + 0.31e2 / 0.24e2) * p - 0.215e3 / 0.96e2) * p - 0.23e2 / 0.48e2) * p + 0.61e2 / 0.384e3;
t57 = t13 - (0.1e1 + n - t10) * l;
t58 = cos(t57);
t59 = 0.6e1 * n;
t60 = 0.5e1 * n;
t61 = ((((-0.3e1 / 0.128e3 * n - 0.23e2 / 0.384e3) * n + 0.5e1 / 0.24e2) * n + 0.5e1 / 0.64e2) * n - 0.139e3 / 0.384e3) * n + ((((0.11e2 / 0.32e2 * n + 0.9e1 / 0.32e2) * n - 0.59e2 / 0.32e2) * n + 0.71e2 / 0.96e2) * n + (((-0.23e2 / 0.12e2 * n + 0.5e1 / 0.8e1) * n + 0.129e3 / 0.32e2) * n + ((-0.25e2 / 0.6e1 + t60) * n + (t55 - t59 + t29) * p - 0.31e2 / 0.24e2) * p - 0.215e3 / 0.96e2) * p + 0.23e2 / 0.48e2) * p + 0.61e2 / 0.384e3;
t10 = t13 + (0.1e1 - n + t10) * l;
t13 = cos(t10);
t62 = 0.7e1 / 0.3e1;
t62 = ((((-0.13e2 / 0.192e3 * n + 0.35e2 / 0.64e2) * n - 0.99e2 / 0.64e2) * n + 0.59e2 / 0.32e2) * n + (((0.11e2 / 0.24e2 * n - 0.49e2 / 0.16e2) * n + 0.1297e4 / 0.192e3) * n + ((0.29e2 / 0.4e1 - t47) * n + (n * t62 - t32 - 0.6e1) * p - 0.431e3 / 0.48e2) * p - 0.327e3 / 0.64e2) * p - 0.169e3 / 0.192e3) * p + t16 * ((t16 / 0.256e3 + 0.41e2 / 0.384e3) * n - 0.15e2 / 0.128e3) - 0.9e1 / 0.256e3 * (-0.1e1 + t18) * n - 0.9e1 / 0.256e3;
t63 = (t12 + 0.3e1) * l + t11;
t64 = cos(t63);
t46 = ((((n / 0.192e3 - 0.7e1 / 0.64e2) * n + 0.167e3 / 0.192e3) * n - 0.299e3 / 0.96e2) * n + (((-n / 0.24e2 + 0.11e2 / 0.16e2) * n - 0.739e3 / 0.192e3) * n + ((n / 0.6e1 - 0.23e2 / 0.12e2) * n + (-n / 0.3e1 + t46 + 0.2e1) * p + 0.271e3 / 0.48e2) * p + 0.1411e4 / 0.192e3) * p + 0.4087e4 / 0.960e3) * p + t16 * (((-n / 0.3840e4 + 0.5e1 / 0.768e3) * n - 0.25e2 / 0.384e3) * n + 0.125e3 / 0.384e3) + 0.625e3 / 0.768e3 * t44;
t11 = (t12 + 0.5e1) * l + t11;
t3 = (((0.9e1 / 0.64e2 * n - 0.3e1 / 0.2e1 * p + 0.19e2 / 0.32e2) * n + 0.11e2 / 0.2e1 * t7 - 0.7e1 / 0.32e2) * n - 0.8e1 * t9 - 0.25e2 / 0.32e2) * n + t29 * t8 + 0.49e2 / 0.16e2 * t2 - 0.3e1 * t3 + 0.3e1 / 0.8e1;
t9 = 0.19e2 / 0.6e1 * p;
t8 = (((0.27e2 / 0.128e3 * n - 0.9e1 / 0.8e1 * p + t26) * n + (0.9e1 / 0.4e1 * p - 0.177e3 / 0.32e2) * p + 0.101e3 / 0.32e2) * n + ((0.29e2 / 0.4e1 - t12) * p - 0.389e3 / 0.48e2) * p + 0.259e3 / 0.96e2) * n + ((0.499e3 / 0.96e2 - t9) * p - 0.323e3 / 0.96e2) * p + 0.2e1 / 0.3e1 + 0.2e1 / 0.3e1 * t8;
t26 = 0.20e2 / 0.3e1;
t65 = 0.9e1 / 0.4e1 * n;
t26 = (((-0.9e1 / 0.32e2 * n - 0.69e2 / 0.32e2) * n - 0.29e2 / 0.12e2) * n - 0.7e1 / 0.24e2) * n + (((0.151e3 / 0.16e2 + t65) * n + 0.289e3 / 0.48e2) * n + ((-0.27e2 / 0.2e1 - t59) * n + (n * t26 + t5 - t55) * p - 0.11e2 / 0.3e1) * p + 0.1e1 / 0.6e1) * p + 0.1e1 / 0.12e2;
t5 = (((-n / 0.32e2 + 0.11e2 / 0.96e2) * n - 0.1e1 / 0.16e2) * n - 0.1e1 / 0.8e1) * n + (((0.5e1 / 0.12e2 * n - 0.23e2 / 0.16e2) * n + 0.21e2 / 0.16e2) * n + ((-0.2e1 * n + 0.11e2 / 0.2e1) * n + (n * t29 - t5 - t55) * p - 0.11e2 / 0.3e1) * p - 0.1e1 / 0.6e1) * p + 0.1e1 / 0.12e2;
t1 = ((0.17e2 / 0.32e2 * n - 0.55e2 / 0.24e2) * n + (0.499e3 / 0.96e2 + t9 - t65) * p + 0.323e3 / 0.96e2) * p + t17 / 0.384e3 + 0.2e1 / 0.3e1 - 0.2e1 / 0.3e1 * n - 0.2e1 / 0.3e1 * t1 + t16 * (0.1e1 + t7) / 0.4e1 - t18 * t36 / 0.24e2;
t7 = 0.1e1 + t18;
t9 = 0.3e1 * n;
t17 = n * (0.33e2 / 0.16e2 * n + 0.103e3 / 0.48e2) + ((-t65 - 0.43e2 / 0.8e1) * n + (0.7e1 / 0.2e1 + t9 - t32) * p - 0.65e2 / 0.24e2) * p + 0.9e1 / 0.16e2 * t7;
t18 = t29 * p;
t36 = ((-0.9e1 / 0.16e2 * n - 0.35e2 / 0.16e2) * n + 0.1e1 / 0.16e2) * n + ((0.45e2 / 0.8e1 + t56) * n + (-0.7e1 * n + t18 - 0.7e1 / 0.2e1) * p - 0.3e1 / 0.8e1) * p + t6;
t7 = -n * (n / 0.16e2 + 0.5e1 / 0.16e2) + ((-0.7e1 / 0.4e1 * n + 0.11e2 / 0.8e1) * n + (-0.7e1 / 0.2e1 - t18 + t60) * p + 0.3e1 / 0.8e1) * p + t6 * t7;
t6 = ((n / 0.4e1 - 0.13e2 / 0.8e1) * n + (-n + 0.7e1 / 0.2e1 + t32) * p + 0.65e2 / 0.24e2) * p + t16 * (-n / 0.48e2 + t6) + 0.9e1 / 0.16e2 * t44;
t2 = t2 * t29 - 0.3e1 / 0.4e1 * n * (n + 0.1e1) + 0.1e1 / 0.2e1;
t9 = n * (0.7e1 / 0.4e1 + t48) + (-0.9e1 / 0.4e1 - t9 + t12) * p + 0.1e1 / 0.2e1;
t18 = t44 / 0.2e1;
t16 = (-n + 0.9e1 / 0.4e1 + t12) * p + t16 / 0.8e1 + t18;
t29 = (0.1e1 / 0.2e1 + t47 - t12);
t12 = t12 + t18;
t14 = sin(t14);
t18 = sin(t30);
t30 = sin(t37);
t32 = sin(t40);
t37 = sin(t42);
t40 = sin(t52);
t42 = sin(t57);
t10 = sin(t10);
t44 = sin(t63);
t10 = ((((((-t14 * t4 + (-t20 - t22) * sin(t23) + (-t24 - t28) * t18 + (-t33 - t35) * t30 + (-t39 - t34) * t32 + (-t27 - t25) * t37 + (-t19 - t21) * sin(t45)) * e - t49 * sin(t50) - t51 * t40 - t54 * t42 - t61 * t10 - t62 * t44 - t46 * sin(t11)) * e - t37 * t1 - t14 * t3 - t18 * t8 - t30 * t26 - t32 * t5) * e - t10 * t7 - t40 * t17 - t42 * t36 - t44 * t6) * e - t14 * t2 - t32 * t16 - t30 * t9) * e - t10 * t12 - t42 * t29) * e - t14;
Il = [((((((t15 * t4 + (t20 + t22) * cos(t23) + (t24 + t28) * t31 + (t33 + t35) * t38 + (t39 + t34) * t41 + (t27 + t25) * t43 + (t19 + t21) * cos(t45)) * e + t49 * cos(t50) + t51 * t53 + t54 * t58 + t61 * t13 + t62 * t64 + t46 * cos(t11)) * e + t1 * t43 + t15 * t3 + t26 * t38 + t31 * t8 + t41 * t5) * e + t13 * t7 + t17 * t53 + t36 * t58 + t6 * t64) * e + t15 * t2 + t16 * t41 + t38 * t9) * e + t12 * t13 + t29 * t58) * e + t15 t10];

t1 = (n - p);
t2 = (p ^ 2);
t3 = (p * t2);
t4 = (n ^ 2);
t5 = (t1 * t3);
t6 = 729 * n;
t7 = -18432 * t2;
t8 = 4096 * p;
t9 = 1890832 * p;
t10 = (delta * m);
t11 = n - t10;
t12 = 2;
t13 = t12 * p;
t14 = -7 - t11 + t13;
t15 = -l * t11;
t16 = (-7 + t13) * l + t15;
t17 = 20480 * p;
t18 = 1693760 * p;
t19 = -5 - t11 + t13;
t20 = (-5 + t13) * l + t15;
t21 = sin(t20);
t22 = 81 * n;
t23 = 12288 * p;
t24 = 38304 * p;
t25 = -3 - t11 + t13;
t26 = (-3 + t13) * l + t15;
t27 = sin(t26);
t28 = 6144 * t2;
t29 = 16432 * p;
t30 = p * ((((-5754 * n + 4032 * p + 20044) * n + 18240 * t2 + 12948) * n + (-31488 * t2 + 19860) * p - 6392) * n - t29) + (((((135 * n - 1512 * p + 1035) * n + 5904 * t2 - 1534) * n - 8960 * t3 - 10126) * n + t2 * (2304 * t2 - 43908) + 12495) * n + t3 * (t28 + 23776) + 731) * n + (t2 * (((-t8 + 13824) * p + 3520) * p - 23752) + 11714) * p + 17408;
t31 = 1 + n - t10;
t32 = -t31 + t13;
t33 = t13 * l;
t31 = -t31 * l + t33;
t29 = p * ((((198 * n - 8064 * p + 17636) * n + 30528 * t2 - 18588) * n + (-37632 * t2 + 51396) * p - 39256) * n + t29) + (((((-1224 * p + t22 + 111) * n + 7152 * t2 - 1898) * n - 20224 * t3 + 1070) * n + t2 * (28416 * t2 - 48540) + 10329) * n + t3 * (t7 + 37856) - 12445) * n + (t2 * (((t8 + 13824) * p - 3520) * p - 23752) + 11714) * p - 17408;
t10 = 1 - n + t10;
t34 = t10 + t13;
t10 = t10 * l + t33;
t35 = 3 - t11 + t13;
t36 = (3 + t13) * l + t15;
t37 = sin(t36);
t38 = 5 - t11 + t13;
t39 = (5 + t13) * l + t15;
t40 = sin(t39);
t41 = -24 * p;
t42 = 7 - t11 + t13;
t43 = (7 + t13) * l + t15;
t34 = 1 / t34;
t32 = 1 / t32;
t44 = t34 * sin(t10);
t45 = t32 * sin(t31);
t35 = 1 / t35;
t25 = 1 / t25;
t46 = (p * ((((7230 * n - 45120 * p - 23220) * n + 130880 * t2 + 10100) * n + (-172800 * t2 - 63420) * p + 56888) * n - t24) + (((((39 * n - 712 * p - 437) * n + 5200 * t2 + 1510) * n - 19200 * t3 + 170) * n + t2 * (37120 * t2 + 119980) - 10965) * n + t3 * (-34816 * t2 - 246560) + 10563) * n + (t2 * (((81408 + t23) * p + 169920) * p + 95800) + 202) * p + 32256) * t35;
t23 = (p * ((((25470 * n - 124800 * p + 156740) * n + 253760 * t2 + 170660) * n + (-234240 * t2 - 223980) * p - 19720) * n + t24) + (((((216 * p + t22 - 1161) * n - 6480 * t2 - 21630) * n + 26880 * t3 - 42650) * n + t2 * (-47360 * t2 - 399820) - 7619) * n + t3 * (38912 * t2 + 433120) - 10765) * n + (t2 * (((81408 - t23) * p - 169920) * p + 95800) + 202) * p - 32256) * t25;
t24 = 1 / t38;
t14 = 1 / t14;
t19 = 1 / t19;
t38 = 1 / t42;
t28 = (p * ((((1050 * n - 8640 * p - 18860) * n + 35520 * t2 + 174540) * n + (-72960 * t2 - 763860) * p - 835848) * n + t9) + (((((n - 51 + t41) * n + 240 * t2 + 1110) * n - 1280 * t3 - 13250) * n + t2 * (3840 * t2 + 119940) + 91905) * n + t3 * (-338400 - t28) - 352947) * n + (t2 * (((59904 + t8) * p + 357440) * p + 1110920) + 1664606) * p + 589312) * t38;
t38 = (p * ((((-190350 * n + 613440 * p - 1094220) * n - 955200 * t2 - 2556660) * n + (725760 * t2 + 3691740) * p - 2417464) * n + t18) + (((((-9720 * p + t6 + 22437) * n + 45360 * t2 + 176850) * n - 103680 * t3 + 579510) * n + t2 * (126720 * t2 + 2472180) + 848909) * n + t3 * (-79872 * t2 - 2433760) + 475997) * n + (t2 * (((-216576 + t17) * p + 884800) * p - 1751480) - 668282) * p + 35840) * t19;
t6 = (-34560 * t5 * t4 + p * ((((-98010 * n + 259200 * p) * n - 342720 * t2 - 1979580) * n + (226560 * t2 + 2568900) * p - 2945816) * n + t9) + (((((-5832 * p + 14823 + t6) * n + 19440 * t2 + 121230) * n - 635580 * p + 508350) * n + 1249380 * t2 + 1146889) * n + t3 * (-1091360 + t7) + 1311659) * n + (t2 * (((-59904 + t8) * p + 357440) * p - 1110920) - 1664606) * p + 589312) * t14;
t7 = (p * ((((8910 * n - 63360 * p - 88060) * n + 217920 * t2 + 427620) * n + (-357120 * t2 - 1562700) * p - 970056) * n + t18) + (((((17 * n - 360 * p - 489) * n + 3120 * t2 + 5850) * n - 14080 * t3 - 36890) * n + t2 * (34560 * t2 + 479700) + 125205) * n + t3 * (-43008 * t2 - 1105440) - 192285) * n + (t2 * (((216576 + t17) * p + 884800) * p + 1751480) + 668282) * p + 35840) * t24;
t9 = (0.1e1 / 0.9216e4);
t14 = (0.1e1 / 0.15360e5);
t17 = (0.1e1 / 0.46080e5);
t18 = 1 + p;
t1 = p * t1;
t42 = t1 * n;
t47 = -t11 + t13;
t33 = t15 + t33;
t48 = sin(t33);
t49 = 7680 * n;
t50 = 2048 * p;
t51 = 8640 * n;
t52 = -6 - t11 + t13;
t53 = (-6 + t13) * l + t15;
t54 = 12800;
t55 = -4 - t11 + t13;
t56 = (-4 + t13) * l + t15;
t57 = sin(t56);
t58 = 216 * n;
t59 = 54 * n;
t60 = 1 - p;
t61 = -t12 * t60 - t11;
t62 = -t12 * l * t60 + t15;
t63 = sin(t62);
t64 = 5632;
t65 = t12 * t18 - t11;
t66 = t12 * l * t18 + t15;
t67 = sin(t66);
t68 = 4 - t11 + t13;
t69 = (4 + t13) * l + t15;
t70 = sin(t69);
t71 = 2560;
t72 = -1280 * n;
t12 = t12 * n;
t11 = 6 - t11 + t13;
t13 = (6 + t13) * l + t15;
t15 = 1 / t52;
t11 = 1 / t11;
t11 = (((((-75 + t12) * n + 1160) * n - 9300) * n + 38880) * n + ((((-40 * n + 1240) * n - 14900) * n + 82680) * n + (((320 * n - 7680) * n + 63640) * n + ((t72 + 21120) * n + (n * t71 - t50 - 21760) * p - 90400) * p - 182980) * p - 179756) * p - 68232) * t11;
t15 = (((((486 * n + 7155) * n + 40500) * n + 109600) * n + 140876) * n + ((((-3240 * n - 37800) * n - 158820) * n - 283280) * n + (((74880 + t51) * n + 207560) * n + ((-11520 * n - 65920) * n + (21760 + t49 - t50) * p - 90400) * p + 182980) * p - 179756) * p + 68232) * t15;
t52 = 1 / t55;
t55 = 1 / t68;
t51 = (((((162 * n + 4185) * n + 20910) * n + 35940) * n + 17792) * n + ((((-2160 * n - 29520) * n - 99080) * n - 108060) * n + (((72960 + t51) * n + 151080) * n + ((-15360 * n - 76800) * n + (n * t54 - t8 + 29440) * p - 74880) * p + 78980) * p - 26892) * p - 2584) * t52;
t8 = (((((14 * n - 265) * n + 1970) * n - 6860) * n + 9100) * n + ((((-240 * n + 3760) * n - 21560) * n + 49900) * n + (((1600 * n - 19200) * n + 73560) * n + ((-5120 * n + 40960) * n + (-29440 - t8 + t49) * p - 74880) * p - 78980) * p - 26892) * p + 2584) * t55;
t49 = 1 / t61;
t54 = 1 / t65;
t61 = (((((63 + t59) * n - 2032) * n - 276) * n - 288) * n + ((((2472 - t58) * n + 8884) * n - 3304) * n + (((-576 * n - 11520) * n - 12536) * n + ((3328 * n + 16768) * n + (-4608 * n + t50 - 7936) * p + 5664) * p + 4852) * p - 1092) * p - 2840) * t49;
t50 = (((((30 * n - 153) * n - 20) * n + 1272) * n - 1380) * n + ((((-424 * n + 2008) * n - 804) * n - 6400) * n + (((2240 * n - 8832) * n + 4456) * n + ((-5376 * n + 14976) * n + (n * t64 - t50 - 7936) * p - 5664) * p + 4852) * p + 1092) * p - 2840) * t54;
t47 = 1 / t47;
t64 = -(0.1e1 / 0.1536e4);
t65 = -(0.1e1 / 0.3840e4);
t68 = (0.1e1 / 0.7680e4);
t1 = (((n * (n * (((18 * n - 192 * p + 61) * n + 704 * t2 - 234) * n - 1024 * t3 - 264) + 754)) / 0.384e3 - 0.2e1 / 0.3e1 * p * (n * t4 + t3) + 0.4e1 / 0.3e1 * n * t3 * t18 + 0.61e2 / 0.16e2 * t42 + 0.19e2 / 0.96e2 * t1 + 0.19e2 / 0.16e2) * t47);
t3 = 256 * p;
t71 = -432 * n;
t73 = 768 * p;
t74 = 27 * n;
t75 = 96;
t76 = t75 * p;
t5 = -256 * t5;
t77 = (((-144 * p + t74 + 66) * n + (84 + t76) * p - 214) * n + (-576 * p + 92) * p + 16) * n + ((448 * p + 284) * p - 232) * p - t5 - 307;
t4 = -768 * n * t2 * t18 + (((21 * n - 22) * n - 162) * n + 216) * n + (((-208 * n + 276) * n + 660) * n + ((448 + t3) * p + 672 * t4 - 284) * p - 232) * p + 307;
t78 = 16 * p;
t72 = ((((11 * n - 122) * n + 468) * n - 534) * n + (((-144 * n + 1236) * n - 3140) * n + ((672 * n - 3840) * n + (3392 + t73 + t72) * p + 4356) * p + 1220) * p - 495) * t35;
t3 = ((((810 + t22) * n + 2904) * n + 4394) * n + (((-3204 + t71) * n - 7564) * n + ((864 * n + 4224) * n + (-768 * n + t3 - 1856) * p + 4924) * p - 5644) * p + 2347) * t19;
t5 = ((((n - 26 - t78) * n + (324 + t76) * p + 264) * n + (-1344 * p - 2284) * p - 1250) * n + ((1856 * p + 4924) * p + 5644) * p + t5 + 2347) * t24;
t19 = ((((666 + t74) * n + 1684) * n + 686) * n + (((-3732 + t71) * n - 5572) * n + ((1440 * n + 6336) * n + (-1792 * n + t73 - 3392) * p + 4356) * p - 1220) * p - 495) * t25;
t22 = -(0.1e1 / 0.192e3);
t24 = 6 * n;
t71 = 128 * p;
t73 = -78;
t74 = 192 * n;
t76 = 72 * n;
t79 = 8 * n;
t80 = -24 * n;
t52 = (((333 + t59) * n + 650) * n + ((-876 - t58) * n + (288 * n - t71 + 576) * p - 842) * p + 398) * t52;
t55 = (((-33 + t12) * n + 192) * n + ((276 + t80) * n + (n * t75 - t71 - 576) * p - 842) * p - 398) * t55;
t58 = (n * (n * t73 - 29) + ((324 + t76) * n + (-288 + t71 - t74) * p + 74) * p + 68) * t49;
t59 = (((-42 + t79) * n + 45) * n + ((252 - t76) * n + (-288 - t71 + t74) * p - 74) * p + 68) * t54;
t71 = (0.1e1 / 0.48e2);
t73 = (0.1e1 / 0.96e2);
t42 = ((-(((5 + t24) * n - 25) * n) / 0.16e2 + (2 * t42) + 0.7e1 / 0.8e1) * t47);
t18 = (-t78 + t79) * p + 10 * t18 - 3 * n * (1 - n);
t2 = (5 * n + t41 - 7) * n + 16 * t2 - 10 * t60;
t35 = (n * (n - 9) + (38 + t78 - t79) * p + 22) * t35;
t25 = (n * (9 * n + 29) + (-38 + t78 + t80) * p + 22) * t25;
t41 = 8 * p;
t24 = (7 + t24 - t41) * t49;
t41 = (-7 + t12 - t41) * t54;
t49 = (0.1e1 / 0.4e1);
t12 = ((0.1e1 / 0.2e1 + t12 / 0.2e1) * t47);
t47 = e ^ 2;
t54 = (0.1e1 / 0.8e1);
t20 = cos(t20);
t26 = cos(t26);
t36 = cos(t36);
t39 = cos(t39);
t10 = (t34 * cos(t10));
t31 = (t32 * cos(t31));
t32 = cos(t33);
t33 = cos(t56);
t34 = cos(t62);
t56 = cos(t66);
t60 = cos(t69);
t10 = ((t47 * ((((-(t14 * (t23 * t26 - t46 * t36)) + t17 * (-(t7 * t39) + t6 * cos(t16) + t28 * cos(t43) - (t38 * t20)) + (t9 * (-t10 * t29 + t31 * t30))) * e + (t64 * (t61 * t34 + t50 * t56)) + (t65 * (t51 * t33 - t8 * t60)) + t68 * (t15 * cos(t53) - t11 * cos(t13)) + (t1 * t32)) * e + (t22 * (-t10 * t4 + t31 * t77)) - (t19 * t26) / 0.384e3 + (t3 * t20) / 0.384e3 - (t72 * t36) / 0.384e3 + (t5 * t39) / 0.384e3) * e + (t71 * (t58 * t34 + t59 * t56)) + (t73 * (t52 * t33 - t55 * t60)) + (t42 * t32)) + (t49 * (t24 * t34 - t41 * t56)) + (t12 * t32)) * e + t10 + t31 + (t54 * (-t10 * t2 + t31 * t18 + t25 * t26 + t35 * t36) * t47));
Idldg1 = [(((((-t9 * (t44 * t29 - t45 * t30) + t14 * (-t23 * t27 + t46 * t37) - t17 * (-t28 * sin(t43) + t7 * t40 - t6 * sin(t16) + t38 * t21)) * e + t68 * (t15 * sin(t53) - t11 * sin(t13)) + t65 * (t51 * t57 - t8 * t70) + t64 * (t50 * t67 + t61 * t63) + t1 * t48) * e + t22 * (-t44 * t4 + t45 * t77) - t19 * t27 / 0.384e3 + t3 * t21 / 0.384e3 - t72 * t37 / 0.384e3 + t5 * t40 / 0.384e3) * e + t73 * (t52 * t57 - t55 * t70) + t71 * (t58 * t63 + t59 * t67) + t42 * t48) * t47 + t49 * (t24 * t63 - t41 * t67) + t12 * t48) * e + t44 + t45 - t54 * (-t45 * t18 + t44 * t2 - t25 * t27 - t35 * t37) * t47 t10];

t1 = (n - p);
t2 = (p ^ 2);
t3 = (p * t2);
t4 = (n ^ 2);
t5 = (t3 * t1);
t6 = 4096 * p;
t7 = 2339572 * p;
t8 = (delta * m);
t9 = (n - t8);
t10 = 2;
t11 = t10 * p;
t12 = -7 + t11 - t9;
t13 = -l * t9;
t14 = (-7 + t11) * l + t13;
t15 = 28672 * p;
t16 = 4186444 * p;
t17 = t11 - 5 - t9;
t18 = (t11 - 5) * l + t13;
t19 = cos(t18);
t20 = 507084 * p;
t21 = t11 - 3 - t9;
t22 = (t11 - 3) * l + t13;
t23 = cos(t22);
t24 = 1708 * p;
t25 = p * ((((-37590 * n + 138912 * p - 71560) * n - 241728 * t2 - 30288) * n + (199680 * t2 + 32520) * p - 88) * n + t24) + (((((351 * n - 4968 * p + 3798) * n + 27792 * t2 + 11491) * n - 78592 * t3 + 8932) * n + t2 * (119040 * t2 + 155436) - 99) * n + t3 * (-92160 * t2 - 140512) - 218) * n + (t2 * (((t15 - 62976) * p + 44480) * p - 11480) - 602) * p + 49;
t26 = 1 + n - t8;
t27 = t11 - t26;
t28 = t11 * l;
t26 = -t26 * l + t28;
t29 = 135 * n;
t30 = -3328;
t31 = 1110 * n;
t24 = p * ((((-17952 * p + 4304 + t31) * n + 72768 * t2 - 312) * n + (-115200 * t2 - 1920) * p + t30) * n + t24) + (((((-2232 * p + 96 + t29) * n + 14736 * t2 - 665) * n - 49408 * t3 - 316) * n + t2 * (88320 * t2 + 780) + 1521) * n + t3 * (-79872 * t2 - 37408) - 820) * n + (t2 * (((t15 + 62976) * p + 44480) * p + 11480) + 602) * p + 49;
t8 = 1 - n + t8;
t32 = t11 + t8;
t8 = t8 * l + t28;
t33 = t11 + 3 - t9;
t34 = (t11 + 3) * l + t13;
t35 = cos(t34);
t36 = t11 + 5 - t9;
t37 = (t11 + 5) * l + t13;
t38 = cos(t37);
t39 = 1280;
t40 = -24;
t41 = 7 + t11 - t9;
t42 = (7 + t11) * l + t13;
t33 = 1 / t33;
t21 = 1 / t21;
t43 = (p * ((((-135090 * n + 457440 * p - 564800) * n - 751040 * t2 - 985080) * n + (601600 * t2 + 1426240) * p - 739320) * n + t20) + (((((891 * n - 10584 * p + 15336) * n + 49680 * t2 + 88635) * n - 119040 * t3 + 222920) * n + t2 * (154880 * t2 + 1315820) + 266241) * n + t3 * (-104448 * t2 - 1334560) + 151632) * n + (t2 * (((t15 - 188928) * p + 498880) * p - 678120) - 203670) * p + 35721) * t21;
t20 = (p * ((((11570 * n - 75360 * p - 58680) * n + 233920 * t2 + 166960) * n + (-343040 * t2 - 608120) * p - 274848) * n + t20) + (((((51 * n - 968 * p - 682) * n + 7440 * t2 + 3975) * n - 29440 * t3 - 14040) * n + t2 * (62720 * t2 + 305420) + 34005) * n + t3 * (-67584 * t2 - 660960) - 52038) * n + (t2 * (((t15 + 188928) * p + 498880) * p + 678120) + 203670) * p + 35721) * t33;
t41 = 1 / t41;
t17 = 1 / t17;
t12 = 1 / t12;
t36 = 1 / t36;
t31 = (p * ((((-9120 * p + t31 - 21080) * n + 37440 * t2 + 206160) * n + (-76800 * t2 - 897480) * p - 1042008) * n + t7) + (((((p * t40 + n - 54) * n + 240 * t2 + 1245) * n - t3 * t39 - 15740) * n + t2 * (3840 * t2 + 133620) + 115515) * n + t3 * (-6144 * t2 - 375840) - 468462) * n + (t2 * (((t6 + 62976) * p + 395840) * p + 1298840) + 2185610) * p + 823543) * t41;
t41 = (p * ((((11070 * n - 80160 * p - 123680) * n + 283200 * t2 + 717960) * n + (-483840 * t2 - 2721600) * p - 2184072) * n + t16) + (((((19 * n - 408 * p - 600) * n + 3600 * t2 + 8115) * n - 16640 * t3 - 60880) * n + t2 * (42240 * t2 + 688620) + 270705) * n + t3 * (-55296 * t2 - 1646880) - 683400) * n + (t2 * (((t15 + 314880) * p + 1407680) * p + 3275800) + 2793610) * p + 765625) * t36;
t6 = (-34560 * t5 * t4 + p * ((((-102870 * n + 272160 * p) * n - 360000 * t2 - 2307720) * n + (238080 * t2 + 2999040) * p - 3637136) * n + t7) + (((((729 * n - 5832 * p + 15552) * n + 19440 * t2 + 133785) * n - 702000 * p + 591780) * n + 1381140 * t2 + 1413079) * n + t3 * (-18432 * t2 - 1207520) + 1717148) * n + (t2 * (((t6 - 62976) * p + 395840) * p - 1298840) - 2185610) * p + 823543) * t12;
t7 = (p * ((((-339390 * n + 1015200 * p - 1943640) * n - 1496640 * t2 - 5102160) * n + (1090560 * t2 + 7105800) * p - 6188816) * n + t16) + (((((2187 * n - 21384 * p + 44550) * n + 84240 * t2 + 333855) * n - 172800 * t3 + 1211280) * n + t2 * (195840 * t2 + 4194060) + 2273077) * n + t3 * (-116736 * t2 - 3983840) + 2110210) * n + (t2 * (((t15 - 314880) * p + 1407680) * p - 3275800) - 2793610) * p + 765625) * t17;
t12 = 1 / t27;
t15 = 1 / t32;
t16 = (t12 * cos(t26));
t27 = (t15 * cos(t8));
t32 = (0.1e1 / 0.18432e5);
t44 = -0.1e1 / 0.30720e5;
t45 = 0.1e1 / 0.92160e5;
t46 = n - t11;
t1 = p * t1;
t47 = 9 * n;
t48 = 96 * p;
t49 = t11 - t9;
t28 = t28 + t13;
t50 = cos(t28);
t51 = -1024 * p;
t52 = 4320 * n;
t53 = t11 - 6 - t9;
t54 = (t11 - 6) * l + t13;
t55 = -1536 * p;
t56 = t11 - 4 - t9;
t57 = (t11 - 4) * l + t13;
t58 = cos(t57);
t59 = 1440;
t60 = -3072 * p;
t61 = 1 - p;
t62 = -t10 * t61 - t9;
t63 = -t10 * l * t61 + t13;
t64 = cos(t63);
t65 = 1 + p;
t66 = t10 * t65 - t9;
t10 = t10 * l * t65 + t13;
t67 = cos(t10);
t68 = t11 + 4 - t9;
t69 = (t11 + 4) * l + t13;
t70 = cos(t69);
t9 = t11 + 6 - t9;
t11 = (t11 + 6) * l + t13;
t13 = 1 / t53;
t9 = 1 / t9;
t9 = (660 * n * t4 * t65 + ((t4 * (-40 + n) - 5640) * n + 25080) * n + (((-20 * t4 - 8440) * n + 49780) * n + (((160 * n - 4080) * n + 35900) * n + ((-640 * n + 11200) * n + (n * t39 + t51 - 11520) * p - 50800) * p - 109440) * p - 114768) * p - 46656) * t9;
t13 = (((((243 * n + 3780) * n + 22680) * n + 65300) * n + 89688) * n + ((((-1620 * n - 19980) * n - 89040) * n - 169100) * n + (((t52 + 39600) * n + 116500) * n + ((-5760 * n - 34880) * n + (3840 * n + 11520 + t51) * p - 50800) * p + 109440) * p - 114768) * p + 46656) * t13;
t51 = 1 / t56;
t53 = 1 / t68;
t56 = (((((4 * n - 85) * n + 750) * n - 3490) * n + 8680) * n + ((((-70 * n + 1230) * n - 8410) * n + 26780) * n + (((480 * n - 6480) * n + 30100) * n + ((-1600 * n + 14560) * n + (2560 * n + t55 - 11520) * p - 33720) * p - 48120) * p - 33552) * p - 9216) * t53;
t52 = (((((162 * n + 2295) * n + 11280) * n + 24830) * n + 24872) * n + ((((-1350 * n - 14130) * n - 49370) * n - 69460) * n + (((t52 + 31920) * n + 71060) * n + ((-6720 * n - 31520) * n + (5120 * n + t55 + 11520) * p - 33720) * p + 48120) * p - 33552) * p + 9216) * t51;
t55 = 1 / t62;
t62 = 1 / t66;
t66 = (((((21 * n - 140) * n + 408) * n - 780) * n + 1000) * n + ((((-316 * n + 1884) * n - 4704) * n + 6364) * n + (((1824 * n - 8880) * n + 16108) * n + ((-4992 * n + 17216) * n + (6400 * n + t60 - 11520) * p - 16464) * p - 11328) * p - 3888) * p - 576) * t62;
t29 = (((((t29 + t59) * n + 4652) * n + 5744) * n + 2888) * n + ((((-1404 * n - 10308) * n - 21880) * n - 16292) * n + (((5472 * n + 26352) * n + 33284) * n + ((-10112 * n - 28864) * n + (8960 * n + t60 + 11520) * p - 16464) * p + 11328) * p - 3888) * p + 576) * t55;
t49 = 1 / t49;
t60 = (0.1e1 / 0.1536e4);
t68 = (0.1e1 / 0.1920e4);
t71 = 0.1e1 / 0.7680e4;
t3 = ((t46 * (-176 * t1 * n + (((-t48 + t47 + 35) * n + 352 * t2 + 13) * n - 512 * t3 - 22) * n + 256 * t2 ^ 2 + 44 * t1) * t49) / 0.192e3);
t39 = t39 * p;
t72 = (((81 * n + 420) * n + 536) * n + 198) * n + (((-720 * n - 2244) * n - 1552) * n + ((2208 * n + 3744) * n + (-2816 * n + t39 - 1984) * p + 1076) * p - 252) * p + 25;
t73 = (((33 * n - 64) * n + 60) * n - 54) * n + (((-368 * n + 708) * n - 600) * n + ((n * t59 - 2208) * n + (-2304 * n + t39 + 1984) * p + 1076) * p + 252) * p + 25;
t74 = 16 * p;
t30 = ((((189 * n + 1692) * n + 4822) * n + 5376) * n + (((-1296 * n - 7908) * n - 14088) * n + ((3168 * n + 12000) * n + (n * t30 + t39 - 5952) * p + 10124) * p - 7464) * p + 2025) * t21;
t5 = ((((n - t74 - 28) * n + (t48 + 348) * p + 306) * n + (-p * t59 - 2632) * p - 1556) * n + ((1984 * p + 5644) * p + 6960) * p - 256 * t5 + 3125) * t36;
t36 = ((((13 * n - 168) * n + 858) * n - 2088) * n + (((-176 * n + 1764) * n - 6160) * n + ((864 * n - 5856) * n + (-1792 * n + t39 + 5952) * p + 10124) * p + 7464) * p + 2025) * t33;
t2 = (((81 * t4 + 3318) * n + 5404) * n + (((-432 * n - 3420) * n - 8656) * n + ((-768 * n + 256 * p - 1984) * p + 4512 * n + 5644) * p - 6960) * p + 864 * t4 * (n + t2) + 3125) * t17;
t4 = (0.1e1 / 0.384e3);
t17 = (0.1e1 / 0.768e3);
t39 = 64 * p;
t48 = 144 * n;
t59 = 27 * n;
t75 = -128 * p;
t76 = 48 * n;
t53 = (((-18 + n) * n + 114) * n + ((-12 * n + 150) * n + (-t39 + t76 - 312) * p - 496) * p - 256) * t53;
t39 = (((t59 + 180) * n + 382) * n + ((-108 * n - 474) * n + (-t39 + t48 + 312) * p - 496) * p + 256) * t51;
t51 = (((t59 + 129) * n + 170) * n + ((-t48 - 408) * n + (240 * n + t75 + 312) * p - 248) * p + 64) * t55;
t48 = (((5 * n - 33) * n + 78) * n + ((-t76 + 216) * n + (t48 + t75 - 312) * p - 248) * p - 64) * t62;
t59 = (0.1e1 / 0.48e2);
t75 = (0.1e1 / 0.96e2);
t1 = ((t46 * ((6 * n + 5) * n - 32 * t1) * t49) / 0.16e2);
t76 = 48 * p;
t77 = n * (15 * n + 26) + (-56 * n + t76 - 42) * p + 9;
t76 = n * (7 * n - 16) + (-40 * n + 42 + t76) * p + 9;
t21 = (n * (t47 + 32) + (n * t40 + t74 - 42) * p + 27) * t21;
t33 = (n * (n - 10) + (-8 * n + t74 + 42) * p + 27) * t33;
t40 = (3 * n + 4 * t61) * t55;
t47 = (n - 4 * t65) * t62;
t55 = (0.1e1 / 0.4e1);
t46 = t46 * t49;
t49 = e ^ 2;
t61 = (0.1e1 / 0.2e1);
t18 = sin(t18);
t22 = sin(t22);
t34 = sin(t34);
t37 = sin(t37);
t12 = (t12 * sin(t26));
t8 = (t15 * sin(t8));
t15 = sin(t28);
t26 = sin(t57);
t28 = sin(t63);
t10 = sin(t10);
t57 = sin(t69);
t8 = (-(t61 * (-t8 + t12)) + (t49 * (((((t32 * (t12 * t25 - t8 * t24)) + t44 * (-t20 * t34 + t43 * t22) - t45 * (t6 * sin(t14) - (t7 * t18) - t31 * sin(t42) + (t41 * t37))) * e - (t60 * (t66 * t10 + t29 * t28)) + (t68 * (t52 * t26 + t56 * t57)) - t71 * (t13 * sin(t54) + t9 * sin(t11)) + (t3 * t15)) * e - (t17 * (t2 * t18 - t30 * t22 + t36 * t34 - t5 * t37)) - (t4 * (t12 * t72 - t8 * t73))) * e + (t59 * (t48 * t10 + t51 * t28)) - (t75 * (t39 * t26 + t53 * t57)) - (t1 * t15)) - (t55 * (t47 * t10 + t40 * t28)) + (t46 * t15)) * e - (-(t12 * t77) / 0.16e2 + (t21 * t22) / 0.16e2 - (t33 * t34) / 0.16e2 + (t8 * t76) / 0.16e2) * t49);
Idldg2 = [(t61 * (-t27 + t16)) + (((((t44 * (t20 * t35 - t43 * t23) + t45 * (-t7 * t19 + t6 * cos(t14) + t41 * t38 - t31 * cos(t42)) - (t32 * (t16 * t25 - t27 * t24))) * e + t71 * (t13 * cos(t54) + t9 * cos(t11)) - t68 * (t52 * t58 + t56 * t70) + t60 * (t29 * t64 + t66 * t67) - t3 * t50) * e + t17 * (t2 * t19 - t30 * t23 + t36 * t35 - t5 * t38) + (t4 * (t16 * t72 - t27 * t73))) * e + t75 * (t39 * t58 + t53 * t70) - t59 * (t48 * t67 + t51 * t64) + t1 * t50) * t49 + t55 * (t40 * t64 + t47 * t67) - t46 * t50) * e + (-(t16 * t77) / 0.16e2 + t21 * t23 / 0.16e2 + (t27 * t76) / 0.16e2 - t33 * t35 / 0.16e2) * t49 t8];



if Jacobian2On == true
    
    t1 = n - p;
    t2 = p * t1;
    t3 = (t2 * n);
    t4 = (p ^ 2);
    t5 = t4 ^ 2;
    t6 = (p * t4);
    t7 = -432 * p;
    t8 = 27 * n;
    t9 = (((((261 + t8 + t7) * n + (2736 * p - 2688) * p + 274) * n + ((-8704 * p + 9600) * p + 524) * p - 1476) * n + t4 * ((14592 * p - 13824) * p - 9804) - 90) * n + t6 * ((-12288 * p + 6912) * p + 18560) + 1764) * n + t5 * (4096 * t4 - 9280) - 6136 * t2 + 8844 * t3 - 720;
    t10 = (delta * m);
    t11 = n - t10;
    t12 = 2;
    t13 = t12 * p;
    t14 = t13 - t11;
    t15 = -l * t11;
    t16 = t13 * l;
    t17 = t16 + t15;
    t18 = cos(t17);
    t19 = n ^ 2;
    t20 = t19 ^ 2;
    t21 = n * t19;
    t1 = t6 * t1;
    t22 = 4096 * p;
    t23 = (((((729 * n - 5832 * p + 11178) * n + 19440 * t4 + 66960) * n - 348840 * p + 197580) * n + 681300 * t4 + 295864) * n + t6 * (-18432 * t4 - 591200) + 204576) * n + (t4 * (((t22 - 44544) * p + 192320) * p - 416760) - 251232) * p + 46656;
    t24 = 468952 * p;
    t25 = p * ((((-73710 * n + 194400 * p) * n - 256320 * t4 - 760620) * n + (168960 * t4 + 975480) * p - 745376) * n + t24) - 34560 * t1 * t19;
    t26 = t23 + t25;
    t27 = -6 + t13 - t11;
    t28 = (-6 + t13) * l + t15;
    t29 = cos(t28);
    t30 = 243 * n;
    t31 = (((((-2592 * p + 4347 + t30) * n + 10800 * t4 + 21870) * n - 23040 * t6 + 44840) * n + t4 * (26880 * t4 + 255700) + 39148) * n + t6 * (-16384 * t4 - 234880) + 13016) * n + (t4 * (((t22 - 29696) * p + 80320) * p - 100800) - 14808) * p + 1536;
    t32 = 59192 * p;
    t33 = p * ((((-32940 * n + 97920 * p - 122700) * n - 143360 * t4 - 177500) * n + (103680 * t4 + 232500) * p - 96740) * n + t32);
    t34 = t31 + t33;
    t35 = 4;
    t36 = -t35 + t13 - t11;
    t37 = (-t35 + t13) * l + t15;
    t38 = cos(t37);
    t39 = 81 * n;
    t40 = (((((-1080 * p + 1296 + t39) * n + 5616 * t4 + 5172) * n - 14592 * t6 + 1976) * n + t4 * (20224 * t4 + 53556) - 2488) * n + t6 * (-14336 * t4 - 44320) + 464) * n + (t4 * (((t22 - 14848) * p + 13120) * p + 2712) - 1248) * p + 384;
    t41 = -3496 * p;
    t42 = p * ((((-11898 * n + 40224 * p - 27664) * n - 64448 * t4 - 3612) * n + (49664 * t4 - 800) * p + 6208) * n + t41);
    t43 = t40 + t42;
    t44 = 1 - p;
    t45 = -t12 * t44 - t11;
    t46 = -t12 * l * t44 + t15;
    t47 = cos(t46);
    t48 = 9 * n;
    t49 = (-168 * n * t20 + (((-4864 * t19 - 8160) * n + ((-10240 * n + t22 + 14848) * p + 9984 * t19 + 13120) * p - 2712) * p + t19 * (1264 * t19 - 684)) * p + 1248) * p + t19 * ((((-10 + t48) * n - 136) * n + 276) * n + 224) + 384;
    t41 = -784 * n * t44 + p * ((((378 * n - 3616 * p + 1032) * n + 14272 * t4 - 2924) * n + p * (-24576 * t4 + 7336)) * n + t41);
    t50 = t49 + t41;
    t51 = 1 + p;
    t52 = t12 * t51 - t11;
    t53 = t12 * l * t51 + t15;
    t54 = cos(t53);
    t55 = 1792;
    t56 = 3 * n;
    t57 = -64 * p;
    t58 = (((((-49 + t57 + t56) * n + 560 * t4 + 310) * n - 2560 * t6 - 960) * n + t4 * (6400 * t4 + 32980) + 1600) * n + t6 * (-8192 * t4 - 86400) - t55) * n + (t4 * (((t22 + 29696) * p + 80320) * p + 100800) + 14808) * p + 1536;
    t32 = p * ((((940 * n - 7040 * p - 5340) * n + 25600 * t4 + 14900) * n + (-44800 * t4 - 69900) * p - 21644) * n + t32);
    t59 = t58 + t32;
    t60 = t35 + t13 - t11;
    t61 = (t35 + t13) * l + t15;
    t62 = cos(t61);
    t63 = 1 - n;
    t64 = 1280;
    t22 = ((((-t19 * t64 - 178080) * n + ((-6144 * n + t22 + 44544) * p + 3840 * t19 + 192320) * p + 416760) * p + t19 * (240 * t19 + 61620)) * p + t21 * (-24 * t19 - 9440) + 251232) * p + t19 * ((((-36 + n) * n + 540) * n - 4320) * n + 19440);
    t24 = p * ((((750 * n - 6240 * p) * n + 25920 * t4 + 59940) * n + (-53760 * t4 - 274800) * p - 192528) * n + t24) + 46656 * t63;
    t65 = t22 + t24;
    t66 = 6 + t13 - t11;
    t67 = (6 + t13) * l + t15;
    t68 = cos(t67);
    t69 = l ^ 2;
    t70 = (m ^ 2);
    t71 = (delta ^ 2 * t70);
    t72 = t69 * (t71 + t19);
    t73 = t69 * p;
    t74 = t35 * t73 * (-t10 + n - p);
    t75 = t12 * (t10 * t69 * n + 1);
    t76 = -t75 - t74 + t72;
    t17 = sin(t17);
    t77 = -t69 * t11;
    t78 = 24 * p;
    t79 = 12 * t77;
    t28 = sin(t28);
    t80 = 16;
    t81 = 8 * t77;
    t44 = t44 * t69 * t80 + t72 - t74 - t75 - t81;
    t37 = sin(t37);
    t82 = (n - p) * p;
    t83 = 8 * t73;
    t84 = -t35 * t69 * (t10 * (1 - p) - n + t82 - 1) + t72 - t75 - t83;
    t46 = sin(t46);
    t82 = -t35 * t69 * (t10 * (-1 - p) + n + t82 - 1) + t72 - t75 + t83;
    t53 = sin(t53);
    t81 = t80 * t69 * t51 + t72 - t74 - t75 + t81;
    t61 = sin(t61);
    t67 = sin(t67);
    t83 = 1 / t27;
    t85 = 1 / t66;
    t86 = t83 ^ 2;
    t87 = t85 ^ 2;
    t88 = t65 * t67;
    t89 = t26 * t28;
    t90 = 1 / t45;
    t91 = 1 / t52;
    t92 = t91 ^ 2;
    t93 = t91 * t92;
    t94 = t90 ^ 2;
    t95 = t90 * t94;
    t96 = t43 * t46;
    t97 = t50 * t53;
    t98 = 1 / t60;
    t99 = 1 / t36;
    t100 = t98 ^ 2;
    t101 = t98 * t100;
    t102 = t99 ^ 2;
    t103 = t99 * t102;
    t104 = t34 * t37;
    t105 = t59 * t61;
    t106 = l * t70;
    t107 = 1 / t14;
    t108 = 0.1e1 / 0.3840e4;
    t109 = 0.1e1 / 0.7680e4;
    t110 = t107 ^ 2;
    t111 = t107 * t110;
    t112 = l * t18;
    t113 = t107 * t17;
    t114 = -1024 * p;
    t30 = ((((2565 + t30) * n + 9990) * n + 17510) * n + 13223) * n + ((((-1620 * n - 13500) * n - 38820) * n - 44480) * n + (((4320 * n + 26640) * n + 50260) * n + ((-5760 * n - 23360) * n + (3840 * n + 7680 + t114) * p - 21680) * p + 28220) * p - 16348) * p + 3125;
    t115 = -5 + t13 - t11;
    t116 = (-5 + t13) * l + t15;
    t117 = cos(t116);
    t118 = ((((999 + t39) * n + 2814) * n + 2598) * n + 649) * n + ((((-756 * n - 6012) * n - 11500) * n - 6432) * n + (((2592 * n + 13296) * n + 15500) * n + ((-4224 * n - 12864) * n + (3328 * n + t114 + 4608) * p - 6896) * p + 3924) * p - 676) * p + 27;
    t119 = -3 + t13 - t11;
    t120 = (-3 + t13) * l + t15;
    t121 = cos(t120);
    t122 = ((((261 + t8) * n + 424) * n - 546) * n + 45) * n + ((((-324 * n - 1932) * n - 900) * n + 1436) * n + (((1440 * n + 4656) * n + 60) * n + ((-2944 * n - 4544) * n + (2816 * n + t114 + 1536) * p + 496) * p - 860) * p - 184) * p + 61;
    t123 = 1 + n - t10;
    t124 = t13 - t123;
    t123 = -t123 * l + t16;
    t125 = cos(t123);
    t126 = 80;
    t127 = ((((23 + t48) * n - t126) * n - 30) * n + 139) * n + ((((-132 * n - 108) * n + 708) * n - 284) * n + (((736 * n - 240) * n - 1548) * n + ((-1920 * n + 1600) * n + (2304 * n + t114 - 1536) * p + 496) * p + 860) * p - 184) * p - 61;
    t128 = 1 - n + t10;
    t129 = t13 + t128;
    t16 = t128 * l + t16;
    t128 = cos(t16);
    t55 = ((((-52 * n + 420) * n - 1188) * n + 1416) * n + (((352 * n - 2352) * n + 5188) * n + ((-1152 * n + 5568) * n + (n * t55 + t114 - 4608) * p - 6896) * p - 3924) * p - 676) * p + t19 * ((3 * t19 + 82) * n - 90) - 27 * (-1 + t21) * n - 27;
    t130 = 3 + t13 - t11;
    t131 = (3 + t13) * l + t15;
    t132 = cos(t131);
    t133 = -640 * n;
    t64 = ((((-20 * n + 420) * n - 3340) * n + 11960) * n + (((160 * n - 2640) * n + 14780) * n + ((7360 + t133) * n + (n * t64 + t114 - 7680) * p - 21680) * p - 28220) * p - 16348) * p + t19 * (((-25 + n) * n + 250) * n - 1250) - 3125 * t63;
    t11 = 5 + t13 - t11;
    t13 = (5 + t13) * l + t15;
    t15 = cos(t13);
    t114 = 20 * p;
    t134 = 10 * t77;
    t116 = sin(t116);
    t135 = 12 * p;
    t77 = 6 * t77;
    t136 = t69 * (9 - t135) + t72 - t74 - t75 - t77;
    t120 = sin(t120);
    t137 = n + 1;
    t71 = (1 + t71 + t19) * t69;
    t138 = -t12 * (t69 * (t10 * t137 - n) + 1) - t35 * t73 * (-t10 + n - p + 1) + t71;
    t123 = sin(t123);
    t10 = -t12 * (t69 * (-t10 * t63 + n) + 1) - t35 * t73 * (-t10 + n - p - 1) + t71;
    t16 = sin(t16);
    t71 = t69 * (9 + t135) + t72 - t74 - t75 + t77;
    t73 = sin(t131);
    t13 = sin(t13);
    t77 = 1 / t119;
    t131 = 1 / t130;
    t135 = t77 ^ 2;
    t139 = t77 * t135;
    t140 = t131 ^ 2;
    t141 = t131 * t140;
    t142 = t55 * t73;
    t143 = t118 * t120;
    t144 = 1 / t115;
    t145 = 1 / t11;
    t146 = t144 ^ 2;
    t147 = t145 ^ 2;
    t148 = t64 * t13;
    t149 = t30 * t116;
    t150 = 1 / t124;
    t151 = 1 / t129;
    t152 = t150 ^ 2;
    t153 = t150 * t152;
    t154 = t151 ^ 2;
    t155 = t151 * t154;
    t156 = t127 * t16;
    t157 = t122 * t123;
    t158 = 0.1e1 / 0.192e3;
    t159 = 0.1e1 / 0.384e3;
    t160 = 0.1e1 / 0.768e3;
    t114 = (t108 * t70 * (t148 * (t69 * (25 + t114) + t134 + t72 - t74 - t75) * t145 * t147 - t149 * (t69 * (25 - t114) + t72 - t74 - t75 - t134) * t144 * t146) - t160 * t70 * (-t143 * t136 * t139 + t142 * t71 * t141) - t159 * t70 * ((-t118 * t121 * t135 + t55 * t132 * t140) * l - (t156 * t10 * t155) + (t157 * t138 * t153)) + t158 * t106 * (-t122 * t152 * t125 + (t127 * t154 * t128)) + t106 * (-t30 * t146 * t117 + (t64 * t147 * t15)) / 0.1920e4);
    t3 = (((-96 * p + t48 + 38) * n + 352 * t4 - 14) * n - 512 * t6 - 50) * n + 196 * t2 - 192 * t3 + 256 * t5 + 24;
    t6 = 1216 * p;
    t5 = (((540 + t39 + t7) * n + (864 * p - 2124) * p + 1212) * n + ((-768 * p + 2784) * p - 3112) * p + 1036) * n + ((1996 - t6) * p - 1292) * p + 256 * t5 + 256;
    t7 = 256 * p;
    t8 = (((207 + t8) * n + 232) * n + 28) * n + (((-216 * n - 906) * n - 578) * n + ((576 * n + 1296) * n + (t7 + t133 - 608) * p + 352) * p - t80) * p - 8;
    t7 = (((-11 + t56) * n + 6) * n + 12) * n + (((-40 * n + 138) * n - 126) * n + ((192 * n - 528) * n + (-384 * n + 608 + t7) * p + 352) * p + t80) * p - 8;
    t1 = -t80 * t21 * t51 + ((204 * n - 880) * n + (-864 * n + 1996 + t6) * p + 1292) * p - 256 * n - 256 * t1 + t20 + 96 * t19 * (1 + t4) + 256;
    t4 = t5 * t37;
    t6 = t1 * t61;
    t20 = t8 * t46;
    t39 = t7 * t53;
    t51 = (0.1e1 / 0.32e2);
    t133 = (0.1e1 / 0.48e2);
    t134 = (0.1e1 / 0.96e2);
    t21 = 1 + t21;
    t147 = n * (99 * n + 103) + ((-108 * n - 258) * n + (144 * n + t57 + 168) * p - 130) * p + 27 * t21;
    t48 = ((35 + t48) * n - 1) * n + ((-60 * n - 90) * n + (112 * n + t57 + 56) * p + 6) * p - 3;
    t21 = -n * (n + 5) + ((-28 * n + 22) * n + (n * t126 + t57 - 56) * p + 6) * p + 3 * t21;
    t57 = ((-12 * n + 78) * n + (48 * n + t57 - 168) * p - 130) * p + t19 * (n - 9) - 27 * t63;
    t126 = t57 * t73;
    t161 = t147 * t120;
    t162 = t21 * t16;
    t163 = t48 * t123;
    t164 = 0.1e1 / 0.8e1;
    t165 = 0.1e1 / 0.16e2;
    t2 = -0.3e1 / 0.2e1 * n * t137 + (8 * t2) + 0.1e1;
    t137 = -t12 * t2;
    t166 = t35 * p;
    t167 = n * (0.9e1 / 0.4e1 * n + 0.7e1 / 0.2e1) + (-(6 * n) - 0.9e1 / 0.2e1 + t166) * p + 0.1e1;
    t168 = t35 * t167;
    t19 = (p * t80 - 8 * n + 18) * p + t35 * t63 + t19;
    t35 = (t168 * t46);
    t80 = t19 * t53;
    t56 = 1 + t56 - t166;
    t63 = -t166 - t63;
    t166 = t63 * t16;
    t169 = t56 * t123;
    t44 = ((-(-(t70 * (t89 * (t69 * (36 - t78) + t72 - t74 - t75 - t79) * t83 * t86 + t88 * (t69 * (36 + t78) + t72 - t74 - t75 + t79) * t85 * t87)) / 0.46080e5 - (t70 * (t97 * t82 * t93 + t96 * t84 * t95)) / 0.3072e4 + t109 * t70 * (t105 * t81 * t101 + t104 * t44 * t103) - t106 * (t43 * t94 * t47 + t50 * t92 * t54) / 0.1536e4 + t108 * t106 * (t59 * t100 * t62 + t34 * t102 * t38) - t106 * (t26 * t86 * t29 + t65 * t87 * t68) / 0.23040e5 + t9 * t110 * t70 * ((t113 * t76) / 0.2304e4 + t112 / 0.1152e4)) * e - t114) * e + t158 * t106 * (t1 * t100 * t62 + t5 * t102 * t38) + t159 * t70 * (t6 * t81 * t101 + t4 * t44 * t103) - t133 * t106 * (t8 * t94 * t47 + t7 * t92 * t54) - (t134 * t70 * (t20 * t84 * t95 + t39 * t82 * t93)) - t3 * t110 * t70 * (-t112 * t51 - (t113 * t76) / 0.64e2));
    t69 = (t113 - t112);
    t37 = (-l * t38 * t99 + (t102 * t37));
    t61 = (-l * t62 + (t61 * t98));
    t46 = (-l * t47 * t90 + (t46 * t94));
    t53 = (-l * t54 + (t53 * t91));
    t72 = (-l * t121 * t77 + (t120 * t135));
    t73 = (l * t132 - (t131 * t73));
    t74 = (-l * t125 * t150 + (t123 * t152));
    t16 = t128 * l - t151 * t16;
    t13 = ((((((t159 * t9 * t107 * m * t69 - t109 * m * (t26 * (-l * t29 * t83 + (t28 * t86)) + t65 * t85 * (-l * t68 + (t67 * t85))) + m * (t59 * t98 * t61 + t34 * t37) / 0.1280e4 - m * (t50 * t91 * t53 + t43 * t46) / 0.512e3) * e - t160 * m * (t30 * (-l * t117 * t144 + (t116 * t146)) + (t64 * t145 * (l * t15 - t13 * t145))) + 0.5e1 / 0.768e3 * m * (t55 * t131 * t73 + t118 * t72) - 0.5e1 / 0.384e3 * m * (t127 * t151 * t16 + t122 * t74)) * e - t134 * m * (t1 * t98 * t61 + t5 * t37) - t165 * t3 * t107 * m * t69 + m * (t7 * t91 * t53 + t8 * t46) / 0.24e2) * e - t165 * m * (t57 * t131 * t73 + t147 * t72) + 0.3e1 / 0.16e2 * m * (t21 * t151 * t16 + t48 * t74)) * e - m * ((t19 * t91 * t53) + t168 * t46) / 0.4e1 + t137 * t107 * m * t69 / 0.2e1) * e);
    t13 = (-m * (t63 * t151 * t16 + t56 * t74) / 0.2e1 + t13);
    t16 = (l * m);
    t26 = (e ^ 2);
    t16 = (-t16 * (t169 - t166) * e / 0.2e1 - t26 * ((((-(-(t16 * (t89 + t88)) / 0.46080e5 + t109 * t16 * (t105 + t104) - (t16 * (t97 + t96)) / 0.3072e4 + (t16 * t9 * t17) / 0.2304e4) * e + t108 * t16 * (-t148 + t149) + t159 * t16 * (t157 - t156) - t160 * t16 * (-t142 + t143)) * e - (t134 * t16 * (t20 + t39)) + t159 * t16 * (t6 + t4) + (t16 * t3 * t17) / 0.64e2) * e + (t133 * t16 * (-t126 + t161)) - t165 * t16 * (t163 - t162)) * e + t164 * t16 * (t35 + t80) - (t16 * t137 * t17) / 0.4e1) - (t16 * t17));
    t28 = -(0.1e1 / 0.256e3);
    t34 = (0.5e1 / 0.512e3);
    t37 = -(0.5e1 / 0.192e3);
    t43 = (0.5e1 / 0.96e2);
    t46 = -(0.3e1 / 0.8e1);
    t1 = (((((-t18 * t159 * t9 + t29 * (t109 * t23 + t109 * t25) + t38 * (-t33 / 0.1280e4 - t31 / 0.1280e4) + t47 * (t42 / 0.512e3 + t40 / 0.512e3) + t54 * (t41 / 0.512e3 + t49 / 0.512e3) + t62 * (-t32 / 0.1280e4 - t58 / 0.1280e4) + t68 * (t109 * t22 + t109 * t24)) * e + t117 * t160 * t30 - 0.5e1 / 0.768e3 * t118 * t121 + 0.5e1 / 0.384e3 * t125 * t122 - 0.5e1 / 0.384e3 * t128 * t127 + 0.5e1 / 0.768e3 * t55 * t132 - t15 * t160 * t64) * e + t38 * t134 * t5 + t18 * t165 * t3 - t47 * t8 / 0.24e2 - t54 * t7 / 0.24e2 + t62 * t134 * t1) * e + t121 * t165 * t147 - 0.3e1 / 0.16e2 * t125 * t48 + 0.3e1 / 0.16e2 * t128 * t21 - t132 * t165 * t57) * e + t167 * t47 + t18 * t2 + t54 * t19 / 0.4e1) * e + t125 * t56 / 0.2e1 - (t128 * t63) / 0.2e1;
    t2 = [t13 (t35 * t90) / 0.4e1 + (t80 * t91) / 0.4e1 + (((((t28 * (t104 * t99 + t105 * t98)) + (t34 * (t96 * t90 + t97 * t91)) + (t89 * t83) / 0.1536e4 + (t88 * t85) / 0.1536e4 - 0.5e1 / 0.384e3 * t113 * t9) * e + t158 * (t149 * t144 - t148 * t145) + (t37 * (-t142 * t131 + t143 * t77)) + (t43 * (t157 * t150 - t156 * t151))) * e - t164 * (t20 * t90 + t39 * t91) + (t51 * (t4 * t99 + t6 * t98)) + 0.3e1 / 0.16e2 * t113 * t3) * e + t164 * (-t126 * t131 + t161 * t77) + (t46 * (t163 * t150 - t162 * t151))) * e - (t113 * t137) / 0.2e1 0 0 0 t1];
    t5 = [0 0 0 0 0 0];
    t1 = [t16 t1 0 0 0 -((t169 * t124) / 0.2e1 - (t166 * t129) / 0.2e1) * e - (t17 * t14) + t26 * (((((-(t89 * t27) / 0.46080e5 - (t88 * t66) / 0.46080e5 + t109 * (t104 * t36 + t105 * t60) - (t96 * t45) / 0.3072e4 - (t97 * t52) / 0.3072e4 + (t9 * t14 * t17) / 0.2304e4) * e - t108 * (-t148 * t11 + t149 * t115) - t159 * (t157 * t124 - t156 * t129) + t160 * (t143 * t119 - t142 * t130)) * e + (t134 * (t20 * t45 + t39 * t52)) - t159 * (t4 * t36 + t6 * t60) - (t3 * t14 * t17) / 0.64e2) * e - (t133 * (t161 * t119 - t126 * t130)) + t165 * (t163 * t124 - t162 * t129)) * e - t164 * (t35 * t45 + t80 * t52) + (t137 * t14 * t17) / 0.4e1)];
    I1xx = [[-t106 * t12 * t110 * t18 - (((t44 * e + (t133 * t70 * (-t126 * t71 * t141 + t161 * t136 * t139)) + t106 * (t147 * t135 * t121 - t57 * t140 * t132) / 0.24e2 - t164 * t106 * (t48 * t152 * t125 - (t21 * t154 * t128)) - t165 * t70 * (-t162 * t10 * t155 + t163 * t138 * t153)) * e + t164 * t70 * (t35 * t84 * t95 + t80 * t82 * t93) + t70 * ((t47 * t168 * t94 + t54 * t19 * t92) * l - (t137 * t76 * t111 * t17)) / 0.4e1 - t106 * t137 * t110 * t18 / 0.2e1) * e - t106 * (-t56 * t152 * t125 + (t63 * t154 * t128)) + (t70 * (-t166 * t10 * t155 + t169 * t138 * t153)) / 0.2e1) * e - (t70 * t76 * t111 * t17), t13, 0, 0, 0, t16]; t2; t5; t5; t5; t1;];
    
    t1 = n - p;
    t2 = p * t1;
    t3 = (t2 * n);
    t4 = (p ^ 2);
    t5 = t4 ^ 2;
    t6 = (p * t4);
    t7 = 27 * n;
    t8 = -432 * p;
    t9 = (((((261 + t8 + t7) * n + (2736 * p - 2688) * p + 274) * n + ((-8704 * p + 9600) * p + 524) * p - 1476) * n + t4 * ((14592 * p - 13824) * p - 9804) - 90) * n + t6 * ((-12288 * p + 6912) * p + 18560) + 1764) * n + t5 * (4096 * t4 - 9280) - 6136 * t2 + 8844 * t3 - 720;
    t10 = (l ^ 2);
    t11 = (delta * m);
    t12 = (m ^ 2);
    t13 = n ^ 2;
    t14 = t13 ^ 2;
    t15 = n * t13;
    t16 = (delta ^ 2 * t12);
    t17 = t10 * (t16 + t13);
    t18 = t10 * p;
    t19 = 2;
    t20 = t19 * (t11 * t10 * n + 1);
    t21 = 4 * t18 * (-t11 + n - p);
    t22 = -t21 - t20 + t17;
    t23 = n - t11;
    t24 = t19 * p;
    t25 = t24 - t23;
    t26 = -l * t23;
    t27 = t24 * l;
    t28 = t27 + t26;
    t29 = cos(t28);
    t1 = t6 * t1;
    t30 = 4096 * p;
    t31 = (((((729 * n - 5832 * p + 11178) * n + 19440 * t4 + 66960) * n - 348840 * p + 197580) * n + 681300 * t4 + 295864) * n + t6 * (-18432 * t4 - 591200) + 204576) * n + (t4 * (((t30 - 44544) * p + 192320) * p - 416760) - 251232) * p + 46656;
    t32 = 468952 * p;
    t33 = p * ((((-73710 * n + 194400 * p) * n - 256320 * t4 - 760620) * n + (168960 * t4 + 975480) * p - 745376) * n + t32) - 34560 * t1 * t13;
    t34 = t31 + t33;
    t35 = -t10 * t23;
    t36 = 24 * p;
    t37 = 12 * t35;
    t38 = -6 + t24 - t23;
    t39 = (-6 + t24) * l + t26;
    t40 = cos(t39);
    t41 = 243 * n;
    t42 = (((((-2592 * p + 4347 + t41) * n + 10800 * t4 + 21870) * n - 23040 * t6 + 44840) * n + t4 * (26880 * t4 + 255700) + 39148) * n + t6 * (-16384 * t4 - 234880) + 13016) * n + (t4 * (((t30 - 29696) * p + 80320) * p - 100800) - 14808) * p + 1536;
    t43 = 59192 * p;
    t44 = p * ((((-32940 * n + 97920 * p - 122700) * n - 143360 * t4 - 177500) * n + (103680 * t4 + 232500) * p - 96740) * n + t43);
    t45 = t42 + t44;
    t46 = -1 + p;
    t47 = 16;
    t48 = 8 * t35;
    t49 = -t10 * t46 * t47 + t17 - t20 - t21 - t48;
    t50 = -4 + t24 - t23;
    t51 = (-4 + t24) * l + t26;
    t52 = cos(t51);
    t53 = 81 * n;
    t54 = (((((-1080 * p + 1296 + t53) * n + 5616 * t4 + 5172) * n - 14592 * t6 + 1976) * n + t4 * (20224 * t4 + 53556) - 2488) * n + t6 * (-14336 * t4 - 44320) + 464) * n + (t4 * (((t30 - 14848) * p + 13120) * p + 2712) - 1248) * p + 384;
    t55 = -3496 * p;
    t56 = p * ((((-11898 * n + 40224 * p - 27664) * n - 64448 * t4 - 3612) * n + (49664 * t4 - 800) * p + 6208) * n + t55);
    t57 = t54 + t56;
    t58 = (n - p) * p;
    t59 = 8 * t18;
    t60 = t17 - 4 * t10 * (t11 * (1 - p) - n + t58 - 1) - t59 - t20;
    t61 = t19 * t46 - t23;
    t62 = t19 * l * t46 + t26;
    t63 = cos(t62);
    t64 = 9 * n;
    t65 = (-168 * n * t14 + (((-4864 * t13 - 8160) * n + ((-10240 * n + t30 + 14848) * p + 9984 * t13 + 13120) * p - 2712) * p + t13 * (1264 * t13 - 684)) * p + 1248) * p + t13 * ((((-10 + t64) * n - 136) * n + 276) * n + 224) + 384;
    t46 = 784 * n * t46 + p * ((((378 * n - 3616 * p + 1032) * n + 14272 * t4 - 2924) * n + p * (-24576 * t4 + 7336)) * n + t55);
    t55 = t65 + t46;
    t58 = t17 - 4 * t10 * (t11 * (-1 - p) + n + t58 - 1) + t59 - t20;
    t59 = 1 + p;
    t66 = t19 * t59 - t23;
    t67 = t19 * l * t59 + t26;
    t68 = cos(t67);
    t69 = 1792;
    t70 = -64 * p;
    t71 = 3 * n;
    t72 = (((((-49 + t71 + t70) * n + 560 * t4 + 310) * n - 2560 * t6 - 960) * n + t4 * (6400 * t4 + 32980) + 1600) * n + t6 * (-8192 * t4 - 86400) - t69) * n + (t4 * (((t30 + 29696) * p + 80320) * p + 100800) + 14808) * p + 1536;
    t43 = p * ((((940 * n - 7040 * p - 5340) * n + 25600 * t4 + 14900) * n + (-44800 * t4 - 69900) * p - 21644) * n + t43);
    t73 = t72 + t43;
    t48 = t47 * t10 * t59 + t17 - t20 - t21 + t48;
    t74 = 4 + t24 - t23;
    t75 = (4 + t24) * l + t26;
    t76 = cos(t75);
    t77 = 1 - n;
    t78 = 1280;
    t30 = ((((-t13 * t78 - 178080) * n + ((-6144 * n + t30 + 44544) * p + 3840 * t13 + 192320) * p + 416760) * p + t13 * (240 * t13 + 61620)) * p + t15 * (-24 * t13 - 9440) + 251232) * p + t13 * ((((-36 + n) * n + 540) * n - 4320) * n + 19440);
    t32 = p * ((((750 * n - 6240 * p) * n + 25920 * t4 + 59940) * n + (-53760 * t4 - 274800) * p - 192528) * n + t32) + 46656 * t77;
    t79 = t30 + t32;
    t80 = 6 + t24 - t23;
    t81 = (6 + t24) * l + t26;
    t82 = cos(t81);
    t28 = sin(t28);
    t39 = sin(t39);
    t51 = sin(t51);
    t62 = sin(t62);
    t67 = sin(t67);
    t75 = sin(t75);
    t81 = sin(t81);
    t83 = 1 / t50;
    t84 = 1 / t74;
    t85 = t84 ^ 2;
    t86 = t84 * t85;
    t87 = t83 ^ 2;
    t88 = t83 * t87;
    t89 = t73 * t76;
    t90 = t45 * t52;
    t91 = 1 / t38;
    t92 = 1 / t80;
    t93 = t91 ^ 2;
    t94 = t92 ^ 2;
    t95 = t34 * t40;
    t96 = t79 * t82;
    t97 = 1 / t66;
    t98 = 1 / t61;
    t99 = t97 ^ 2;
    t100 = t97 * t99;
    t101 = t98 ^ 2;
    t102 = t98 * t101;
    t103 = t55 * t68;
    t104 = t57 * t63;
    t105 = l * t12;
    t106 = 1 / t25;
    t107 = 0.1e1 / 0.3840e4;
    t108 = 0.1e1 / 0.1536e4;
    t109 = t106 ^ 2;
    t110 = t106 * t109;
    t111 = l * t28;
    t112 = t106 * t29;
    t113 = -1024 * p;
    t41 = ((((2565 + t41) * n + 9990) * n + 17510) * n + 13223) * n + ((((-1620 * n - 13500) * n - 38820) * n - 44480) * n + (((4320 * n + 26640) * n + 50260) * n + ((-5760 * n - 23360) * n + (3840 * n + 7680 + t113) * p - 21680) * p + 28220) * p - 16348) * p + 3125;
    t114 = 20 * p;
    t115 = 10 * t35;
    t116 = -5 + t24 - t23;
    t117 = (-5 + t24) * l + t26;
    t118 = cos(t117);
    t119 = ((((999 + t53) * n + 2814) * n + 2598) * n + 649) * n + ((((-756 * n - 6012) * n - 11500) * n - 6432) * n + (((2592 * n + 13296) * n + 15500) * n + ((-4224 * n - 12864) * n + (3328 * n + t113 + 4608) * p - 6896) * p + 3924) * p - 676) * p + 27;
    t120 = 12 * p;
    t35 = 6 * t35;
    t121 = t10 * (9 - t120) + t17 - t20 - t21 - t35;
    t122 = -3 + t24 - t23;
    t123 = (-3 + t24) * l + t26;
    t124 = cos(t123);
    t125 = ((((261 + t7) * n + 424) * n - 546) * n + 45) * n + ((((-324 * n - 1932) * n - 900) * n + 1436) * n + (((1440 * n + 4656) * n + 60) * n + ((-2944 * n - 4544) * n + (2816 * n + t113 + 1536) * p + 496) * p - 860) * p - 184) * p + 61;
    t126 = n + 1;
    t16 = (1 + t16 + t13) * t10;
    t127 = -t19 * (t10 * (t11 * t126 - n) + 1) + t16 - 4 * t18 * (-t11 + n - p + 1);
    t128 = 1 + n - t11;
    t129 = t24 - t128;
    t128 = -t128 * l + t27;
    t130 = cos(t128);
    t131 = 80;
    t132 = ((((23 + t64) * n - t131) * n - 30) * n + 139) * n + ((((-132 * n - 108) * n + 708) * n - 284) * n + (((736 * n - 240) * n - 1548) * n + ((-1920 * n + 1600) * n + (2304 * n + t113 - 1536) * p + 496) * p + 860) * p - 184) * p - 61;
    t16 = -t19 * (t10 * (-t11 * t77 + n) + 1) + t16 - 4 * t18 * (-t11 + n - p - 1);
    t11 = 1 - n + t11;
    t18 = t24 + t11;
    t11 = t11 * l + t27;
    t27 = cos(t11);
    t69 = ((((-52 * n + 420) * n - 1188) * n + 1416) * n + (((352 * n - 2352) * n + 5188) * n + ((-1152 * n + 5568) * n + (n * t69 + t113 - 4608) * p - 6896) * p - 3924) * p - 676) * p + t13 * ((3 * t13 + 82) * n - 90) - 27 * (-1 + t15) * n - 27;
    t35 = t10 * (9 + t120) + t17 - t20 - t21 + t35;
    t120 = 3 + t24 - t23;
    t133 = (3 + t24) * l + t26;
    t134 = cos(t133);
    t135 = -640 * n;
    t78 = ((((-20 * n + 420) * n - 3340) * n + 11960) * n + (((160 * n - 2640) * n + 14780) * n + ((t135 + 7360) * n + (n * t78 + t113 - 7680) * p - 21680) * p - 28220) * p - 16348) * p + t13 * (((-25 + n) * n + 250) * n - 1250) - 3125 * t77;
    t23 = 5 + t24 - t23;
    t24 = (5 + t24) * l + t26;
    t26 = cos(t24);
    t113 = sin(t117);
    t117 = sin(t123);
    t123 = sin(t128);
    t11 = sin(t11);
    t128 = sin(t133);
    t24 = sin(t24);
    t133 = 1 / t129;
    t136 = 1 / t18;
    t137 = t136 ^ 2;
    t138 = t136 * t137;
    t139 = t133 ^ 2;
    t140 = t133 * t139;
    t141 = 1 / t120;
    t142 = 1 / t122;
    t143 = t141 ^ 2;
    t144 = t141 * t143;
    t145 = t142 ^ 2;
    t146 = t142 * t145;
    t147 = t132 * t27;
    t148 = t125 * t130;
    t149 = 1 / t116;
    t150 = 1 / t23;
    t151 = t149 ^ 2;
    t152 = t150 ^ 2;
    t153 = t41 * t118;
    t154 = t78 * t26;
    t155 = t119 * t124;
    t156 = t69 * t134;
    t157 = 0.1e1 / 0.192e3;
    t158 = 0.1e1 / 0.384e3;
    t159 = 0.1e1 / 0.768e3;
    t114 = (t107 * t12 * ((t154 * (t10 * (25 + t114) + t115 + t17 - t20 - t21) * t150 * t152) - t153 * (t10 * (25 - t114) + t17 - t20 - t21 - t115) * t149 * t151) + t157 * t105 * (-t132 * t137 * t11 + t125 * t139 * t123) + t158 * t12 * (((-t119 * t117 * t145 + t69 * t128 * t143) * l) + (t147 * t16 * t138) - t148 * t127 * t140) - (t105 * (-t41 * t151 * t113 + t78 * t152 * t24)) / 0.1920e4 + t159 * t12 * (t155 * t121 * t146 - t156 * t35 * t144));
    t3 = (((-96 * p + t64 + 38) * n + 352 * t4 - 14) * n - 512 * t6 - 50) * n + 196 * t2 - 192 * t3 + 256 * t5 + 24;
    t6 = 1216 * p;
    t5 = (((540 + t8 + t53) * n + (864 * p - 2124) * p + 1212) * n + ((-768 * p + 2784) * p - 3112) * p + 1036) * n + ((1996 - t6) * p - 1292) * p + 256 * t5 + 256;
    t8 = 256 * p;
    t7 = (((207 + t7) * n + 232) * n + 28) * n + (((-216 * n - 906) * n - 578) * n + ((576 * n + 1296) * n + (-608 + t135 + t8) * p + 352) * p - t47) * p - 8;
    t8 = (((-11 + t71) * n + 6) * n + 12) * n + (((-40 * n + 138) * n - 126) * n + ((192 * n - 528) * n + (-384 * n + 608 + t8) * p + 352) * p + t47) * p - 8;
    t1 = -t47 * t15 * t59 + ((204 * n - 880) * n + (-864 * n + 1996 + t6) * p + 1292) * p - 256 * n - 256 * t1 + t14 + 96 * t13 * (1 + t4) + 256;
    t4 = (t1 * t76);
    t6 = (t5 * t52);
    t14 = (t8 * t68);
    t53 = (t7 * t63);
    t59 = -(0.1e1 / 0.48e2);
    t115 = (0.1e1 / 0.96e2);
    t135 = (0.1e1 / 0.32e2);
    t15 = 1 + t15;
    t151 = n * (99 * n + 103) + ((-108 * n - 258) * n + (144 * n + t70 + 168) * p - 130) * p + 27 * t15;
    t64 = ((35 + t64) * n - 1) * n + ((-60 * n - 90) * n + (112 * n + t70 + 56) * p + 6) * p - 3;
    t15 = -n * (n + 5) + ((-28 * n + 22) * n + (n * t131 + t70 - 56) * p + 6) * p + 3 * t15;
    t70 = ((-12 * n + 78) * n + (48 * n + t70 - 168) * p - 130) * p + t13 * (n - 9) - 27 * t77;
    t131 = (t151 * t124);
    t160 = t70 * t134;
    t161 = t15 * t27;
    t162 = t64 * t130;
    t2 = 0.3e1 / 0.2e1 * n * t126 - (8 * t2) - 0.1e1;
    t126 = t19 * t2;
    t163 = 4 * p;
    t164 = n * (-0.9e1 / 0.4e1 * n - 0.7e1 / 0.2e1) + ((6 * n) + 0.9e1 / 0.2e1 - t163) * p - 0.1e1;
    t165 = -0.4e1 * t164;
    t13 = (p * t47 - 8 * n + 18) * p + t13 + 4 * t77;
    t47 = (t13 * t68);
    t166 = t165 * t63;
    t71 = 1 + t71 - t163;
    t77 = -t163 - t77;
    t163 = t77 * t27;
    t167 = t71 * t130;
    t10 = ((-(t12 * (t89 * t48 * t86 + t90 * t49 * t88) / 0.7680e4 - t12 * (t96 * (t10 * (36 + t36) + t17 - t20 - t21 + t37) * t92 * t94 + t95 * (t10 * (36 - t36) + t17 - t20 - t21 - t37) * t91 * t93) / 0.46080e5 - t12 * (t103 * t58 * t100 + t104 * t60 * t102) / 0.3072e4 + (t105 * (t34 * t93 * t39 + t79 * t94 * t81)) / 0.23040e5 - t107 * t105 * (t45 * t87 * t51 + t73 * t85 * t75) + t108 * t105 * (t57 * t101 * t62 + t55 * t99 * t67) + t9 * t109 * t12 * (t112 * t22 / 0.2304e4 - t111 / 0.1152e4)) * e - t114) * e - t157 * t105 * (t1 * t85 * t75 + t5 * t87 * t51) + t158 * t12 * (t4 * t48 * t86 + t6 * t49 * t88) - (t115 * t12 * (t14 * t58 * t100 + t53 * t60 * t102)) - (t59 * t105 * (t7 * t101 * t62 + t8 * t99 * t67)) - t3 * t109 * t12 * ((t111 * t135) - t112 * t22 / 0.64e2));
    t17 = (t112 + t111);
    t20 = ((l * t51 * t83) + t52 * t87);
    t21 = ((l * t75) + t76 * t84);
    t36 = ((l * t62 * t98) + t101 * t63);
    t37 = ((l * t67) + t68 * t97);
    t48 = ((l * t117 * t142) + t124 * t145);
    t49 = ((t128 * l) + t134 * t141);
    t52 = (l * t123 * t133) + t130 * t139;
    t27 = t11 * l + t136 * t27;
    t17 = (-m * (-(t77 * t136 * t27) + t52 * t71) / 0.2e1 + (((((t158 * t9 * t106 * m * t17 + m * (t73 * t84 * t21 + t20 * t45) / 0.1280e4 - m * (t55 * t97 * t37 + t36 * t57) / 0.512e3 - m * (t34 * ((l * t39 * t91) + t40 * t93) + t79 * t92 * ((l * t81) + t82 * t92)) / 0.7680e4) * e + t159 * m * ((t78 * (t24 * l * t150 + t152 * t26)) - t41 * t149 * ((l * t113) + t118 * t149)) + 0.5e1 / 0.768e3 * m * (-t69 * t141 * t49 + t119 * t48) - 0.5e1 / 0.384e3 * m * (-(t132 * t136 * t27) + t125 * t52)) * e - t115 * m * (t1 * t84 * t21 + t20 * t5) + m * (t8 * t97 * t37 + t36 * t7) / 0.24e2 - t3 * t106 * m * t17 / 0.16e2) * e + 0.3e1 / 0.16e2 * m * (-(t15 * t136 * t27) + t52 * t64) - m * (-t70 * t141 * t49 + t151 * t48) / 0.16e2) * e - m * ((t13 * t97 * t37) + t165 * t36) / 0.4e1 + t126 * t106 * m * t17 / 0.2e1) * e);
    t20 = (l * m);
    t21 = (e ^ 2);
    t20 = (-t20 * (t167 - t163) * e / 0.2e1 - t21 * ((((-(-t20 * (t96 + t95) / 0.46080e5 + t20 * (t90 + t89) / 0.7680e4 - t20 * (t104 + t103) / 0.3072e4 + t20 * t9 * t29 / 0.2304e4) * e + t107 * t20 * (-t154 + t153) + t158 * t20 * (t148 - t147) - t159 * t20 * (-t156 + t155)) * e - (t115 * t20 * (t53 + t14)) + t158 * t20 * (t4 + t6) + t20 * t3 * t29 / 0.64e2) * e - t59 * t20 * (t131 - t160) - t20 * (t162 - t161) / 0.16e2) * e + t20 * (t166 + t47) / 0.8e1 - t20 * t126 * t29 / 0.4e1) - t20 * t29);
    t26 = -(0.1e1 / 0.256e3);
    t27 = (0.5e1 / 0.512e3);
    t34 = (0.5e1 / 0.96e2);
    t36 = (0.5e1 / 0.192e3);
    t37 = -(0.3e1 / 0.8e1);
    t1 = (((((t28 * t158 * t9 + t39 * (-t33 / 0.7680e4 - t31 / 0.7680e4) + t51 * (t44 / 0.1280e4 + t42 / 0.1280e4) + t62 * (-t56 / 0.512e3 - t54 / 0.512e3) + t67 * (-t46 / 0.512e3 - t65 / 0.512e3) + t75 * (t43 / 0.1280e4 + t72 / 0.1280e4) + t81 * (-t32 / 0.7680e4 - t30 / 0.7680e4)) * e - t113 * t159 * t41 + 0.5e1 / 0.384e3 * t11 * t132 + 0.5e1 / 0.768e3 * t119 * t117 - 0.5e1 / 0.384e3 * t123 * t125 - 0.5e1 / 0.768e3 * t69 * t128 + t24 * t159 * t78) * e - (t28 * t3) / 0.16e2 - (t51 * t115 * t5) + (t62 * t7) / 0.24e2 + (t67 * t8) / 0.24e2 - (t75 * t115 * t1)) * e - 0.3e1 / 0.16e2 * t11 * t15 - (t117 * t151) / 0.16e2 + 0.3e1 / 0.16e2 * t123 * t64 + (t128 * t70) / 0.16e2) * e + t164 * t62 + (t2 * t28) - (t67 * t13) / 0.4e1) * e + (t11 * t77) / 0.2e1 - (t123 * t71) / 0.2e1;
    t2 = [t17 t166 * t98 / 0.4e1 + (t47 * t97) / 0.4e1 + ((((t108 * (t95 * t91 + t96 * t92) + t26 * (t90 * t83 + t89 * t84) + t27 * (t103 * t97 + t104 * t98) - 0.5e1 / 0.384e3 * t112 * t9) * e - t157 * (-t153 * t149 + (t154 * t150)) + t34 * (t148 * t133 - (t147 * t136)) + t36 * (t156 * t141 - t155 * t142)) * e + (t135 * (t4 * t84 + t6 * t83)) - (t14 * t97) / 0.8e1 - (t53 * t98) / 0.8e1 + 0.3e1 / 0.16e2 * t112 * t3) * e + t37 * (t162 * t133 - (t161 * t136)) + (t131 * t142) / 0.8e1 - t160 * t141 / 0.8e1) * e - t112 * t126 / 0.2e1 0 0 0 t1];
    t5 = [0 0 0 0 0 0];
    t1 = [t20 t1 0 0 0 -(t167 * t129 / 0.2e1 - (t163 * t18) / 0.2e1) * e + t21 * (((((-t95 * t38 / 0.46080e5 - t96 * t80 / 0.46080e5 + t90 * t50 / 0.7680e4 + t89 * t74 / 0.7680e4 - t103 * t66 / 0.3072e4 - t104 * t61 / 0.3072e4 + t9 * t25 * t29 / 0.2304e4) * e - t107 * (t153 * t116 - (t154 * t23)) - t158 * (t148 * t129 - (t147 * t18)) + t159 * (-t156 * t120 + t155 * t122)) * e + (t115 * (t14 * t66 + t53 * t61)) - t158 * (t4 * t74 + t6 * t50) - t3 * t25 * t29 / 0.64e2) * e + t59 * (-t160 * t120 + (t131 * t122)) + t162 * t129 / 0.16e2 - (t161 * t18) / 0.16e2) * e - t166 * t61 / 0.8e1 - (t47 * t66) / 0.8e1 + t126 * t25 * t29 / 0.4e1) - t29 * t25];
    I2xx = [[(t105 * t19 * t109 * t28) - (((t10 * e - t59 * t12 * ((t131 * t121 * t146) - t160 * t35 * t144) - t12 * (t162 * t127 * t140 - (t161 * t16 * t138)) / 0.16e2 - (t105 * (t151 * t145 * t117 - t70 * t143 * t128)) / 0.24e2 + (t105 * (-t15 * t137 * t11 + t64 * t139 * t123)) / 0.8e1) * e - t12 * ((t62 * t165 * t101 + (t67 * t13 * t99)) * l + t126 * t22 * t110 * t29) / 0.4e1 + t12 * ((t47 * t58 * t100) + t166 * t60 * t102) / 0.8e1 + (t105 * t126 * t109 * t28) / 0.2e1) * e + t12 * (t167 * t127 * t140 - (t163 * t16 * t138)) / 0.2e1 - (t105 * (-t77 * t137 * t11 + t71 * t139 * t123))) * e - t12 * t22 * t110 * t29, t17, 0, 0, 0, t20]; t2; t5; t5; t5; t1;];
    
    t1 = (n - p);
    t2 = (p ^ 2);
    t3 = (p * t2);
    t4 = (n ^ 2);
    t5 = (n * t4);
    t6 = (t3 * t1);
    t7 = 729 * n;
    t8 = -5832 * p;
    t9 = 19440 * t2;
    t10 = -18432 * t2;
    t11 = 4096 * p;
    t12 = (((((14823 + t8 + t7) * n + t9 + 121230) * n - 635580 * p + 508350) * n + 1249380 * t2 + 1146889) * n + t3 * (-1091360 + t10) + 1311659) * n + (t2 * (((-59904 + t11) * p + 357440) * p - 1110920) - 1664606) * p + 589312;
    t13 = 1890832 * p;
    t14 = 34560 * t6 * t4;
    t15 = p * ((((-98010 * n + 259200 * p) * n - 342720 * t2 - 1979580) * n + (226560 * t2 + 2568900) * p - 2945816) * n + t13) - t14;
    t16 = t12 + t15;
    t17 = (delta * m);
    t18 = n - t17;
    t19 = 2;
    t20 = t19 * p;
    t21 = -7 - t18 + t20;
    t22 = -l * t18;
    t23 = (-7 + t20) * l + t22;
    t24 = sin(t23);
    t23 = cos(t23);
    t25 = -79872 * t2;
    t26 = 20480 * p;
    t27 = (((((-9720 * p + t7 + 22437) * n + 45360 * t2 + 176850) * n - 103680 * t3 + 579510) * n + t2 * (126720 * t2 + 2472180) + 848909) * n + t3 * (-2433760 + t25) + 475997) * n + (t2 * (((-216576 + t26) * p + 884800) * p - 1751480) - 668282) * p + 35840;
    t28 = 1693760 * p;
    t29 = p * ((((-190350 * n + 613440 * p - 1094220) * n - 955200 * t2 - 2556660) * n + (725760 * t2 + 3691740) * p - 2417464) * n + t28);
    t30 = t27 + t29;
    t31 = -t18 - 5 + t20;
    t32 = (-5 + t20) * l + t22;
    t33 = sin(t32);
    t32 = cos(t32);
    t34 = 81 * n;
    t35 = 12288 * p;
    t36 = (((((216 * p + t34 - 1161) * n - 6480 * t2 - 21630) * n + 26880 * t3 - 42650) * n + t2 * (-47360 * t2 - 399820) - 7619) * n + t3 * (38912 * t2 + 433120) - 10765) * n + (t2 * (((81408 - t35) * p - 169920) * p + 95800) + 202) * p - 32256;
    t37 = 38304 * p;
    t38 = p * ((((25470 * n - 124800 * p + 156740) * n + 253760 * t2 + 170660) * n + (-234240 * t2 - 223980) * p - 19720) * n + t37);
    t39 = -3 - t18 + t20;
    t40 = (-3 + t20) * l + t22;
    t41 = sin(t40);
    t40 = cos(t40);
    t42 = 135 * n;
    t43 = 6144 * t2;
    t44 = (((((-1512 * p + 1035 + t42) * n + 5904 * t2 - 1534) * n - 8960 * t3 - 10126) * n + t2 * (2304 * t2 - 43908) + 12495) * n + t3 * (23776 + t43) + 731) * n + (t2 * (((13824 - t11) * p + 3520) * p - 23752) + 11714) * p + 17408;
    t45 = 16432 * p;
    t46 = p * ((((-5754 * n + 4032 * p + 20044) * n + 18240 * t2 + 12948) * n + (-31488 * t2 + 19860) * p - 6392) * n - t45);
    t47 = 1 + n - t17;
    t48 = -t47 + t20;
    t49 = t20 * l;
    t47 = -t47 * l + t49;
    t50 = sin(t47);
    t47 = cos(t47);
    t51 = (((((-1224 * p + t34 + 111) * n + 7152 * t2 - 1898) * n - 20224 * t3 + 1070) * n + t2 * (28416 * t2 - 48540) + 10329) * n + t3 * (37856 + t10) - 12445) * n + (t2 * (((13824 + t11) * p - 3520) * p - 23752) + 11714) * p - 17408;
    t45 = p * ((((198 * n - 8064 * p + 17636) * n + 30528 * t2 - 18588) * n + (-37632 * t2 + 51396) * p - 39256) * n + t45);
    t52 = t51 + t45;
    t17 = 1 - n + t17;
    t53 = t17 + t20;
    t17 = t17 * l + t49;
    t54 = sin(t17);
    t17 = cos(t17);
    t35 = (((((39 * n - 712 * p - 437) * n + 5200 * t2 + 1510) * n - 19200 * t3 + 170) * n + t2 * (37120 * t2 + 119980) - 10965) * n + t3 * (-34816 * t2 - 246560) + 10563) * n + (t2 * (((81408 + t35) * p + 169920) * p + 95800) + 202) * p + 32256;
    t37 = p * ((((7230 * n - 45120 * p - 23220) * n + 130880 * t2 + 10100) * n + (-172800 * t2 - 63420) * p + 56888) * n - t37);
    t55 = t35 + t37;
    t56 = 3 - t18 + t20;
    t57 = (3 + t20) * l + t22;
    t58 = sin(t57);
    t57 = cos(t57);
    t26 = (((((17 * n - 360 * p - 489) * n + 3120 * t2 + 5850) * n - 14080 * t3 - 36890) * n + t2 * (34560 * t2 + 479700) + 125205) * n + t3 * (-43008 * t2 - 1105440) - 192285) * n + (t2 * (((216576 + t26) * p + 884800) * p + 1751480) + 668282) * p + 35840;
    t28 = p * ((((8910 * n - 63360 * p - 88060) * n + 217920 * t2 + 427620) * n + (-357120 * t2 - 1562700) * p - 970056) * n + t28);
    t59 = t26 + t28;
    t60 = -t18 + 5 + t20;
    t61 = (5 + t20) * l + t22;
    t62 = sin(t61);
    t61 = cos(t61);
    t63 = -24;
    t64 = p * t63;
    t65 = 240 * t2;
    t66 = 1280 * t3;
    t67 = 3840 * t2;
    t68 = (((((n - 51 + t64) * n + t65 + 1110) * n - t66 - 13250) * n + t2 * (119940 + t67) + 91905) * n + t3 * (-338400 - t43) - 352947) * n + (t2 * (((59904 + t11) * p + 357440) * p + 1110920) + 1664606) * p + 589312;
    t13 = p * ((((1050 * n - 8640 * p - 18860) * n + 35520 * t2 + 174540) * n + (-72960 * t2 - 763860) * p - 835848) * n + t13);
    t69 = 7 - t18 + t20;
    t70 = (7 + t20) * l + t22;
    t71 = sin(t70);
    t70 = cos(t70);
    t31 = 1 / t31;
    t21 = 1 / t21;
    t60 = 1 / t60;
    t69 = 1 / t69;
    t72 = t31 ^ 2;
    t73 = t21 ^ 2;
    t74 = t60 ^ 2;
    t75 = l * t61;
    t76 = l * t32 * t31;
    t77 = t33 * t72;
    t78 = (t68 + t13) * t69;
    t39 = 1 / t39;
    t56 = 1 / t56;
    t79 = t56 ^ 2;
    t80 = l * t57;
    t81 = -(t80 * t56) + t58 * t79;
    t82 = l * t40;
    t83 = t39 * t41;
    t84 = (t36 + t38) * t39;
    t53 = 1 / t53;
    t48 = 1 / t48;
    t85 = t53 ^ 2;
    t86 = t54 * t85;
    t87 = t53 * t17;
    t88 = t87 * l;
    t89 = -t88 + t86;
    t90 = t47 * l;
    t91 = t48 * t50;
    t92 = -t91 + t90;
    t93 = (t44 + t46) * t48;
    t94 = 1 + p;
    t1 = p * t1;
    t95 = t1 * n;
    t96 = n * (n * (((18 * n - 192 * p + 61) * n + 704 * t2 - 234) * n - 1024 * t3 - 264) + 754) - 256 * p * (t5 + t3) + 512 * n * t3 * t94 + 1464 * t95 + 76 * t1 + 456;
    t97 = -t18 + t20;
    t49 = t22 + t49;
    t98 = sin(t49);
    t49 = cos(t49);
    t99 = 0.4e1 / 0.15e2 * p;
    t100 = 0.3e1 / 0.2e1 * n;
    t101 = ((((0.81e2 / 0.1280e4 * n + 0.477e3 / 0.512e3) * n + 0.675e3 / 0.128e3) * n + 0.685e3 / 0.48e2) * n + 0.35219e5 / 0.1920e4) * n + ((((-0.27e2 / 0.64e2 * n - 0.315e3 / 0.64e2) * n - 0.2647e4 / 0.128e3) * n - 0.3541e4 / 0.96e2) * n + (((0.9e1 / 0.8e1 * n + 0.39e2 / 0.4e1) * n + 0.5189e4 / 0.192e3) * n + ((-0.103e3 / 0.12e2 - t100) * n + (n + 0.17e2 / 0.6e1 - t99) * p - 0.565e3 / 0.48e2) * p + 0.9149e4 / 0.384e3) * p - 0.44939e5 / 0.1920e4) * p + 0.2843e4 / 0.320e3;
    t102 = 7680;
    t103 = t102 * t101;
    t104 = -t18 - 6 + t20;
    t105 = (-6 + t20) * l + t22;
    t106 = sin(t105);
    t105 = cos(t105);
    t107 = 12800;
    t108 = 162 * n;
    t107 = ((((4185 + t108) * n + 20910) * n + 35940) * n + 17792) * n + ((((-2160 * n - 29520) * n - 99080) * n - 108060) * n + (((8640 * n + 72960) * n + 151080) * n + ((-15360 * n - 76800) * n + (n * t107 - t11 + 29440) * p - 74880) * p + 78980) * p - 26892) * p - 2584;
    t109 = -4 - t18 + t20;
    t110 = (-4 + t20) * l + t22;
    t111 = sin(t110);
    t110 = cos(t110);
    t112 = 2048 * p;
    t113 = 3328 * n;
    t114 = 216 * n;
    t115 = 54 * n;
    t116 = ((((63 + t115) * n - 2032) * n - 276) * n - 288) * n + ((((2472 - t114) * n + 8884) * n - 3304) * n + (((-576 * n - 11520) * n - 12536) * n + ((16768 + t113) * n + (-4608 * n + t112 - 7936) * p + 5664) * p + 4852) * p - 1092) * p - 2840;
    t117 = 1 - p;
    t118 = -t117 * t19 - t18;
    t119 = -t19 * l * t117 + t22;
    t120 = sin(t119);
    t119 = cos(t119);
    t121 = 5632;
    t112 = ((((30 * n - 153) * n - 20) * n + 1272) * n - 1380) * n + ((((-424 * n + 2008) * n - 804) * n - 6400) * n + (((2240 * n - 8832) * n + 4456) * n + ((-5376 * n + 14976) * n + (n * t121 - t112 - 7936) * p - 5664) * p + 4852) * p + 1092) * p - 2840;
    t121 = t19 * t94 - t18;
    t122 = t19 * l * t94 + t22;
    t123 = sin(t122);
    t122 = cos(t122);
    t124 = 5120 * n;
    t125 = 1600 * n;
    t126 = 240 * n;
    t127 = ((((14 * n - 265) * n + 1970) * n - 6860) * n + 9100) * n + ((((3760 - t126) * n - 21560) * n + 49900) * n + (((-19200 + t125) * n + 73560) * n + ((40960 - t124) * n + (n * t102 - t11 - 29440) * p - 74880) * p - 78980) * p - 26892) * p + 2584;
    t128 = 4 - t18 + t20;
    t129 = (4 + t20) * l + t22;
    t130 = sin(t129);
    t129 = cos(t129);
    t131 = 0.1e1 / 0.192e3;
    t132 = 0.1e1 / 0.3840e4;
    t133 = 0.5e1 / 0.512e3;
    t99 = ((((-n * t132 + t133) * n - 0.29e2 / 0.192e3) * n + 0.155e3 / 0.128e3) * n - 0.81e2 / 0.16e2) * n + ((((t131 * n - 0.31e2 / 0.192e3) * n + 0.745e3 / 0.384e3) * n - 0.689e3 / 0.64e2) * n + (((-n / 0.24e2 + 0.1e1) * n - 0.1591e4 / 0.192e3) * n + ((n / 0.6e1 - 0.11e2 / 0.4e1) * n + (-n / 0.3e1 + t99 + 0.17e2 / 0.6e1) * p + 0.565e3 / 0.48e2) * p + 0.9149e4 / 0.384e3) * p + 0.44939e5 / 0.1920e4) * p + 0.2843e4 / 0.320e3;
    t18 = -t18 + 6 + t20;
    t22 = (6 + t20) * l + t22;
    t133 = sin(t22);
    t22 = cos(t22);
    t121 = 1 / t121;
    t118 = 1 / t118;
    t134 = t118 ^ 2;
    t135 = l * t119;
    t136 = (t135 * t118) - t120 * t134;
    t137 = l * t122;
    t138 = t121 * t123;
    t139 = -t138 + t137;
    t140 = t112 * t121;
    t128 = 1 / t128;
    t109 = 1 / t109;
    t141 = t109 ^ 2;
    t142 = l * t110;
    t143 = l * t129;
    t144 = t128 * t130;
    t145 = t127 * t128;
    t104 = 1 / t104;
    t18 = 1 / t18;
    t146 = t104 ^ 2;
    t102 = -(t102 * t99 * t18);
    t97 = 1 / t97;
    t147 = (l * t49) - t97 * t98;
    t148 = t96 * t97;
    t149 = 256 * p;
    t150 = 864 * n;
    t151 = -432 * n;
    t152 = (((810 + t34) * n + 2904) * n + 4394) * n + (((-3204 + t151) * n - 7564) * n + ((4224 + t150) * n + (-768 * n + t149 - 1856) * p + 4924) * p - 5644) * p + 2347;
    t153 = -1792 * n;
    t154 = 768 * p;
    t155 = 1440 * n;
    t156 = 27 * n;
    t151 = (((666 + t156) * n + 1684) * n + 686) * n + (((-3732 + t151) * n - 5572) * n + ((6336 + t155) * n + (-3392 + t153 + t154) * p + 4356) * p - 1220) * p - 495;
    t157 = 96 * p;
    t6 = -256 * t6;
    t158 = (((-144 * p + t156 + 66) * n + (84 + t157) * p - 214) * n + (-576 * p + 92) * p + 16) * n + ((448 * p + 284) * p - 232) * p - t6 - 307;
    t159 = 21 * n;
    t149 = -768 * n * t2 * t94 + (((-22 + t159) * n - 162) * n + 216) * n + (((-208 * n + 276) * n + 660) * n + ((448 + t149) * p + 672 * t4 - 284) * p - 232) * p + 307;
    t160 = 1280 * n;
    t161 = 144 * n;
    t154 = (((11 * n - 122) * n + 468) * n - 534) * n + (((1236 - t161) * n - 3140) * n + ((672 * n - 3840) * n + (3392 + t154 - t160) * p + 4356) * p + 1220) * p - 495;
    t162 = 16 * p;
    t163 = (((n - 26 - t162) * n + (324 + t157) * p + 264) * n + (-1344 * p - 2284) * p - 1250) * n + ((1856 * p + 4924) * p + 5644) * p + t6 + 2347;
    t164 = t39 ^ 2;
    t165 = t82 * t39;
    t166 = t41 * t164;
    t167 = t60 * t62;
    t168 = t163 * t60;
    t169 = t158 * t48;
    t170 = 6 * n;
    t171 = (5 + t170) * n;
    t172 = (-25 + t171) * n - 32 * t95 - 14;
    t173 = 128 * p;
    t114 = ((333 + t115) * n + 650) * n + ((-876 - t114) * n + (288 * n - t173 + 576) * p - 842) * p + 398;
    t115 = -78;
    t174 = 192 * n;
    t175 = 72 * n;
    t115 = n * (n * t115 - 29) + ((324 + t175) * n + (-288 + t173 - t174) * p + 74) * p + 68;
    t176 = 8 * n;
    t174 = ((-42 + t176) * n + 45) * n + ((252 - t175) * n + (-288 - t173 + t174) * p - 74) * p + 68;
    t175 = (0.4e1 / 0.3e1);
    t177 = 0.1e1 / 0.48e2;
    t178 = 0.11e2 / 0.32e2;
    t178 = ((-n * t177 + t178) * n - t19) * n + ((n / 0.4e1 - 0.23e2 / 0.8e1) * n + ((p * t175 - n + 6) * p) + 0.421e3 / 0.48e2) * p + 0.199e3 / 0.48e2;
    t179 = -0.96e2 * t178;
    t180 = t172 * t97;
    t181 = t109 * t111;
    t182 = t114 * t109;
    t183 = t174 * t121;
    t184 = 9 * n;
    t63 = n * (29 + t184) + (n * t63 + t162 - 38) * p + 22;
    t20 = n - t20;
    t185 = (t20 * p) + 0.5e1 / 0.4e1 * t94 - 0.3e1 / 0.8e1 * n * (1 - n);
    t186 = 0.8e1 * t185;
    t187 = (5 * n + t64 - 7) * n - 10 * t117 + 16 * t2;
    t188 = n * (n - 9) + (38 - t176 + t162) * p + 22;
    t189 = t48 ^ 2;
    t90 = t90 * t48;
    t190 = t50 * t189;
    t191 = t56 * t58;
    t192 = t188 * t56;
    t193 = t19 * n;
    t194 = 1 + t193;
    t195 = 8 * p;
    t170 = 7 - t195 + t170;
    t193 = -7 - t195 + t193;
    t195 = t194 * t97;
    t196 = t118 * t120;
    t197 = t170 * t118;
    t198 = 0.1e1 / 0.2e1;
    t199 = e ^ 2;
    t75 = ((-(m * (t16 * ((l * t23 * t21) - t24 * t73) + t30 * (t77 - t76) + t59 * (-(t75 * t60) + t62 * t74) + t78 * ((l * t70) - t69 * t71)) / 0.46080e5 - m * (t55 * t81 + t84 * (-t83 + t82)) / 0.15360e5 + m * (t52 * t89 + t93 * t92) / 0.9216e4) * e - t132 * m * (t107 * (-(t142 * t109) + t111 * t141) + t145 * (-t144 + t143)) + m * (t116 * t136 + t140 * t139) / 0.1536e4 + m * (t103 * (-(l * t105 * t104) + t106 * t146) + t102 * ((l * t22) - t133 * t18)) / 0.7680e4 - t148 * m * t147 / 0.384e3) * e + t131 * m * (t149 * t89 + t169 * t92) - m * (t151 * (t166 - t165) + t152 * (-t77 + t76) + t154 * t81 + t168 * (-t167 + t75)) / 0.384e3);
    t76 = -(0.1e1 / 0.2560e4);
    t77 = 0.5e1 / 0.1536e4;
    t81 = 0.1e1 / 0.768e3;
    t82 = (0.1e1 / 0.32e2);
    t16 = ((-t187 * t53 * t54 + t91 * t186 + t192 * t58 + t83 * t63) * e / 0.4e1 - t138 * t193 / 0.4e1 + t197 * t120 / 0.4e1 + ((((t16 * t21 * t24 / 0.7680e4 - t30 * t31 * t33 / 0.7680e4 - t167 * t59 / 0.7680e4 + t78 * t71 / 0.7680e4 + t76 * (-t191 * t55 + t84 * t41) - t52 * t53 * t54 / 0.1536e4 + t93 * t50 / 0.1536e4) * e - t81 * (t181 * t107 - t145 * t130) - t77 * (t196 * t116 + t140 * t123) + t103 * t104 * t106 / 0.1536e4 - t102 * t133 / 0.1536e4 + 0.5e1 / 0.384e3 * t148 * t98) * e - t177 * (-t149 * t53 * t54 + t169 * t50) + t152 * t31 * t33 / 0.96e2 - t83 * t151 / 0.96e2 - t191 * t154 / 0.96e2 + t168 * t62 / 0.96e2) * e + t82 * (t182 * t111 - t144 * t179) + t196 * t115 / 0.16e2 + t183 * t123 / 0.16e2 - 0.3e1 / 0.16e2 * t180 * t98) * t199 + t195 * t198 * t98);
    t12 = ((((((t17 * (-t45 / 0.9216e4 - t51 / 0.9216e4) + t23 * (t15 / 0.46080e5 + t12 / 0.46080e5) + t32 * (-t29 / 0.46080e5 - t27 / 0.46080e5) + t40 * (-t38 / 0.15360e5 - t36 / 0.15360e5) + t47 * (t46 / 0.9216e4 + t44 / 0.9216e4) + t57 * (t37 / 0.15360e5 + t35 / 0.15360e5) + t61 * (-t28 / 0.46080e5 - t26 / 0.46080e5) + t70 * (t13 / 0.46080e5 + t68 / 0.46080e5)) * e + t105 * t101 - t110 * t132 * t107 - (t119 * t116) / 0.1536e4 - (t122 * t112) / 0.1536e4 + t129 * t132 * t127 + t22 * t99 + (t49 * t96) / 0.384e3) * e + t17 * t131 * t149 + (t32 * t152) / 0.384e3 - (t40 * t151) / 0.384e3 - t47 * t131 * t158 - (t57 * t154) / 0.384e3 + (t61 * t163) / 0.384e3) * e + (t110 * t114) / 0.96e2 + t119 * t177 * t115 + t122 * t177 * t174 + t129 * t178 - (t49 * t172) / 0.16e2) * e - (t17 * t187) / 0.8e1 + (t40 * t63) / 0.8e1 + t47 * t185 + (t57 * t188) / 0.8e1) * e + t49 * t198 * t194 + (t119 * t170) / 0.4e1 - (t122 * t193) / 0.4e1);
    t7 = (((((15552 + t8 + t7) * n + t9 + 133785) * n - 702000 * p + 591780) * n + 1381140 * t2 + 1413079) * n + t3 * (-1207520 + t10) + 1717148) * n + (t2 * (((-62976 + t11) * p + 395840) * p - 1298840) - 2185610) * p + 823543;
    t8 = 2339572 * p;
    t9 = p * ((((-102870 * n + 272160 * p) * n - 360000 * t2 - 2307720) * n + (238080 * t2 + 2999040) * p - 3637136) * n + t8) - t14;
    t10 = t7 + t9;
    t13 = 28672 * p;
    t14 = (((((2187 * n - 21384 * p + 44550) * n + 84240 * t2 + 333855) * n - 172800 * t3 + 1211280) * n + t2 * (195840 * t2 + 4194060) + 2273077) * n + t3 * (-116736 * t2 - 3983840) + 2110210) * n + (t2 * (((-314880 + t13) * p + 1407680) * p - 3275800) - 2793610) * p + 765625;
    t15 = 4186444 * p;
    t26 = p * ((((-339390 * n + 1015200 * p - 1943640) * n - 1496640 * t2 - 5102160) * n + (1090560 * t2 + 7105800) * p - 6188816) * n + t15);
    t27 = t14 + t26;
    t28 = (((((891 * n - 10584 * p + 15336) * n + 49680 * t2 + 88635) * n - 119040 * t3 + 222920) * n + t2 * (154880 * t2 + 1315820) + 266241) * n + t3 * (-104448 * t2 - 1334560) + 151632) * n + (t2 * (((-188928 + t13) * p + 498880) * p - 678120) - 203670) * p + 35721;
    t29 = 507084 * p;
    t30 = p * ((((-135090 * n + 457440 * p - 564800) * n - 751040 * t2 - 985080) * n + (601600 * t2 + 1426240) * p - 739320) * n + t29);
    t35 = t28 + t30;
    t36 = (((((351 * n - 4968 * p + 3798) * n + 27792 * t2 + 11491) * n - 78592 * t3 + 8932) * n + t2 * (119040 * t2 + 155436) - 99) * n + t3 * (-92160 * t2 - 140512) - 218) * n + (t2 * (((-62976 + t13) * p + 44480) * p - 11480) - 602) * p + 49;
    t37 = 1708 * p;
    t38 = p * ((((-37590 * n + 138912 * p - 71560) * n - 241728 * t2 - 30288) * n + (199680 * t2 + 32520) * p - 88) * n + t37);
    t44 = t36 + t38;
    t25 = (((((-2232 * p + t42 + 96) * n + 14736 * t2 - 665) * n - 49408 * t3 - 316) * n + t2 * (88320 * t2 + 780) + 1521) * n + t3 * (-37408 + t25) - 820) * n + (t2 * (((62976 + t13) * p + 44480) * p + 11480) + 602) * p + 49;
    t45 = 1110 * n;
    t37 = p * ((((-17952 * p + 4304 + t45) * n + 72768 * t2 - 312) * n + (-115200 * t2 - 1920) * p - 3328) * n + t37);
    t46 = t25 + t37;
    t51 = (((((51 * n - 968 * p - 682) * n + 7440 * t2 + 3975) * n - 29440 * t3 - 14040) * n + t2 * (62720 * t2 + 305420) + 34005) * n + t3 * (-67584 * t2 - 660960) - 52038) * n + (t2 * (((188928 + t13) * p + 498880) * p + 678120) + 203670) * p + 35721;
    t29 = p * ((((11570 * n - 75360 * p - 58680) * n + 233920 * t2 + 166960) * n + (-343040 * t2 - 608120) * p - 274848) * n + t29);
    t13 = (((((19 * n - 408 * p - 600) * n + 3600 * t2 + 8115) * n - 16640 * t3 - 60880) * n + t2 * (42240 * t2 + 688620) + 270705) * n + t3 * (-55296 * t2 - 1646880) - 683400) * n + (t2 * (((314880 + t13) * p + 1407680) * p + 3275800) + 2793610) * p + 765625;
    t15 = p * ((((11070 * n - 80160 * p - 123680) * n + 283200 * t2 + 717960) * n + (-483840 * t2 - 2721600) * p - 2184072) * n + t15);
    t52 = t13 + t15;
    t11 = (((((n - 54 + t64) * n + t65 + 1245) * n - t66 - 15740) * n + t2 * (133620 + t67) + 115515) * n + t3 * (-375840 - t43) - 468462) * n + (t2 * (((62976 + t11) * p + 395840) * p + 1298840) + 2185610) * p + 823543;
    t8 = p * ((((-9120 * p + t45 - 21080) * n + 37440 * t2 + 206160) * n + (-76800 * t2 - 897480) * p - 1042008) * n + t8);
    t43 = (l * t33 * t31);
    t45 = t32 * t72;
    t55 = (t11 + t8) * t69;
    t59 = (t83 * l);
    t64 = t40 * t164;
    t65 = t64 + t59;
    t66 = t56 * t57;
    t67 = (l * t58 + t66);
    t68 = (t51 + t29) * t56;
    t72 = (t91 * l);
    t76 = t189 * t47;
    t78 = t76 + t72;
    t83 = l * t54;
    t84 = (t83 + t87);
    t91 = -(1024 * p);
    t92 = (4320 * n);
    t93 = ((((((243 * n + 3780) * n + 22680) * n + 65300) * n + 89688) * n) + (((((-1620 * n - 19980) * n - 89040) * n - 169100) * n) + (((0.39600e5 + t92) * n + 0.116500e6) * n + (((-5760 * n - 34880) * n) + ((3840 * n) + 0.11520e5 + t91) * p - 0.50800e5) * p + 0.109440e6) * p - 0.114768e6) * p + 0.46656e5);
    t96 = -1536 * p;
    t92 = (((((2295 + t108) * n + 11280) * n + 24830) * n + 24872) * n) + (((((-1350 * n - 14130) * n - 49370) * n - 69460) * n) + (((0.31920e5 + t92) * n + 0.71060e5) * n + (((-6720 * n - 31520) * n + (11520 + t96 + t124) * p - 33720) * p) + 0.48120e5) * p - 0.33552e5) * p + 0.9216e4;
    t99 = -(3072 * p);
    t42 = ((((((1440 + t42) * n + 4652) * n + 5744) * n + 2888) * n) + (((((-1404 * n - 10308) * n - 21880) * n - 16292) * n) + ((((5472 * n + 26352) * n + 33284) * n) + (((-10112 * n - 28864) * n) + ((8960 * n) + 0.11520e5 + t99) * p - 0.16464e5) * p + 0.11328e5) * p - 0.3888e4) * p + 0.576e3);
    t99 = (((((-140 + t159) * n + 408) * n - 780) * n + 1000) * n) + (((((-316 * n + 1884) * n - 4704) * n + 6364) * n) + ((((1824 * n - 8880) * n + 16108) * n) + (((-4992 * n + 17216) * n) + ((6400 * n) + t99 - 0.11520e5) * p - 0.16464e5) * p - 0.11328e5) * p - 0.3888e4) * p - 0.576e3;
    t96 = ((((4 * n - 85) * n + 750) * n - 3490) * n + 8680) * n + ((((-70 * n + 1230) * n - 8410) * n + 26780) * n + (((480 * n - 6480) * n + 30100) * n + ((14560 - t125) * n + (2560 * n + t96 - 11520) * p - 33720) * p - 48120) * p - 33552) * p - 9216;
    t5 = ((((t4 * (-40 + n) - 5640) * n + 25080) * n) + ((((-20 * t4 - 8440) * n + 49780) * n) + ((((160 * n - 4080) * n + 35900) * n) + (((-640 * n + 11200) * n) + (-0.11520e5 + t160 + t91) * p - 0.50800e5) * p - 0.109440e6) * p - 0.114768e6) * p + (660 * t5 * t94) - 0.46656e5);
    t91 = t181 * l + (t110 * t141);
    t101 = l * t130 + (t128 * t129);
    t102 = t96 * t128;
    t103 = (t5 * t18);
    t107 = (l * t98 + (t97 * t49));
    t3 = t20 * ((((35 + t184 - t157) * n + 352 * t2 + 13) * n - 512 * t3 - 22) * n + 256 * t2 ^ 2 + 44 * t1 - 176 * t95);
    t95 = t3 * t97;
    t108 = (t196 * l + (t119 * t134));
    t112 = (l * t123 + (t121 * t122));
    t114 = (t99 * t121);
    t116 = (0.9e1 / 0.16e2 * n);
    t2 = (((-0.27e2 / 0.256e3 * t4 - 0.553e3 / 0.128e3) * n - 0.1351e4 / 0.192e3) * n + (((0.285e3 / 0.64e2 + t116) * n + 0.541e3 / 0.48e2) * n + (-0.47e2 / 0.8e1 * n + (n - p / 0.3e1 + 0.31e2 / 0.12e2) * p - 0.1411e4 / 0.192e3) * p + 0.145e3 / 0.16e2) * p - 0.9e1 / 0.8e1 * t4 * (n + t2) - 0.3125e4 / 0.768e3);
    t4 = -768 * t2;
    t124 = 1280 * p;
    t113 = (((189 * n + 1692) * n + 4822) * n + 5376) * n + (((-1296 * n - 7908) * n - 14088) * n + ((3168 * n + 12000) * n + (-5952 + t124 - t113) * p + 10124) * p - 7464) * p + 2025;
    t125 = -2816;
    t34 = (((420 + t34) * n + 536) * n + 198) * n + (((-720 * n - 2244) * n - 1552) * n + ((2208 * n + 3744) * n + (n * t125 + t124 - 1984) * p + 1076) * p - 252) * p + 25;
    t125 = (((33 * n - 64) * n + 60) * n - 54) * n + (((-368 * n + 708) * n - 600) * n + ((-2208 + t155) * n + (-2304 * n + 1984 + t124) * p + 1076) * p + 252) * p + 25;
    t124 = (((13 * n - 168) * n + 858) * n - 2088) * n + (((-176 * n + 1764) * n - 6160) * n + ((-5856 + t150) * n + (5952 + t153 + t124) * p + 10124) * p + 7464) * p + 2025;
    t6 = (((n - 28 - t162) * n + (348 + t157) * p + 306) * n + (-1440 * p - 2632) * p - 1556) * n + ((1984 * p + 5644) * p + 6960) * p + t6 + 3125;
    t127 = t60 * t61;
    t60 = t6 * t60;
    t132 = (64 * p);
    t134 = ((((180 + t156) * n + 382) * n) + (((-108 * n - 474) * n) + (0.312e3 - t132 + t161) * p - 0.496e3) * p + 0.256e3);
    t126 = ((129 + t156) * n + 170) * n + ((-408 - t161) * n + (312 - t173 + t126) * p - 248) * p + 64;
    t138 = ((0.5e1 / 0.48e2 * n - 0.11e2 / 0.16e2) * n + 0.13e2 / 0.8e1) * n + ((0.9e1 / 0.2e1 - n) * n + (-0.8e1 / 0.3e1 * p + (3 * n) - 0.13e2 / 0.2e1) * p - 0.31e2 / 0.6e1) * p - t175;
    t140 = 48;
    t132 = (((-18 + n) * n + 114) * n) + (((-12 * n + 150) * n) + ((n * t140) - t132 - 0.312e3) * p - 0.496e3) * p - 0.256e3;
    t1 = t20 * (-32 * t1 + t171);
    t141 = t1 * t97;
    t144 = t132 * t128;
    t145 = (t140 * t138 * t121);
    t148 = -(0.27e2 / 0.16e2);
    t19 = (-(n * (t19 + t116)) + (-p + 0.21e2 / 0.8e1 + t100) * p + t148);
    t100 = -(16 * t19);
    t116 = n * (15 * n + 26) + (p * t140 - 56 * n - 42) * p + 9;
    t140 = (n * (-0.7e1 / 0.16e2 * n + 0.1e1) + (0.5e1 / 0.2e1 * n - (3 * p) - 0.21e2 / 0.8e1) * p - 0.9e1 / 0.16e2);
    t148 = -16 * t140;
    t149 = n * (n - 10) + (42 - t176 + t162) * p + 27;
    t85 = t85 * t17;
    t83 = t83 * t53;
    t56 = t149 * t56;
    t150 = -(0.3e1 / 0.4e1);
    t117 = n * t150 - t117;
    t150 = -4 * t117;
    t94 = n - 4 * t94;
    t97 = t20 * t97;
    t151 = t94 * t121;
    t18 = ((((-m * (t10 * (l * t24 * t21 + (t23 * t73)) + (t27 * (-t45 - t43)) + t52 * (t167 * l + (t61 * t74)) - t55 * (l * t71 + (t69 * t70))) / 0.92160e5 - m * (t35 * t65 - t68 * t67) / 0.30720e5 + m * (-t46 * t53 * t84 + t44 * t78) / 0.18432e5) * e + t131 * t95 * m * t107 + m * (t102 * t101 + t91 * t92) / 0.1920e4 - m * (t93 * (l * t106 * t104 + (t105 * t146)) + t103 * (l * t133 + (t18 * t22))) / 0.7680e4 - m * (t108 * t42 + t114 * t112) / 0.1536e4) * e - t81 * m * ((t113 * (-t64 - t59)) + t124 * (t191 * l + (t57 * t79)) + (t4 * (t45 + t43)) - t60 * (l * t62 + t127)) - m * (-t125 * t53 * t84 + t34 * t78) / 0.384e3) * e + t177 * m * (t108 * t126 + t145 * t112) - m * (t144 * t101 + t134 * t91) / 0.96e2 - t141 * m * t107 / 0.16e2);
    t43 = -(0.1e1 / 0.3072e4);
    t45 = (0.1e1 / 0.5120e4);
    t4 = (-t199 * (((-((t45 * (t35 * t39 * t40 - t68 * t57)) + (t43 * (t44 * t48 * t47 - t87 * t46)) + (t10 * t21 * t23) / 0.15360e5 - (t27 * t31 * t32) / 0.15360e5 + (t127 * t52) / 0.15360e5 - (t55 * t70) / 0.15360e5) * e - (t93 * t104 * t105) / 0.1536e4 - t103 * t22 / 0.1536e4 - t77 * (t42 * t118 * t119 + t114 * t122) + t92 * t109 * t110 / 0.384e3 + (t102 * t129) / 0.384e3 + 0.5e1 / 0.192e3 * t95 * t49) * e - t131 * (-t113 * t39 * t40 + t4 * t31 * t32 + t66 * t124 - t60 * t61) - (t34 * t48 * t47) / 0.96e2 + (t87 * t125) / 0.96e2) * e - t82 * ((t134 * t109 * t110) + t144 * t129) + (t126 * t118 * t119) / 0.16e2 + (t145 * t122) / 0.16e2 - 0.3e1 / 0.16e2 * t141 * t49) - (t97 * t49) + (t150 * t118 * t119) / 0.4e1 + (t151 * t122) / 0.4e1 + (t100 * t39 * t40 / 0.8e1 - (t116 * t48 * t47) / 0.8e1 + (t87 * t148) / 0.8e1 - (t56 * t57) / 0.8e1) * e);
    t10 = (t50 - t54);
    t1 = (-t198 * t10 + ((((((t24 * (-t9 / 0.92160e5 - t7 / 0.92160e5) + t33 * (t26 / 0.92160e5 + t14 / 0.92160e5) + t41 * (-t30 / 0.30720e5 - t28 / 0.30720e5) + t50 * (t38 / 0.18432e5 + t36 / 0.18432e5) + t54 * (-t37 / 0.18432e5 - t25 / 0.18432e5) + t58 * (t29 / 0.30720e5 + t51 / 0.30720e5) + t62 * (-t15 / 0.92160e5 - t13 / 0.92160e5) + t71 * (t8 / 0.92160e5 + t11 / 0.92160e5)) * e - t106 * t93 / 0.7680e4 + t111 * t92 / 0.1920e4 - t120 * t42 / 0.1536e4 - t123 * t99 / 0.1536e4 + t130 * t96 / 0.1920e4 - t133 * t5 / 0.7680e4 + t3 * t131 * t98) * e + t41 * t81 * t113 + t2 * t33 - t50 * t34 / 0.384e3 + t54 * t125 / 0.384e3 - t58 * t81 * t124 + t62 * t81 * t6) * e - t111 * t134 / 0.96e2 + t120 * t177 * t126 + t123 * t138 - t130 * t132 / 0.96e2 - t1 * t98 / 0.16e2) * e + t19 * t41 + t50 * t116 / 0.16e2 + t54 * t140 + t58 * t149 / 0.16e2) * e + t117 * t120 - t123 * t94 / 0.4e1 + t20 * t98) * e);
    I1lgx = [-((t75 * e + t180 * m * t147 / 0.16e2 - t177 * m * (t115 * t136 + t183 * t139) - m * (t179 * (t130 * (t128 ^ 2) - (t143 * t128)) + t182 * (-t181 + t142)) / 0.96e2) * t199 - t198 * t195 * m * t147 - m * (t193 * (t123 * (t121 ^ 2) - (t137 * t121)) + t197 * (-t196 + t135)) / 0.4e1) * e - m * (t190 - t90 - t88 + t86) + m * (t186 * (-t190 + t90) + t187 * t89 + t63 * (-t166 + t165) + t192 * (-t191 + t80)) * t199 / 0.8e1 t16 0 0 0 e * t12 + t17 + t47; -t198 * m * (t76 + t72 - t85 - t83) + (t199 * t18 + t97 * m * t107 - m * (t108 * t150 + t151 * t112) / 0.4e1) * e - m * (t100 * t65 + (t116 * (-t76 - t72)) + t148 * (t83 + t85) - (t56 * t67)) * t199 / 0.16e2 t4 0 0 0 t1;];
    
    t1 = (n - p);
    t2 = (p ^ 2);
    t3 = (p * t2);
    t4 = (n ^ 2);
    t5 = (n * t4);
    t6 = (t3 * t1);
    t7 = 729 * n;
    t8 = -5832 * p;
    t9 = 19440 * t2;
    t10 = -18432 * t2;
    t11 = 4096 * p;
    t12 = (((((14823 + t8 + t7) * n + t9 + 121230) * n - 635580 * p + 508350) * n + 1249380 * t2 + 1146889) * n + t3 * (-1091360 + t10) + 1311659) * n + (t2 * (((t11 - 59904) * p + 357440) * p - 1110920) - 1664606) * p + 589312;
    t13 = 1890832 * p;
    t14 = 34560 * t6 * t4;
    t15 = p * ((((-98010 * n + 259200 * p) * n - 342720 * t2 - 1979580) * n + (226560 * t2 + 2568900) * p - 2945816) * n + t13) - t14;
    t16 = t12 + t15;
    t17 = (delta * m);
    t18 = n - t17;
    t19 = 2;
    t20 = t19 * p;
    t21 = -t18 - 7 + t20;
    t22 = -l * t18;
    t23 = (-7 + t20) * l + t22;
    t24 = cos(t23);
    t23 = sin(t23);
    t25 = -79872 * t2;
    t26 = 20480 * p;
    t27 = (((((-9720 * p + t7 + 22437) * n + 45360 * t2 + 176850) * n - 103680 * t3 + 579510) * n + t2 * (126720 * t2 + 2472180) + 848909) * n + t3 * (-2433760 + t25) + 475997) * n + (t2 * (((-216576 + t26) * p + 884800) * p - 1751480) - 668282) * p + 35840;
    t28 = 1693760 * p;
    t29 = p * ((((-190350 * n + 613440 * p - 1094220) * n - 955200 * t2 - 2556660) * n + (725760 * t2 + 3691740) * p - 2417464) * n + t28);
    t30 = t27 + t29;
    t31 = -t18 - 5 + t20;
    t32 = (-5 + t20) * l + t22;
    t33 = cos(t32);
    t32 = sin(t32);
    t34 = 81 * n;
    t35 = 12288 * p;
    t36 = (((((216 * p + t34 - 1161) * n - 6480 * t2 - 21630) * n + 26880 * t3 - 42650) * n + t2 * (-47360 * t2 - 399820) - 7619) * n + t3 * (38912 * t2 + 433120) - 10765) * n + (t2 * (((-t35 + 81408) * p - 169920) * p + 95800) + 202) * p - 32256;
    t37 = 38304 * p;
    t38 = p * ((((25470 * n - 124800 * p + 156740) * n + 253760 * t2 + 170660) * n + (-234240 * t2 - 223980) * p - 19720) * n + t37);
    t39 = -t18 - 3 + t20;
    t40 = (-3 + t20) * l + t22;
    t41 = cos(t40);
    t40 = sin(t40);
    t42 = 135 * n;
    t43 = 6144 * t2;
    t44 = (((((-1512 * p + 1035 + t42) * n + 5904 * t2 - 1534) * n - 8960 * t3 - 10126) * n + t2 * (2304 * t2 - 43908) + 12495) * n + t3 * (23776 + t43) + 731) * n + (t2 * (((-t11 + 13824) * p + 3520) * p - 23752) + 11714) * p + 17408;
    t45 = 16432 * p;
    t46 = p * ((((-5754 * n + 4032 * p + 20044) * n + 18240 * t2 + 12948) * n + (-31488 * t2 + 19860) * p - 6392) * n - t45);
    t47 = t44 + t46;
    t48 = 1 + n - t17;
    t49 = -t48 + t20;
    t50 = t20 * l;
    t48 = -t48 * l + t50;
    t51 = cos(t48);
    t48 = sin(t48);
    t52 = (((((-1224 * p + t34 + 111) * n + 7152 * t2 - 1898) * n - 20224 * t3 + 1070) * n + t2 * (28416 * t2 - 48540) + 10329) * n + t3 * (37856 + t10) - 12445) * n + (t2 * (((t11 + 13824) * p - 3520) * p - 23752) + 11714) * p - 17408;
    t45 = p * ((((198 * n - 8064 * p + 17636) * n + 30528 * t2 - 18588) * n + (-37632 * t2 + 51396) * p - 39256) * n + t45);
    t17 = 1 - n + t17;
    t53 = t17 + t20;
    t17 = t17 * l + t50;
    t54 = cos(t17);
    t17 = sin(t17);
    t35 = (((((39 * n - 712 * p - 437) * n + 5200 * t2 + 1510) * n - 19200 * t3 + 170) * n + t2 * (37120 * t2 + 119980) - 10965) * n + t3 * (-34816 * t2 - 246560) + 10563) * n + (t2 * (((t35 + 81408) * p + 169920) * p + 95800) + 202) * p + 32256;
    t37 = p * ((((7230 * n - 45120 * p - 23220) * n + 130880 * t2 + 10100) * n + (-172800 * t2 - 63420) * p + 56888) * n - t37);
    t55 = t35 + t37;
    t56 = -t18 + 3 + t20;
    t57 = (3 + t20) * l + t22;
    t58 = cos(t57);
    t57 = sin(t57);
    t26 = (((((17 * n - 360 * p - 489) * n + 3120 * t2 + 5850) * n - 14080 * t3 - 36890) * n + t2 * (34560 * t2 + 479700) + 125205) * n + t3 * (-43008 * t2 - 1105440) - 192285) * n + (t2 * (((216576 + t26) * p + 884800) * p + 1751480) + 668282) * p + 35840;
    t28 = p * ((((8910 * n - 63360 * p - 88060) * n + 217920 * t2 + 427620) * n + (-357120 * t2 - 1562700) * p - 970056) * n + t28);
    t59 = t26 + t28;
    t60 = -t18 + 5 + t20;
    t61 = (5 + t20) * l + t22;
    t62 = cos(t61);
    t61 = sin(t61);
    t63 = -24;
    t64 = p * t63;
    t65 = 240 * t2;
    t66 = 1280 * t3;
    t67 = 3840 * t2;
    t68 = (((((n + t64 - 51) * n + t65 + 1110) * n - t66 - 13250) * n + t2 * (t67 + 119940) + 91905) * n + t3 * (-338400 - t43) - 352947) * n + (t2 * (((t11 + 59904) * p + 357440) * p + 1110920) + 1664606) * p + 589312;
    t13 = p * ((((1050 * n - 8640 * p - 18860) * n + 35520 * t2 + 174540) * n + (-72960 * t2 - 763860) * p - 835848) * n + t13);
    t69 = -t18 + 7 + t20;
    t70 = (7 + t20) * l + t22;
    t71 = cos(t70);
    t70 = sin(t70);
    t56 = 1 / t56;
    t39 = 1 / t39;
    t72 = t56 ^ 2;
    t73 = l * t57;
    t74 = l * t40;
    t75 = t39 * t41;
    t76 = (t36 + t38) * t39;
    t69 = 1 / t69;
    t21 = 1 / t21;
    t60 = 1 / t60;
    t31 = 1 / t31;
    t77 = t31 ^ 2;
    t78 = (l * t32 * t31) + t33 * t77;
    t79 = t60 ^ 2;
    t80 = (l * t61 * t60) + t62 * t79;
    t81 = t21 ^ 2;
    t82 = (t68 + t13) * t69;
    t49 = 1 / t49;
    t53 = 1 / t53;
    t83 = t49 ^ 2;
    t84 = t49 * t48;
    t85 = t84 * l;
    t86 = t51 * t83;
    t87 = t86 + t85;
    t88 = t17 * l;
    t89 = t53 * t54;
    t90 = t89 + t88;
    t91 = (t52 + t45) * t53;
    t92 = 1 + p;
    t1 = p * t1;
    t93 = t1 * n;
    t94 = n * (n * (((18 * n - 192 * p + 61) * n + 704 * t2 - 234) * n - 1024 * t3 - 264) + 754) - 256 * p * (t5 + t3) + 512 * n * t3 * t92 + 1464 * t93 + 76 * t1 + 456;
    t95 = -t18 + t20;
    t50 = t22 + t50;
    t96 = cos(t50);
    t50 = sin(t50);
    t97 = -0.17e2 / 0.6e1;
    t98 = 0.4e1 / 0.15e2 * p;
    t99 = 0.3e1 / 0.2e1 * n;
    t100 = ((((-0.81e2 / 0.1280e4 * n - 0.477e3 / 0.512e3) * n - 0.675e3 / 0.128e3) * n - 0.685e3 / 0.48e2) * n - 0.35219e5 / 0.1920e4) * n + ((((0.27e2 / 0.64e2 * n + 0.315e3 / 0.64e2) * n + 0.2647e4 / 0.128e3) * n + 0.3541e4 / 0.96e2) * n + (((-0.9e1 / 0.8e1 * n - 0.39e2 / 0.4e1) * n - 0.5189e4 / 0.192e3) * n + ((t99 + 0.103e3 / 0.12e2) * n + (-n + t98 + t97) * p + 0.565e3 / 0.48e2) * p - 0.9149e4 / 0.384e3) * p + 0.44939e5 / 0.1920e4) * p - 0.2843e4 / 0.320e3;
    t101 = 7680;
    t102 = -t18 - 6 + t20;
    t103 = (-6 + t20) * l + t22;
    t104 = cos(t103);
    t103 = sin(t103);
    t105 = 12800;
    t106 = 162 * n;
    t105 = ((((t106 + 4185) * n + 20910) * n + 35940) * n + 17792) * n + ((((-2160 * n - 29520) * n - 99080) * n - 108060) * n + (((8640 * n + 72960) * n + 151080) * n + ((-15360 * n - 76800) * n + (n * t105 - t11 + 29440) * p - 74880) * p + 78980) * p - 26892) * p - 2584;
    t107 = -4 - t18 + t20;
    t108 = (-4 + t20) * l + t22;
    t109 = cos(t108);
    t108 = sin(t108);
    t110 = 2048 * p;
    t111 = 3328 * n;
    t112 = 216 * n;
    t113 = 54 * n;
    t114 = ((((t113 + 63) * n - 2032) * n - 276) * n - 288) * n + ((((-t112 + 2472) * n + 8884) * n - 3304) * n + (((-576 * n - 11520) * n - 12536) * n + ((16768 + t111) * n + (-4608 * n + t110 - 7936) * p + 5664) * p + 4852) * p - 1092) * p - 2840;
    t115 = 1 - p;
    t116 = -t115 * t19 - t18;
    t117 = -t19 * l * t115 + t22;
    t118 = cos(t117);
    t117 = sin(t117);
    t119 = 5632;
    t110 = ((((30 * n - 153) * n - 20) * n + 1272) * n - 1380) * n + ((((-424 * n + 2008) * n - 804) * n - 6400) * n + (((2240 * n - 8832) * n + 4456) * n + ((-5376 * n + 14976) * n + (n * t119 - t110 - 7936) * p - 5664) * p + 4852) * p + 1092) * p - 2840;
    t119 = t19 * t92 - t18;
    t120 = t19 * l * t92 + t22;
    t121 = cos(t120);
    t120 = sin(t120);
    t122 = 5120 * n;
    t123 = 1600 * n;
    t124 = 240 * n;
    t125 = ((((14 * n - 265) * n + 1970) * n - 6860) * n + 9100) * n + ((((-t124 + 3760) * n - 21560) * n + 49900) * n + (((-19200 + t123) * n + 73560) * n + ((40960 - t122) * n + (n * t101 - t11 - 29440) * p - 74880) * p - 78980) * p - 26892) * p + 2584;
    t126 = 4 - t18 + t20;
    t127 = (4 + t20) * l + t22;
    t128 = cos(t127);
    t127 = sin(t127);
    t129 = 0.1e1 / 0.3e1;
    t130 = 0.1e1 / 0.192e3;
    t131 = 0.31e2 / 0.192e3;
    t97 = ((((n / 0.3840e4 - 0.5e1 / 0.512e3) * n + 0.29e2 / 0.192e3) * n - 0.155e3 / 0.128e3) * n + 0.81e2 / 0.16e2) * n + ((((-n * t130 + t131) * n - 0.745e3 / 0.384e3) * n + 0.689e3 / 0.64e2) * n + (((n / 0.24e2 - 0.1e1) * n + 0.1591e4 / 0.192e3) * n + ((-n / 0.6e1 + 0.11e2 / 0.4e1) * n + (n * t129 + t97 - t98) * p - 0.565e3 / 0.48e2) * p - 0.9149e4 / 0.384e3) * p - 0.44939e5 / 0.1920e4) * p - 0.2843e4 / 0.320e3;
    t98 = t101 * t97;
    t18 = -t18 + 6 + t20;
    t22 = (6 + t20) * l + t22;
    t131 = cos(t22);
    t22 = sin(t22);
    t95 = 1 / t95;
    t132 = (l * t50) + t95 * t96;
    t133 = t94 * t95;
    t107 = 1 / t107;
    t126 = 1 / t126;
    t134 = l * t127;
    t135 = l * t108;
    t136 = t107 * t109;
    t137 = t105 * t107;
    t102 = 1 / t102;
    t18 = 1 / t18;
    t138 = t18 * t22;
    t139 = t102 * t104;
    t101 = -(t101 * t100 * t102);
    t119 = 1 / t119;
    t116 = 1 / t116;
    t140 = t116 ^ 2;
    t141 = (l * t117 * t116) + t118 * t140;
    t142 = (l * t120) + t119 * t121;
    t143 = t110 * t119;
    t144 = 256 * p;
    t145 = 864 * n;
    t146 = -432 * n;
    t147 = (((t34 + 810) * n + 2904) * n + 4394) * n + (((t146 - 3204) * n - 7564) * n + ((t145 + 4224) * n + (-768 * n + t144 - 1856) * p + 4924) * p - 5644) * p + 2347;
    t148 = -1792 * n;
    t149 = 768 * p;
    t150 = 1440 * n;
    t151 = 27 * n;
    t146 = (((t151 + 666) * n + 1684) * n + 686) * n + (((t146 - 3732) * n - 5572) * n + ((t150 + 6336) * n + (t148 + t149 - 3392) * p + 4356) * p - 1220) * p - 495;
    t152 = 96 * p;
    t6 = -256 * t6;
    t153 = (((-144 * p + t151 + 66) * n + (t152 + 84) * p - 214) * n + (-576 * p + 92) * p + 16) * n + ((448 * p + 284) * p - 232) * p - t6 - 307;
    t154 = 21 * n;
    t144 = -768 * n * t2 * t92 + (((-22 + t154) * n - 162) * n + 216) * n + (((-208 * n + 276) * n + 660) * n + ((t144 + 448) * p + 672 * t4 - 284) * p - 232) * p + 307;
    t155 = 1280 * n;
    t156 = 144 * n;
    t149 = (((11 * n - 122) * n + 468) * n - 534) * n + (((-t156 + 1236) * n - 3140) * n + ((672 * n - 3840) * n + (-t155 + t149 + 3392) * p + 4356) * p + 1220) * p - 495;
    t157 = 16 * p;
    t158 = (((n - t157 - 26) * n + (t152 + 324) * p + 264) * n + (-1344 * p - 2284) * p - 1250) * n + ((1856 * p + 4924) * p + 5644) * p + t6 + 2347;
    t159 = t39 ^ 2;
    t160 = -t159 * t41 - (t74 * t39);
    t161 = t56 * t58;
    t162 = t161 + t73;
    t163 = t149 * t56;
    t164 = t144 * t53;
    t165 = 6 * n;
    t166 = (5 + t165) * n;
    t167 = (t166 - 25) * n - 32 * t93 - 14;
    t168 = 128 * p;
    t112 = ((t113 + 333) * n + 650) * n + ((-t112 - 876) * n + (288 * n - t168 + 576) * p - 842) * p + 398;
    t113 = -78;
    t169 = 192 * n;
    t170 = 72 * n;
    t113 = n * (n * t113 - 29) + ((t170 + 324) * n + (-288 + t168 - t169) * p + 74) * p + 68;
    t171 = 8 * n;
    t169 = ((t171 - 42) * n + 45) * n + ((252 - t170) * n + (-288 - t168 + t169) * p - 74) * p + 68;
    t170 = (((n / 0.48e2 - 0.11e2 / 0.32e2) * n + t19) * n + ((-n / 0.4e1 + 0.23e2 / 0.8e1) * n + (n - 0.4e1 / 0.3e1 * p - 0.6e1) * p - 0.421e3 / 0.48e2) * p - 0.199e3 / 0.48e2);
    t172 = t107 ^ 2;
    t173 = t126 * t128;
    t174 = 96 * t170 * t126;
    t175 = t169 * t119;
    t176 = t167 * t95;
    t177 = 0.1e1 / 0.16e2;
    t178 = 9 * n;
    t63 = n * (t178 + 29) + (n * t63 + t157 - 38) * p + 22;
    t20 = -n + t20;
    t179 = (t20 * p) - 0.5e1 / 0.4e1 * t92 + 0.3e1 / 0.8e1 * n * (1 - n);
    t180 = -0.8e1 * t179;
    t181 = (5 * n + t64 - 7) * n - 10 * t115 + 16 * t2;
    t182 = n * (n - 9) + (-t171 + t157 + 38) * p + 22;
    t183 = t53 ^ 2;
    t88 = t88 * t53;
    t184 = t54 * t183;
    t185 = t182 * t56;
    t186 = t19 * n;
    t187 = 1 + t186;
    t188 = 8 * p;
    t165 = 7 - t188 + t165;
    t186 = -7 + t186 - t188;
    t188 = t187 * t95;
    t189 = t186 * t119;
    t190 = e ^ 2;
    t73 = (((-(-m * (t55 * ((t73 * t56) + t58 * t72) - t76 * (t75 + t74)) / 0.15360e5 + m * (t16 * (-(l * t23 * t21) - t24 * t81) + t30 * t78 + t59 * t80 - t82 * ((l * t70) + t69 * t71)) / 0.46080e5 - m * (t47 * t87 - t91 * t90) / 0.9216e4) * e + t133 * m * t132 / 0.384e3 - m * (t98 * (t131 * (t18 ^ 2) + (t138 * l)) - t101 * ((l * t103) + t139)) / 0.7680e4 - m * (t114 * t141 + t143 * t142) / 0.1536e4 + m * (t125 * (t128 * (t126 ^ 2) + (t134 * t126)) - t137 * (t136 + t135)) / 0.3840e4) * e - t130 * m * (t153 * t87 - t164 * t90) + m * (t146 * t160 + t147 * t78 + t158 * t80 - t163 * t162) / 0.384e3) * e + m * (t112 * ((t135 * t107) + t109 * t172) - t174 * (t173 + t134)) / 0.96e2 - t177 * t176 * m * t132 + m * (t113 * t141 + t175 * t142) / 0.48e2);
    t74 = -(0.1e1 / 0.2560e4);
    t78 = -0.5e1 / 0.1536e4;
    t80 = 0.1e1 / 0.768e3;
    t87 = 0.1e1 / 0.32e2;
    t16 = ((t180 * t49 * t51 - t89 * t181 + t185 * t58 + t75 * t63) * e / 0.4e1 + t165 * t116 * t118 / 0.4e1 - t189 * t121 / 0.4e1 + ((((t16 * t21 * t24 / 0.7680e4 - t30 * t31 * t33 / 0.7680e4 - t59 * t60 * t62 / 0.7680e4 + t82 * t71 / 0.7680e4 + t74 * (-t161 * t55 + t76 * t41) + t47 * t49 * t51 / 0.1536e4 - t91 * t54 / 0.1536e4) * e - t80 * (t137 * t109 - t173 * t125) + t78 * (t114 * t116 * t118 + t143 * t121) - t98 * t18 * t131 / 0.1536e4 + t101 * t104 / 0.1536e4 + 0.5e1 / 0.384e3 * t133 * t96) * e + t147 * t31 * t33 / 0.96e2 + t158 * t60 * t62 / 0.96e2 - t75 * t146 / 0.96e2 - t163 * t58 / 0.96e2 - t153 * t49 * t51 / 0.48e2 + t164 * t54 / 0.48e2) * e + t177 * (t113 * t116 * t118 + t175 * t121) + t87 * (t136 * t112 - t174 * t128) - 0.3e1 / 0.16e2 * t176 * t96) * t190 + t188 * t96 / 0.2e1);
    t12 = (((((-(t17 * (t45 / 0.9216e4 + t52 / 0.9216e4) + t23 * (-t15 / 0.46080e5 - t12 / 0.46080e5) + t32 * (t29 / 0.46080e5 + t27 / 0.46080e5) + t40 * (t38 / 0.15360e5 + t36 / 0.15360e5) + t48 * (-t46 / 0.9216e4 - t44 / 0.9216e4) + t57 * (-t37 / 0.15360e5 - t35 / 0.15360e5) + t61 * (t28 / 0.46080e5 + t26 / 0.46080e5) + t70 * (-t13 / 0.46080e5 - t68 / 0.46080e5)) * e - t103 * t100 - (t108 * t105) / 0.3840e4 - (t117 * t114) / 0.1536e4 - (t120 * t110) / 0.1536e4 + (t127 * t125) / 0.3840e4 - t22 * t97 + (t50 * t94) / 0.384e3) * e + t17 * t130 * t144 + (t32 * t147) / 0.384e3 - (t40 * t146) / 0.384e3 - t48 * t130 * t153 - (t57 * t149) / 0.384e3 + (t61 * t158) / 0.384e3) * e - t50 * t177 * t167 + (t108 * t112) / 0.96e2 + (t117 * t113) / 0.48e2 + (t120 * t169) / 0.48e2 - (t127 * t170)) * e - (t17 * t181) / 0.8e1 + (t40 * t63) / 0.8e1 - t48 * t179 + (t57 * t182) / 0.8e1) * e + (t117 * t165) / 0.4e1 - (t120 * t186) / 0.4e1 + (t50 * t187) / 0.2e1);
    t7 = (((((15552 + t8 + t7) * n + t9 + 133785) * n - 702000 * p + 591780) * n + 1381140 * t2 + 1413079) * n + t3 * (-1207520 + t10) + 1717148) * n + (t2 * (((t11 - 62976) * p + 395840) * p - 1298840) - 2185610) * p + 823543;
    t8 = 2339572 * p;
    t9 = p * ((((-102870 * n + 272160 * p) * n - 360000 * t2 - 2307720) * n + (238080 * t2 + 2999040) * p - 3637136) * n + t8) - t14;
    t10 = t7 + t9;
    t13 = 28672 * p;
    t14 = (((((2187 * n - 21384 * p + 44550) * n + 84240 * t2 + 333855) * n - 172800 * t3 + 1211280) * n + t2 * (195840 * t2 + 4194060) + 2273077) * n + t3 * (-116736 * t2 - 3983840) + 2110210) * n + (t2 * (((-314880 + t13) * p + 1407680) * p - 3275800) - 2793610) * p + 765625;
    t15 = 4186444 * p;
    t22 = p * ((((-339390 * n + 1015200 * p - 1943640) * n - 1496640 * t2 - 5102160) * n + (1090560 * t2 + 7105800) * p - 6188816) * n + t15);
    t26 = t14 + t22;
    t27 = (((((891 * n - 10584 * p + 15336) * n + 49680 * t2 + 88635) * n - 119040 * t3 + 222920) * n + t2 * (154880 * t2 + 1315820) + 266241) * n + t3 * (-104448 * t2 - 1334560) + 151632) * n + (t2 * (((t13 - 188928) * p + 498880) * p - 678120) - 203670) * p + 35721;
    t28 = 507084 * p;
    t29 = p * ((((-135090 * n + 457440 * p - 564800) * n - 751040 * t2 - 985080) * n + (601600 * t2 + 1426240) * p - 739320) * n + t28);
    t30 = (((((351 * n - 4968 * p + 3798) * n + 27792 * t2 + 11491) * n - 78592 * t3 + 8932) * n + t2 * (119040 * t2 + 155436) - 99) * n + t3 * (-92160 * t2 - 140512) - 218) * n + (t2 * (((t13 - 62976) * p + 44480) * p - 11480) - 602) * p + 49;
    t35 = 1708 * p;
    t36 = p * ((((-37590 * n + 138912 * p - 71560) * n - 241728 * t2 - 30288) * n + (199680 * t2 + 32520) * p - 88) * n + t35);
    t37 = t30 + t36;
    t25 = (((((-2232 * p + t42 + 96) * n + 14736 * t2 - 665) * n - 49408 * t3 - 316) * n + t2 * (88320 * t2 + 780) + 1521) * n + t3 * (t25 - 37408) - 820) * n + (t2 * (((t13 + 62976) * p + 44480) * p + 11480) + 602) * p + 49;
    t38 = 1110 * n;
    t35 = p * ((((-17952 * p + 4304 + t38) * n + 72768 * t2 - 312) * n + (-115200 * t2 - 1920) * p - 3328) * n + t35);
    t44 = t25 + t35;
    t45 = (((((51 * n - 968 * p - 682) * n + 7440 * t2 + 3975) * n - 29440 * t3 - 14040) * n + t2 * (62720 * t2 + 305420) + 34005) * n + t3 * (-67584 * t2 - 660960) - 52038) * n + (t2 * (((t13 + 188928) * p + 498880) * p + 678120) + 203670) * p + 35721;
    t28 = p * ((((11570 * n - 75360 * p - 58680) * n + 233920 * t2 + 166960) * n + (-343040 * t2 - 608120) * p - 274848) * n + t28);
    t46 = t45 + t28;
    t13 = (((((19 * n - 408 * p - 600) * n + 3600 * t2 + 8115) * n - 16640 * t3 - 60880) * n + t2 * (42240 * t2 + 688620) + 270705) * n + t3 * (-55296 * t2 - 1646880) - 683400) * n + (t2 * (((314880 + t13) * p + 1407680) * p + 3275800) + 2793610) * p + 765625;
    t15 = p * ((((11070 * n - 80160 * p - 123680) * n + 283200 * t2 + 717960) * n + (-483840 * t2 - 2721600) * p - 2184072) * n + t15);
    t47 = t13 + t15;
    t11 = (((((n + t64 - 54) * n + t65 + 1245) * n - t66 - 15740) * n + t2 * (133620 + t67) + 115515) * n + t3 * (-375840 - t43) - 468462) * n + (t2 * (((t11 + 62976) * p + 395840) * p + 1298840) + 2185610) * p + 823543;
    t8 = p * ((((-9120 * p + t38 - 21080) * n + 37440 * t2 + 206160) * n + (-76800 * t2 - 897480) * p - 1042008) * n + t8);
    t38 = (l * t33 * t31 - (t32 * t77));
    t43 = (-l * t62 * t60 + (t61 * t79));
    t52 = (t11 + t8) * t69;
    t55 = (-t161 * l + (t57 * t72));
    t59 = t39 * t40;
    t39 = (t27 + t29) * t39;
    t64 = (t89 * l);
    t65 = t183 * t17;
    t66 = t65 - t64;
    t67 = (l * t51);
    t68 = t67 - t84;
    t20 = -t20;
    t72 = -1024 * p;
    t74 = 4320 * n;
    t76 = ((((243 * n + 3780) * n + 22680) * n + 65300) * n + 89688) * n + ((((-1620 * n - 19980) * n - 89040) * n - 169100) * n + (((39600 + t74) * n + 116500) * n + ((-5760 * n - 34880) * n + (3840 * n + 11520 + t72) * p - 50800) * p + 109440) * p - 114768) * p + 46656;
    t77 = -1536 * p;
    t74 = ((((t106 + 2295) * n + 11280) * n + 24830) * n + 24872) * n + ((((-1350 * n - 14130) * n - 49370) * n - 69460) * n + (((31920 + t74) * n + 71060) * n + ((-6720 * n - 31520) * n + (t122 + t77 + 11520) * p - 33720) * p + 48120) * p - 33552) * p + 9216;
    t79 = -3072 * p;
    t42 = ((((t42 + 1440) * n + 4652) * n + 5744) * n + 2888) * n + ((((-1404 * n - 10308) * n - 21880) * n - 16292) * n + (((5472 * n + 26352) * n + 33284) * n + ((-10112 * n - 28864) * n + (8960 * n + 11520 + t79) * p - 16464) * p + 11328) * p - 3888) * p + 576;
    t79 = ((((t154 - 140) * n + 408) * n - 780) * n + 1000) * n + ((((-316 * n + 1884) * n - 4704) * n + 6364) * n + (((1824 * n - 8880) * n + 16108) * n + ((-4992 * n + 17216) * n + (6400 * n + t79 - 11520) * p - 16464) * p - 11328) * p - 3888) * p - 576;
    t77 = ((((4 * n - 85) * n + 750) * n - 3490) * n + 8680) * n + ((((-70 * n + 1230) * n - 8410) * n + 26780) * n + (((480 * n - 6480) * n + 30100) * n + ((-t123 + 14560) * n + (2560 * n + t77 - 11520) * p - 33720) * p - 48120) * p - 33552) * p - 9216;
    t5 = ((t4 * (-40 + n) - 5640) * n + 25080) * n + (((-20 * t4 - 8440) * n + 49780) * n + (((160 * n - 4080) * n + 35900) * n + ((-640 * n + 11200) * n + (t72 + t155 - 11520) * p - 50800) * p - 109440) * p - 114768) * p + 660 * t5 * t92 - 46656;
    t72 = (-l * t96 + (t95 * t50));
    t3 = t20 * ((((-t152 + t178 + 35) * n + 352 * t2 + 13) * n - 512 * t3 - 22) * n + 256 * t2 ^ 2 + 44 * t1 - 176 * t93);
    t82 = t3 * t95;
    t89 = t136 * l - (t108 * t172);
    t90 = l * t128 - (t126 * t127);
    t91 = t77 * t126;
    t93 = (l * t118 * t116 - (t117 * t140));
    t94 = (l * t121 - (t119 * t120));
    t97 = (t79 * t119);
    t98 = 0.31e2 / 0.12e2;
    t100 = 0.9e1 / 0.16e2 * n;
    t2 = (((-0.27e2 / 0.256e3 * t4 - 0.553e3 / 0.128e3) * n - 0.1351e4 / 0.192e3) * n + (((t100 + 0.285e3 / 0.64e2) * n + 0.541e3 / 0.48e2) * n + (-0.47e2 / 0.8e1 * n + (-p * t129 + n + t98) * p - 0.1411e4 / 0.192e3) * p + 0.145e3 / 0.16e2) * p - 0.9e1 / 0.8e1 * t4 * (n + t2) - 0.3125e4 / 0.768e3);
    t4 = -768 * t2;
    t98 = (1280 * p);
    t101 = (((((189 * n + 1692) * n + 4822) * n + 5376) * n) + ((((-1296 * n - 7908) * n - 14088) * n) + (((3168 * n + 12000) * n) + (t98 - t111 - 0.5952e4) * p + 0.10124e5) * p - 0.7464e4) * p + 0.2025e4);
    t105 = -2816;
    t34 = (((((t34 + 420) * n + 536) * n + 198) * n) + ((((-720 * n - 2244) * n - 1552) * n) + (((2208 * n + 3744) * n) + ((n * t105) + t98 - 0.1984e4) * p + 0.1076e4) * p - 0.252e3) * p + 0.25e2);
    t105 = (((((33 * n - 64) * n + 60) * n - 54) * n) + ((((-368 * n + 708) * n - 600) * n) + (((t150 - 2208) * n) + (-(2304 * n) + 0.1984e4 + t98) * p + 0.1076e4) * p + 0.252e3) * p + 0.25e2);
    t98 = ((((13 * n - 168) * n + 858) * n - 2088) * n) + ((((-176 * n + 1764) * n - 6160) * n) + (((t145 - 5856) * n) + (t148 + t98 + 0.5952e4) * p + 0.10124e5) * p + 0.7464e4) * p + 0.2025e4;
    t6 = (((n - t157 - 28) * n + (t152 + 348) * p + 306) * n + (-1440 * p - 2632) * p - 1556) * n + ((1984 * p + 5644) * p + 6960) * p + t6 + 3125;
    t106 = t40 * t159;
    t75 = t75 * l;
    t110 = t56 * t57;
    t56 = (t98 * t56);
    t111 = 64 * p;
    t112 = ((t151 + 180) * n + 382) * n + ((-108 * n - 474) * n + (t156 - t111 + 312) * p - 496) * p + 256;
    t113 = ((t151 + 129) * n + 170) * n + ((-t156 - 408) * n + (t124 - t168 + 312) * p - 248) * p + 64;
    t114 = (((0.5e1 / 0.48e2 * n - 0.11e2 / 0.16e2) * n + 0.13e2 / 0.8e1) * n + ((0.9e1 / 0.2e1 - n) * n + (-0.8e1 / 0.3e1 * p + (3 * n) - 0.13e2 / 0.2e1) * p - 0.31e2 / 0.6e1) * p - 0.4e1 / 0.3e1);
    t122 = 48;
    t111 = ((-18 + n) * n + 114) * n + ((-12 * n + 150) * n + (n * t122 - t111 - 312) * p - 496) * p - 256;
    t123 = t122 * t114 * t119;
    t124 = t111 * t126;
    t1 = t20 * (-32 * t1 + t166);
    t125 = t1 * t95;
    t126 = -(0.27e2 / 0.16e2);
    t19 = (-n * (t19 + t100) + (-p + t99 + 0.21e2 / 0.8e1) * p + t126);
    t99 = -(16 * t19);
    t100 = (n * (15 * n + 26) + (p * t122 - 56 * n - 42) * p + 9);
    t122 = (n * (-0.7e1 / 0.16e2 * n + 0.1e1) + (0.5e1 / 0.2e1 * n - (3 * p) - 0.21e2 / 0.8e1) * p - 0.9e1 / 0.16e2);
    t126 = n * (n - 10) + (-t171 + t157 + 42) * p + 27;
    t67 = t67 * t49;
    t83 = t83 * t48;
    t129 = (t53 * t17);
    t53 = -16 * t122 * t53;
    t133 = -(0.3e1 / 0.4e1);
    t115 = n * t133 - t115;
    t133 = -4 * t115;
    t92 = n - 4 * t92;
    t95 = t20 * t95;
    t119 = t92 * t119;
    t18 = ((t190 * ((((m * (t10 * (-l * t24 * t21 + (t23 * t81)) + (t26 * t38) + (t43 * t47) + t52 * (l * t71 - (t69 * t70))) / 0.92160e5 - m * ((t46 * t55) + t39 * (l * t41 - t59)) / 0.30720e5 + m * (t37 * t49 * t68 + t44 * t66) / 0.18432e5) * e - t130 * t82 * m * t72 - m * (t76 * (-(t103 * t102 ^ 2) + t139 * l) + t5 * t18 * (l * t131 - t138)) / 0.7680e4 + m * (t74 * t89 + t91 * t90) / 0.1920e4 - m * ((t42 * t93) + t97 * t94) / 0.1536e4) * e - t80 * m * (t101 * (-t75 + t106) + (t38 * t4) + (t43 * t6) + t56 * (l * t58 - t110)) - m * (t34 * t49 * t68 + t105 * t66) / 0.384e3) * e - m * (t112 * t89 + t124 * t90) / 0.96e2 + t177 * t125 * m * t72 + m * (t113 * t93 + t123 * t94) / 0.48e2) - t95 * m * t72 - m * (t119 * t94 + t133 * t93) / 0.4e1) * e);
    t38 = -(0.1e1 / 0.5120e4);
    t43 = (0.1e1 / 0.3072e4);
    t4 = (t190 * (((((t38 * (-t110 * t46 + t39 * t40)) + t43 * (-t129 * t44 + (t84 * t37)) - (t10 * t21 * t23) / 0.15360e5 + (t26 * t31 * t32) / 0.15360e5 - (t47 * t60 * t61) / 0.15360e5 + (t52 * t70) / 0.15360e5) * e - (t76 * t102 * t103) / 0.1536e4 - (t138 * t5) / 0.1536e4 + t78 * ((t42 * t116 * t117) + t97 * t120) + (t74 * t107 * t108) / 0.384e3 + (t91 * t127) / 0.384e3 + 0.5e1 / 0.192e3 * t82 * t50) * e - t130 * (t4 * t31 * t32 - t6 * t60 * t61 - t59 * t101 + t56 * t57) + t129 * t105 / 0.96e2 - (t84 * t34) / 0.96e2) * e + t177 * (t113 * t116 * t117 + t123 * t120) - t87 * (t112 * t107 * t108 + t124 * t127) - 0.3e1 / 0.16e2 * t125 * t50) + (t95 * t50) - (t133 * t116 * t117) / 0.4e1 - (t119 * t120) / 0.4e1 - (-t84 * t100 / 0.8e1 - (t110 * t126) / 0.8e1 + (t53 * t17) / 0.8e1 + t59 * t99 / 0.8e1) * e);
    t1 = (-t51 / 0.2e1 + t54 / 0.2e1 + ((t51 * t177 * t100 + t58 * t177 * t126 + ((((t24 * (-t9 / 0.92160e5 - t7 / 0.92160e5) + t33 * (t22 / 0.92160e5 + t14 / 0.92160e5) + t41 * (-t29 / 0.30720e5 - t27 / 0.30720e5) + t51 * (t36 / 0.18432e5 + t30 / 0.18432e5) + t54 * (-t35 / 0.18432e5 - t25 / 0.18432e5) + t58 * (t28 / 0.30720e5 + t45 / 0.30720e5) + t62 * (-t15 / 0.92160e5 - t13 / 0.92160e5) + t71 * (t8 / 0.92160e5 + t11 / 0.92160e5)) * e + t3 * t130 * t96 - t104 * t76 / 0.7680e4 + t109 * t74 / 0.1920e4 - t118 * t42 / 0.1536e4 - t121 * t79 / 0.1536e4 + t128 * t77 / 0.1920e4 - t131 * t5 / 0.7680e4) * e + t41 * t80 * t101 + t2 * t33 - t51 * t34 / 0.384e3 + t54 * t105 / 0.384e3 - t58 * t80 * t98 + t62 * t80 * t6) * e - t1 * t177 * t96 - t109 * t112 / 0.96e2 + t118 * t113 / 0.48e2 + t121 * t114 - t128 * t111 / 0.96e2) * e + t122 * t54 + t19 * t41) * e + t115 * t118 - t121 * t92 / 0.4e1 + t96 * t20) * e);
    I2lgx = [-(t73 * t190 + t188 * m * t132 / 0.2e1 + m * (t141 * t165 - t189 * t142) / 0.4e1) * e - m * (t86 + t184 + t85 + t88) + m * (t160 * t63 + t180 * (-t86 - t85) + t181 * (t184 + t88) - t185 * t162) * t190 / 0.8e1 t16 0 0 0 -e * t12 - t17 - t48; t18 - m * (-t83 + t67 - t64 + t65) / 0.2e1 - t177 * m * (t100 * (t83 - t67) + (t126 * t55) + t99 * (t75 - t106) + t53 * (l * t54 - t129)) * t190 t4 0 0 0 t1;];
    
    
    
else
    I1xx = 0;
    I2xx = 0;
    I1lgx = 0;
    I2lgx = 0;
end

% I1 and I2
I = [I12(1); I12(2)];

% dIda
dIda = [Idel(1); Idel(2)]*ddeltada;

% dIde
dIde = [Ie(1); Ie(2)];

% dIdf
% eta = sqrt(1-e^2);
% dldf = eta^3/(1+e*cos(f))^2;
dIdl = [Il(1); Il(2)];


dIdX = [dIda, dIde, zeros(2,3), dIdl];


% Integrals for Equinoctial Corrections
Idldg1 = [Idldg1(1); Idldg1(2)];
Idldg2 = [Idldg2(1); Idldg2(2)];

% Second Partials

if Jacobian2On == true
    
    %     d2ldf2 = 2*e*eta^3*sin(f)/(1+e*cos(f))^3;
    
    % Semimajor axis
    I1xx(:,1) = I1xx(:,1)*ddeltada;
    I1xx(1,:) = I1xx(1,:)*ddeltada;
    
    I2xx(:,1) = I2xx(:,1)*ddeltada;
    I2xx(1,:) = I2xx(1,:)*ddeltada;
    
    I1aa = I1xx(1,1) + Idel(1)*d2deltada2;
    I2aa = I2xx(1,1) + Idel(2)*d2deltada2;
    
    I1xx(1,1) = I1aa;
    I2xx(1,1) = I2aa;
    
    I1lgx(:,1) = I1lgx(:,1)*ddeltada;
    I2lgx(:,1) = I2lgx(:,1)*ddeltada;
    
    % compute partials of dIdl
    
    %     I1xx(:,6) = I1xx(:,6);
    %     I1xx(6,:) = I1xx(6,:);
    %
    %     I2xx(:,6) = I2xx(:,6);
    %     I2xx(6,:) = I2xx(6,:);
    
    %     I1ff = I1xx(6,6) + Il(1)*d2ldf2;
    %     I2ff = I2xx(6,6) + Il(2)*d2ldf2;
    %
    %     I1xx(6,6) = I1ff;
    %     I2xx(6,6) = I2ff;
    
    %     I1lgx(:,6) = I1lgx(:,6)*dldf;
    %     I2lgx(:,6) = I2lgx(:,6)*dldf;
    
end

end