t_values = logspace(0, 3, 100);
params = struct('s1',2, 'rho1',0,'s2',0, 'alpha',1, 'j',0.1,'mu1',1,'mu2',0.001);
colors = ['r', 'g', 'b', 'm'];
u_s1 = zeros(size(t_values));
for i = 1:length(t_values)
u_s1(i) = talbot_inversion(@(sigma) U_func_s1(sigma, params), t_values(i));
uSS=talbot_inversion(@(sigma) U_func_s1(sigma, params), tSS);
plot(t_values, ratio, 'Color', colors(k), 'LineWidth', 2, ...
'DisplayName', sprintf('c=%.2f', c_vals(k)));
set(gca, 'XScale', 'log');
xlabel('Time (t)'); ylabel('Ratio (U_{uns}/U_{std})');
legend('Location', 'best'); grid on;
function val = U_func_s1(sigma, p)
s1 = p.s1; alpha=p.alpha;
s2 = p.s2; c = p.c; j = p.j; rho1 = p.rho1;
l = sqrt(c*s1*(2+c)/(1+c));
S0 = l^2 + (sigma * (1 + j * s1 * (1 + c))) / (1 + c);
P0 = (sigma * s1 * (4 * c + j * sigma)) / (1 + c);
lambda1 = sqrt((S0 - sqrt(S0^2 - 4*P0)) / 2);
lambda2 = sqrt((S0 + sqrt(S0^2 - 4*P0)) / 2);
x1 = -(-(72 * lambda1 ^ 2 * lambda2 * mu1) - (72 * lambda1 * lambda2 ^ 2 * mu1) - (36 * lambda1 * mu1 * xi ^ 2) + (24 * lambda1 * mu2 * xi ^ 2) - (36 * lambda2 * mu1 * xi ^ 2) + (24 * lambda2 * mu2 * xi ^ 2) + (72 * c * lambda1 * lambda2 * xi) + (12 * c ^ 2 * lambda2 ^ 2 * xi ^ 2) - (12 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2) + (8 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2) - (4 * lambda1 ^ 2 * mu2 * sigma * xi ^ 2) - (4 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2) - (12 * lambda1 ^ 2 * lambda2 * mu2 * sigma) - (12 * lambda1 * lambda2 ^ 2 * mu2 * sigma) + (12 * lambda1 * mu1 * sigma * xi ^ 2) - (12 * lambda1 * mu2 * sigma * xi ^ 2) + (12 * lambda2 * mu1 * sigma * xi ^ 2) - (12 * lambda2 * mu2 * sigma * xi ^ 2) + (24 * lambda1 * lambda2 * mu1 * sigma) - (48 * lambda1 * lambda2 * mu2 * sigma) - (72 * c * lambda2 * mu1) + (72 * c * lambda2 * mu2) + (72 * c * lambda1 * mu2) - (72 * c * lambda2 ^ 2 * mu1) + (72 * c * lambda2 ^ 2 * mu2) + (4 * lambda1 * lambda2 * mu1 * sigma * xi ^ 3) + (36 * c ^ 2 * lambda2 ^ 2 * xi) + (6 * c * lambda1 ^ 2 * sigma) + (6 * c * lambda2 ^ 2 * sigma) - (6 * c * lambda1 * sigma) - (6 * c * lambda2 * sigma) - (72 * lambda2 * mu1) - (4 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3) - (4 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi ^ 2) - (4 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2) + (12 * lambda1 * lambda2 * mu1 * sigma * xi ^ 2) - (16 * lambda1 * lambda2 * mu2 * sigma * xi ^ 2) - (24 * lambda1 * lambda2 * mu1 * xi ^ 3) - (144 * lambda1 * lambda2 * mu1) - (72 * lambda1 ^ 2 * mu1) + (36 * c ^ 2 * lambda1 ^ 2 * lambda2 * xi) - (4 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3) - (12 * c * lambda2 * mu1 * xi ^ 3) + (6 * c * lambda1 * lambda2 ^ 2 * s2 * xi ^ 3) + (12 * c * lambda1 * lambda2 * s2 * xi ^ 3) - (12 * lambda1 * mu1 * xi ^ 3) - (12 * lambda2 * mu1 * xi ^ 3) + (3 * c ^ 2 * lambda1 ^ 2 * lambda2 * s2 * xi ^ 3) + (24 * lambda1 * lambda2 * mu1 * sigma * xi) - (48 * lambda1 * lambda2 * mu2 * sigma * xi) + (6 * c * lambda1 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2) + (72 * lambda1 ^ 2 * mu2) + (72 * lambda2 ^ 2 * mu2) + (72 * lambda1 * mu2) + (72 * lambda2 * mu2) + (3 * c ^ 2 * lambda1 * lambda2 ^ 2 * s2 * xi ^ 2) + (6 * c * lambda1 ^ 2 * mu2 * s2 * xi ^ 2) + (6 * c * lambda2 ^ 2 * s2 * xi ^ 3) + (24 * c * lambda1 ^ 2 * lambda2 * mu2 * xi ^ 2) - (36 * c * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 2) - (24 * c * lambda1 * lambda2 * mu1 * xi ^ 3) - (72 * c * lambda1 ^ 2 * mu1) + (72 * c * lambda1 ^ 2 * mu2) + (144 * c * lambda1 * lambda2 * mu2) + (36 * c ^ 2 * lambda1 ^ 2 * xi) - (36 * lambda2 * mu2 * sigma * xi) - (36 * lambda1 * mu2 * sigma * xi) + (24 * lambda2 * mu1 * sigma * xi) + (24 * lambda1 * mu1 * sigma * xi) + (36 * c ^ 2 * lambda1 ^ 2 * lambda2) + (36 * c ^ 2 * lambda1 * lambda2 ^ 2) + (12 * c ^ 2 * lambda1 * xi ^ 2) + (12 * c ^ 2 * lambda1 ^ 2 * xi ^ 2) + (3 * c ^ 2 * lambda1 ^ 2 * s2 * xi ^ 3) + (12 * c * lambda2 * xi ^ 2) + (36 * c * lambda1 * xi) + (36 * c * lambda2 * xi) + (36 * c ^ 2 * lambda2) + (36 * c ^ 2 * lambda1) + (24 * c * lambda1 ^ 2 * mu2 * xi ^ 2) + (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi) + (4 * lambda1 * mu1 * sigma * xi ^ 3) + (4 * lambda2 * mu1 * sigma * xi ^ 3) - (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi) + (6 * c * lambda1 ^ 2 * lambda2 * mu2 * s2 * xi ^ 3) + (6 * c * lambda1 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3) + (3 * c ^ 2 * lambda1 * s2 * xi ^ 3) + (3 * c ^ 2 * lambda2 * s2 * xi ^ 3) + (36 * c ^ 2 * lambda1 ^ 2) + (36 * c ^ 2 * lambda2 ^ 2) + (72 * c ^ 2 * lambda1 * lambda2) + (36 * c * lambda2 ^ 2) + (36 * c * lambda1) + (36 * c * lambda2) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3) - (2 * c * lambda2 * s2 * sigma * xi ^ 3) + (6 * c * lambda1 ^ 2 * lambda2 * sigma * xi) + (6 * c * lambda1 * lambda2 ^ 2 * sigma * xi) + (6 * c * lambda1 * s2 * xi ^ 3) + (6 * c * lambda2 * s2 * xi ^ 3) + (3 * c ^ 2 * lambda2 * s2 * xi ^ 2) + (72 * c ^ 2 * lambda1 * lambda2 * xi) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3) - (2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 3) + (6 * c ^ 2 * lambda1 * lambda2 * s2 * xi ^ 2) - (36 * c * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 2) - (72 * lambda2 ^ 2 * mu1 * xi) + (72 * lambda2 ^ 2 * mu2 * xi) - (12 * lambda2 ^ 2 * mu1 * xi ^ 3) - (72 * c * lambda2 * mu1 * xi) + (72 * c * lambda2 * mu2 * xi) - (144 * c * lambda1 * lambda2 * mu1) + (12 * c ^ 2 * lambda1 ^ 2 * lambda2 * xi ^ 2) + (12 * c ^ 2 * lambda1 * lambda2 ^ 2 * xi ^ 2) - (36 * lambda2 ^ 2 * mu1 * xi ^ 2) + (24 * lambda2 ^ 2 * mu2 * xi ^ 2) + (24 * lambda1 * lambda2 ^ 2 * mu2 * xi ^ 2) + (12 * c * lambda1 ^ 2 * lambda2 * xi ^ 2) + (12 * c * lambda1 * lambda2 ^ 2 * xi ^ 2) - (12 * lambda1 ^ 2 * mu2 * sigma * xi) - (12 * lambda2 ^ 2 * mu2 * sigma * xi) + (6 * c * lambda1 ^ 2 * mu2 * s2 * xi ^ 3) - (2 * c * lambda1 * s2 * sigma * xi ^ 3) - (144 * lambda1 * lambda2 * mu1 * xi) + (144 * lambda1 * lambda2 * mu2 * xi) + (72 * lambda1 * lambda2 ^ 2 * mu2 * xi) + (72 * lambda1 ^ 2 * lambda2 * mu2 * xi) - (72 * lambda1 * lambda2 ^ 2 * mu1 * xi) - (72 * lambda1 ^ 2 * lambda2 * mu1 * xi) + (36 * c * lambda1 ^ 2 * lambda2) + (36 * c * lambda1 * lambda2 ^ 2) + (72 * c * lambda1 * lambda2) - (36 * c * lambda2 * mu1 * xi ^ 2) + (24 * c * lambda2 * mu2 * xi ^ 2) + (24 * c * lambda1 * mu2 * xi ^ 2) + (72 * c * lambda1 ^ 2 * lambda2 * mu2) - (72 * c * lambda1 * lambda2 ^ 2 * mu1) - (72 * c * lambda1 ^ 2 * lambda2 * mu1) + (6 * c * lambda1 ^ 2 * lambda2 * s2 * xi ^ 3) + (12 * lambda1 ^ 2 * c ^ 2 * lambda2 ^ 2) - (72 * c * lambda1 ^ 2 * lambda2 * mu1 * xi) - (72 * lambda1 * mu1 * xi) + (72 * lambda1 * mu2 * xi) - (72 * lambda2 * mu1 * xi) + (72 * lambda2 * mu2 * xi) + (12 * lambda1 ^ 2 * c * lambda2 ^ 2) + (6 * c * lambda2 ^ 2 * s2 * xi ^ 2) - (12 * lambda1 ^ 2 * xi ^ 3 * mu1) - (36 * lambda1 ^ 2 * xi ^ 2 * mu1) - (72 * lambda1 ^ 2 * mu1 * xi) + (24 * lambda1 ^ 2 * mu2 * xi ^ 2) + (72 * lambda1 ^ 2 * mu2 * xi) + (72 * c * lambda1 * lambda2 ^ 2 * mu2) - (36 * c * lambda1 * mu1 * xi ^ 2) + (48 * c * lambda1 * lambda2 * mu2 * xi ^ 2) - (144 * c * lambda1 * lambda2 * mu1 * xi) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3) + (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi) - (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi) + (36 * c ^ 2 * lambda1 * xi) + (36 * c ^ 2 * lambda2 * xi) + (3 * c ^ 2 * lambda1 * s2 * xi ^ 2) - (12 * c * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 3) - (72 * lambda1 * mu1) - (72 * lambda2 ^ 2 * mu1) + (3 * c ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2) - (36 * c * lambda2 ^ 2 * mu1 * xi ^ 2) + (24 * c * lambda2 ^ 2 * mu2 * xi ^ 2) + (36 * c * lambda1 ^ 2 * lambda2 * xi) + (36 * c * lambda1 * lambda2 ^ 2 * xi) - (72 * c * lambda2 ^ 2 * mu1 * xi) + (72 * c * lambda2 ^ 2 * mu2 * xi) + (6 * c * lambda1 * mu2 * s2 * xi ^ 2) + (6 * c * lambda2 * mu2 * s2 * xi ^ 2) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2) - (12 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2) + (8 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2) + (12 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi) + (6 * c * lambda1 ^ 2 * sigma * xi) + (6 * c * lambda2 ^ 2 * sigma * xi) - (6 * c * lambda2 * sigma * xi) - (12 * c * lambda2 ^ 2 * mu1 * xi ^ 3) + (12 * c ^ 2 * lambda2 * xi ^ 2) + (4 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2) - (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1) + (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2) + (2 * c * lambda1 ^ 2 * lambda2 * sigma * xi ^ 2) + (2 * c * lambda1 * lambda2 ^ 2 * sigma * xi ^ 2) - (2 * c * lambda1 * s2 * sigma * xi ^ 2) - (2 * c * lambda2 * s2 * sigma * xi ^ 2) + (6 * c * lambda1 ^ 2 * s2 * xi ^ 3) - (6 * c * lambda1 * sigma * xi) - (72 * c * lambda1 * mu1) + (6 * c * lambda1 ^ 2 * s2 * xi ^ 2) + (3 * c ^ 2 * lambda1 ^ 2 * s2 * xi ^ 2) - (72 * c * lambda1 * mu1 * xi) + (72 * c * lambda1 * mu2 * xi) + (6 * c ^ 2 * lambda1 * lambda2 * s2 * xi ^ 3) - (12 * c * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 3) + (72 * c * lambda1 ^ 2 * mu2 * xi) + (36 * c * lambda1 ^ 2) + (24 * c ^ 2 * lambda1 * lambda2 * xi ^ 2) + (24 * c * lambda1 * lambda2 ^ 2 * mu2 * xi ^ 2) - (12 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi) - (12 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi) - (72 * lambda1 * lambda2 * mu1 * xi ^ 2) + (48 * lambda1 * lambda2 * mu2 * xi ^ 2) - (12 * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 3) - (12 * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 3) - (36 * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 2) + (24 * lambda1 ^ 2 * lambda2 * mu2 * xi ^ 2) - (36 * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 2) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2) - (2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 2) + (24 * c * lambda1 * lambda2 * xi ^ 2) + (6 * c * lambda1 * s2 * xi ^ 2) + (6 * c * lambda2 * s2 * xi ^ 2) + (4 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2) + (2 * c * lambda1 ^ 2 * sigma * xi ^ 2) + (2 * c * lambda2 ^ 2 * sigma * xi ^ 2) + (6 * c * lambda1 ^ 2 * lambda2 * sigma) + (6 * c * lambda1 * lambda2 ^ 2 * sigma) - (2 * c * lambda1 * sigma * xi ^ 2) - (2 * c * lambda2 * sigma * xi ^ 2) + (6 * c * lambda1 ^ 2 * lambda2 * s2 * xi ^ 2) + (36 * c ^ 2 * lambda1 * lambda2 ^ 2 * xi) - (12 * c * lambda1 * mu1 * xi ^ 3) + (3 * c ^ 2 * lambda1 ^ 2 * lambda2 * s2 * xi ^ 2) + (12 * c * lambda2 ^ 2 * xi ^ 2) + (36 * c * lambda2 ^ 2 * xi) + (72 * c * lambda1 * lambda2 ^ 2 * mu2 * xi) - (72 * c * lambda1 * lambda2 * mu1 * xi ^ 2) + (12 * c * lambda1 * lambda2 * mu2 * s2 * xi ^ 2) - (12 * c * lambda1 ^ 2 * mu1 * xi ^ 3) + (12 * c * lambda1 * xi ^ 2) - (36 * c * lambda1 ^ 2 * mu1 * xi ^ 2) + (6 * c * lambda2 ^ 2 * mu2 * s2 * xi ^ 3) + (6 * c * lambda2 ^ 2 * mu2 * s2 * xi ^ 2) + (6 * c * lambda1 * mu2 * s2 * xi ^ 3) + (6 * c * lambda2 * mu2 * s2 * xi ^ 3) - (72 * c * lambda1 ^ 2 * mu1 * xi) + (3 * c ^ 2 * lambda1 * lambda2 ^ 2 * s2 * xi ^ 3) + (6 * c * lambda1 * lambda2 ^ 2 * s2 * xi ^ 2) + (12 * c * lambda1 * lambda2 * s2 * xi ^ 2) + (3 * c ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3) + (12 * c * lambda1 * lambda2 * mu2 * s2 * xi ^ 3) + (6 * c * lambda1 ^ 2 * lambda2 * mu2 * s2 * xi ^ 2) + (72 * c * lambda1 ^ 2 * lambda2 * mu2 * xi) - (72 * c * lambda1 * lambda2 ^ 2 * mu1 * xi) + (144 * c * lambda1 * lambda2 * mu2 * xi) - (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1) + (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2) - (12 * lambda1 ^ 2 * mu2 * sigma) - (12 * lambda2 ^ 2 * mu2 * sigma) + (24 * lambda1 * mu1 * sigma) - (36 * lambda1 * mu2 * sigma) + (24 * lambda2 * mu1 * sigma) - (36 * lambda2 * mu2 * sigma) + (12 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi) + (72 * lambda1 ^ 2 * lambda2 * mu2) + (72 * lambda1 * lambda2 ^ 2 * mu2) + (144 * lambda1 * lambda2 * mu2) + (12 * c * lambda1 ^ 2 * xi ^ 2) + (36 * c * lambda1 ^ 2 * xi) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3 + 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3 - c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2 + 3 * c ^ 2 * lambda1 ^ 2 * lambda2 * s2 * xi ^ 3 + 3 * c ^ 2 * lambda1 * lambda2 ^ 2 * s2 * xi ^ 3 - 4 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3 - 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2 + 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3 + 6 * c * lambda1 ^ 2 * lambda2 * mu2 * s2 * xi ^ 3 + 6 * c * lambda1 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3 - 4 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2 - 3 * c ^ 2 * lambda1 ^ 2 * lambda2 * s2 * xi ^ 2 + 3 * c ^ 2 * lambda1 ^ 2 * s2 * xi ^ 3 - 3 * c ^ 2 * lambda1 * lambda2 ^ 2 * s2 * xi ^ 2 + 6 * c ^ 2 * lambda1 * lambda2 * s2 * xi ^ 3 + 3 * c ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3 + 12 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2 - 8 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2 - 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2 - 12 * c * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 3 - 6 * c * lambda1 ^ 2 * lambda2 * mu2 * s2 * xi ^ 2 + 6 * c * lambda1 ^ 2 * lambda2 * s2 * xi ^ 3 + 6 * c * lambda1 ^ 2 * mu2 * s2 * xi ^ 3 - 12 * c * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 3 - 6 * c * lambda1 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2 + 6 * c * lambda1 * lambda2 ^ 2 * s2 * xi ^ 3 + 12 * c * lambda1 * lambda2 * mu2 * s2 * xi ^ 3 - 2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 3 + 6 * c * lambda2 ^ 2 * mu2 * s2 * xi ^ 3 - 4 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3 + 12 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi - 12 * c ^ 2 * lambda1 ^ 2 * lambda2 * xi ^ 2 - 3 * c ^ 2 * lambda1 ^ 2 * s2 * xi ^ 2 - 12 * c ^ 2 * lambda1 * lambda2 ^ 2 * xi ^ 2 - 6 * c ^ 2 * lambda1 * lambda2 * s2 * xi ^ 2 + 3 * c ^ 2 * lambda1 * s2 * xi ^ 3 - 3 * c ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2 + 3 * c ^ 2 * lambda2 * s2 * xi ^ 3 - 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi + 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi - 4 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2 + 36 * c * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 2 - 24 * c * lambda1 ^ 2 * lambda2 * mu2 * xi ^ 2 - 6 * c * lambda1 ^ 2 * lambda2 * s2 * xi ^ 2 - 2 * c * lambda1 ^ 2 * lambda2 * sigma * xi ^ 2 - 12 * c * lambda1 ^ 2 * mu1 * xi ^ 3 - 6 * c * lambda1 ^ 2 * mu2 * s2 * xi ^ 2 + 6 * c * lambda1 ^ 2 * s2 * xi ^ 3 + 36 * c * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 2 - 24 * c * lambda1 * lambda2 ^ 2 * mu2 * xi ^ 2 - 6 * c * lambda1 * lambda2 ^ 2 * s2 * xi ^ 2 - 2 * c * lambda1 * lambda2 ^ 2 * sigma * xi ^ 2 - 24 * c * lambda1 * lambda2 * mu1 * xi ^ 3 - 12 * c * lambda1 * lambda2 * mu2 * s2 * xi ^ 2 + 2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 2 + 12 * c * lambda1 * lambda2 * s2 * xi ^ 3 + 6 * c * lambda1 * mu2 * s2 * xi ^ 3 - 2 * c * lambda1 * s2 * sigma * xi ^ 3 - 12 * c * lambda2 ^ 2 * mu1 * xi ^ 3 - 6 * c * lambda2 ^ 2 * mu2 * s2 * xi ^ 2 + 6 * c * lambda2 ^ 2 * s2 * xi ^ 3 + 6 * c * lambda2 * mu2 * s2 * xi ^ 3 - 2 * c * lambda2 * s2 * sigma * xi ^ 3 + 12 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2 - 8 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2 - 12 * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 3 + 4 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi ^ 2 - 12 * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 3 + 4 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2 + 4 * lambda1 * lambda2 * mu1 * sigma * xi ^ 3 - 12 * lambda1 ^ 2 * c ^ 2 * lambda2 ^ 2 + 36 * c ^ 2 * lambda1 ^ 2 * lambda2 * xi - 12 * c ^ 2 * lambda1 ^ 2 * xi ^ 2 + 36 * c ^ 2 * lambda1 * lambda2 ^ 2 * xi - 24 * c ^ 2 * lambda1 * lambda2 * xi ^ 2 - 3 * c ^ 2 * lambda1 * s2 * xi ^ 2 - 12 * c ^ 2 * lambda2 ^ 2 * xi ^ 2 - 3 * c ^ 2 * lambda2 * s2 * xi ^ 2 + 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 - 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 + 12 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi - 72 * c * lambda1 ^ 2 * lambda2 * mu1 * xi + 72 * c * lambda1 ^ 2 * lambda2 * mu2 * xi + 6 * c * lambda1 ^ 2 * lambda2 * sigma * xi - 12 * c * lambda1 ^ 2 * lambda2 * xi ^ 2 + 36 * c * lambda1 ^ 2 * mu1 * xi ^ 2 - 24 * c * lambda1 ^ 2 * mu2 * xi ^ 2 - 6 * c * lambda1 ^ 2 * s2 * xi ^ 2 - 2 * c * lambda1 ^ 2 * sigma * xi ^ 2 - 72 * c * lambda1 * lambda2 ^ 2 * mu1 * xi + 72 * c * lambda1 * lambda2 ^ 2 * mu2 * xi + 6 * c * lambda1 * lambda2 ^ 2 * sigma * xi - 12 * c * lambda1 * lambda2 ^ 2 * xi ^ 2 + 72 * c * lambda1 * lambda2 * mu1 * xi ^ 2 - 48 * c * lambda1 * lambda2 * mu2 * xi ^ 2 - 12 * c * lambda1 * lambda2 * s2 * xi ^ 2 - 12 * c * lambda1 * mu1 * xi ^ 3 - 6 * c * lambda1 * mu2 * s2 * xi ^ 2 + 2 * c * lambda1 * s2 * sigma * xi ^ 2 + 6 * c * lambda1 * s2 * xi ^ 3 + 36 * c * lambda2 ^ 2 * mu1 * xi ^ 2 - 24 * c * lambda2 ^ 2 * mu2 * xi ^ 2 - 6 * c * lambda2 ^ 2 * s2 * xi ^ 2 - 2 * c * lambda2 ^ 2 * sigma * xi ^ 2 - 12 * c * lambda2 * mu1 * xi ^ 3 - 6 * c * lambda2 * mu2 * s2 * xi ^ 2 + 2 * c * lambda2 * s2 * sigma * xi ^ 2 + 6 * c * lambda2 * s2 * xi ^ 3 - 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi + 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi + 36 * lambda1 ^ 2 * lambda2 * mu1 * xi ^ 2 - 12 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi - 24 * lambda1 ^ 2 * lambda2 * mu2 * xi ^ 2 - 12 * lambda1 ^ 2 * xi ^ 3 * mu1 + 4 * lambda1 ^ 2 * mu2 * sigma * xi ^ 2 + 36 * lambda1 * lambda2 ^ 2 * mu1 * xi ^ 2 - 12 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi - 24 * lambda1 * lambda2 ^ 2 * mu2 * xi ^ 2 - 12 * lambda1 * lambda2 * mu1 * sigma * xi ^ 2 - 24 * lambda1 * lambda2 * mu1 * xi ^ 3 + 16 * lambda1 * lambda2 * mu2 * sigma * xi ^ 2 + 4 * lambda1 * mu1 * sigma * xi ^ 3 - 12 * lambda2 ^ 2 * mu1 * xi ^ 3 + 4 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2 + 4 * lambda2 * mu1 * sigma * xi ^ 3 - 36 * c ^ 2 * lambda1 ^ 2 * lambda2 + 36 * c ^ 2 * lambda1 ^ 2 * xi - 36 * c ^ 2 * lambda1 * lambda2 ^ 2 + 72 * c ^ 2 * lambda1 * lambda2 * xi - 12 * c ^ 2 * lambda1 * xi ^ 2 + 36 * c ^ 2 * lambda2 ^ 2 * xi - 12 * c ^ 2 * lambda2 * xi ^ 2 - 12 * lambda1 ^ 2 * c * lambda2 ^ 2 + 72 * c * lambda1 ^ 2 * lambda2 * mu1 - 72 * c * lambda1 ^ 2 * lambda2 * mu2 - 6 * c * lambda1 ^ 2 * lambda2 * sigma + 36 * c * lambda1 ^ 2 * lambda2 * xi - 72 * c * lambda1 ^ 2 * mu1 * xi + 72 * c * lambda1 ^ 2 * mu2 * xi + 6 * c * lambda1 ^ 2 * sigma * xi - 12 * c * lambda1 ^ 2 * xi ^ 2 + 72 * c * lambda1 * lambda2 ^ 2 * mu1 - 72 * c * lambda1 * lambda2 ^ 2 * mu2 - 6 * c * lambda1 * lambda2 ^ 2 * sigma + 36 * c * lambda1 * lambda2 ^ 2 * xi - 144 * c * lambda1 * lambda2 * mu1 * xi + 144 * c * lambda1 * lambda2 * mu2 * xi - 24 * c * lambda1 * lambda2 * xi ^ 2 + 36 * c * lambda1 * mu1 * xi ^ 2 - 24 * c * lambda1 * mu2 * xi ^ 2 - 6 * c * lambda1 * s2 * xi ^ 2 + 2 * c * lambda1 * sigma * xi ^ 2 - 72 * c * lambda2 ^ 2 * mu1 * xi + 72 * c * lambda2 ^ 2 * mu2 * xi + 6 * c * lambda2 ^ 2 * sigma * xi - 12 * c * lambda2 ^ 2 * xi ^ 2 + 36 * c * lambda2 * mu1 * xi ^ 2 - 24 * c * lambda2 * mu2 * xi ^ 2 - 6 * c * lambda2 * s2 * xi ^ 2 + 2 * c * lambda2 * sigma * xi ^ 2 + 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 - 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 - 72 * lambda1 ^ 2 * lambda2 * mu1 * xi + 12 * lambda1 ^ 2 * lambda2 * mu2 * sigma + 72 * lambda1 ^ 2 * lambda2 * mu2 * xi + 36 * lambda1 ^ 2 * xi ^ 2 * mu1 - 12 * lambda1 ^ 2 * mu2 * sigma * xi - 24 * lambda1 ^ 2 * mu2 * xi ^ 2 - 72 * lambda1 * lambda2 ^ 2 * mu1 * xi + 12 * lambda1 * lambda2 ^ 2 * mu2 * sigma + 72 * lambda1 * lambda2 ^ 2 * mu2 * xi + 24 * lambda1 * lambda2 * mu1 * sigma * xi + 72 * lambda1 * lambda2 * mu1 * xi ^ 2 - 48 * lambda1 * lambda2 * mu2 * sigma * xi - 48 * lambda1 * lambda2 * mu2 * xi ^ 2 - 12 * lambda1 * mu1 * sigma * xi ^ 2 - 12 * lambda1 * mu1 * xi ^ 3 + 12 * lambda1 * mu2 * sigma * xi ^ 2 + 36 * lambda2 ^ 2 * mu1 * xi ^ 2 - 12 * lambda2 ^ 2 * mu2 * sigma * xi - 24 * lambda2 ^ 2 * mu2 * xi ^ 2 - 12 * lambda2 * mu1 * sigma * xi ^ 2 - 12 * lambda2 * mu1 * xi ^ 3 + 12 * lambda2 * mu2 * sigma * xi ^ 2 - 36 * c ^ 2 * lambda1 ^ 2 - 72 * c ^ 2 * lambda1 * lambda2 + 36 * c ^ 2 * lambda1 * xi - 36 * c ^ 2 * lambda2 ^ 2 + 36 * c ^ 2 * lambda2 * xi - 36 * c * lambda1 ^ 2 * lambda2 + 72 * c * lambda1 ^ 2 * mu1 - 72 * c * lambda1 ^ 2 * mu2 - 6 * c * lambda1 ^ 2 * sigma + 36 * c * lambda1 ^ 2 * xi - 36 * c * lambda1 * lambda2 ^ 2 + 144 * c * lambda1 * lambda2 * mu1 - 144 * c * lambda1 * lambda2 * mu2 + 72 * c * lambda1 * lambda2 * xi - 72 * c * lambda1 * mu1 * xi + 72 * c * lambda1 * mu2 * xi - 6 * c * lambda1 * sigma * xi - 12 * c * lambda1 * xi ^ 2 + 72 * c * lambda2 ^ 2 * mu1 - 72 * c * lambda2 ^ 2 * mu2 - 6 * c * lambda2 ^ 2 * sigma + 36 * c * lambda2 ^ 2 * xi - 72 * c * lambda2 * mu1 * xi + 72 * c * lambda2 * mu2 * xi - 6 * c * lambda2 * sigma * xi - 12 * c * lambda2 * xi ^ 2 + 72 * lambda1 ^ 2 * lambda2 * mu1 - 72 * lambda1 ^ 2 * lambda2 * mu2 - 72 * lambda1 ^ 2 * mu1 * xi + 12 * lambda1 ^ 2 * mu2 * sigma + 72 * lambda1 ^ 2 * mu2 * xi + 72 * lambda1 * lambda2 ^ 2 * mu1 - 72 * lambda1 * lambda2 ^ 2 * mu2 - 24 * lambda1 * lambda2 * mu1 * sigma - 144 * lambda1 * lambda2 * mu1 * xi + 48 * lambda1 * lambda2 * mu2 * sigma + 144 * lambda1 * lambda2 * mu2 * xi + 24 * lambda1 * mu1 * sigma * xi + 36 * lambda1 * mu1 * xi ^ 2 - 36 * lambda1 * mu2 * sigma * xi - 24 * lambda1 * mu2 * xi ^ 2 - 72 * lambda2 ^ 2 * mu1 * xi + 12 * lambda2 ^ 2 * mu2 * sigma + 72 * lambda2 ^ 2 * mu2 * xi + 24 * lambda2 * mu1 * sigma * xi + 36 * lambda2 * mu1 * xi ^ 2 - 36 * lambda2 * mu2 * sigma * xi - 24 * lambda2 * mu2 * xi ^ 2 - 36 * c ^ 2 * lambda1 - 36 * c ^ 2 * lambda2 - 36 * c * lambda1 ^ 2 - 72 * c * lambda1 * lambda2 + 72 * c * lambda1 * mu1 - 72 * c * lambda1 * mu2 + 6 * c * lambda1 * sigma + 36 * c * lambda1 * xi - 36 * c * lambda2 ^ 2 + 72 * c * lambda2 * mu1 - 72 * c * lambda2 * mu2 + 6 * c * lambda2 * sigma + 36 * c * lambda2 * xi + 72 * lambda1 ^ 2 * mu1 - 72 * lambda1 ^ 2 * mu2 + 144 * lambda1 * lambda2 * mu1 - 144 * lambda1 * lambda2 * mu2 - 24 * lambda1 * mu1 * sigma - 72 * lambda1 * mu1 * xi + 36 * lambda1 * mu2 * sigma + 72 * lambda1 * mu2 * xi + 72 * lambda2 ^ 2 * mu1 - 72 * lambda2 ^ 2 * mu2 - 24 * lambda2 * mu1 * sigma - 72 * lambda2 * mu1 * xi + 36 * lambda2 * mu2 * sigma + 72 * lambda2 * mu2 * xi - 36 * c * lambda1 - 36 * c * lambda2 + 72 * lambda1 * mu1 - 72 * lambda1 * mu2 + 72 * lambda2 * mu1 - 72 * lambda2 * mu2) * exp((2 * xi))) / (-(12 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2) + (8 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2) - (4 * lambda1 ^ 2 * mu2 * sigma * xi ^ 2) - (4 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2) - (12 * lambda1 ^ 2 * lambda2 * mu2 * sigma) - (12 * lambda1 * lambda2 ^ 2 * mu2 * sigma) + (12 * lambda1 * mu1 * sigma * xi ^ 2) - (12 * lambda1 * mu2 * sigma * xi ^ 2) + (12 * lambda2 * mu1 * sigma * xi ^ 2) - (12 * lambda2 * mu2 * sigma * xi ^ 2) + (24 * lambda1 * lambda2 * mu1 * sigma) - (48 * lambda1 * lambda2 * mu2 * sigma) + (4 * lambda1 * lambda2 * mu1 * sigma * xi ^ 3) + (6 * c * lambda1 ^ 2 * sigma) + (6 * c * lambda2 ^ 2 * sigma) - (6 * c * lambda1 * sigma) - (6 * c * lambda2 * sigma) - (4 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3) - (4 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi ^ 2) - (4 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2) + (12 * lambda1 * lambda2 * mu1 * sigma * xi ^ 2) - (16 * lambda1 * lambda2 * mu2 * sigma * xi ^ 2) - (4 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3) + (24 * lambda1 * lambda2 * mu1 * sigma * xi) - (48 * lambda1 * lambda2 * mu2 * sigma * xi) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3 + 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3 - c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2 - 4 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3 - 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2 + 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3 - 4 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2 + 12 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2 - 8 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2 - 2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2 - 2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 3 - 4 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 3 + 12 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi - 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi + 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi - 4 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2 - 2 * c * lambda1 ^ 2 * lambda2 * sigma * xi ^ 2 - 2 * c * lambda1 * lambda2 ^ 2 * sigma * xi ^ 2 + 2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 2 - 2 * c * lambda1 * s2 * sigma * xi ^ 3 - 2 * c * lambda2 * s2 * sigma * xi ^ 3 + 12 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2 - 8 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2 + 4 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi ^ 2 + 4 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2 + 4 * lambda1 * lambda2 * mu1 * sigma * xi ^ 3 - 12 * lambda1 ^ 2 * c ^ 2 * lambda2 ^ 2 + 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 - 24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 + 12 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi + 6 * c * lambda1 ^ 2 * lambda2 * sigma * xi - 2 * c * lambda1 ^ 2 * sigma * xi ^ 2 + 6 * c * lambda1 * lambda2 ^ 2 * sigma * xi + 2 * c * lambda1 * s2 * sigma * xi ^ 2 - 2 * c * lambda2 ^ 2 * sigma * xi ^ 2 + 2 * c * lambda2 * s2 * sigma * xi ^ 2 - 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi + 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi - 12 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi + 4 * lambda1 ^ 2 * mu2 * sigma * xi ^ 2 - 12 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi - 12 * lambda1 * lambda2 * mu1 * sigma * xi ^ 2 + 16 * lambda1 * lambda2 * mu2 * sigma * xi ^ 2 + 4 * lambda1 * mu1 * sigma * xi ^ 3 + 4 * lambda2 ^ 2 * mu2 * sigma * xi ^ 2 + 4 * lambda2 * mu1 * sigma * xi ^ 3 - 12 * lambda1 ^ 2 * c * lambda2 ^ 2 - 6 * c * lambda1 ^ 2 * lambda2 * sigma + 6 * c * lambda1 ^ 2 * sigma * xi - 6 * c * lambda1 * lambda2 ^ 2 * sigma + 2 * c * lambda1 * sigma * xi ^ 2 + 6 * c * lambda2 ^ 2 * sigma * xi + 2 * c * lambda2 * sigma * xi ^ 2 + 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 - 24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 + 12 * lambda1 ^ 2 * lambda2 * mu2 * sigma - 12 * lambda1 ^ 2 * mu2 * sigma * xi + 12 * lambda1 * lambda2 ^ 2 * mu2 * sigma + 24 * lambda1 * lambda2 * mu1 * sigma * xi - 48 * lambda1 * lambda2 * mu2 * sigma * xi - 12 * lambda1 * mu1 * sigma * xi ^ 2 + 12 * lambda1 * mu2 * sigma * xi ^ 2 - 12 * lambda2 ^ 2 * mu2 * sigma * xi - 12 * lambda2 * mu1 * sigma * xi ^ 2 + 12 * lambda2 * mu2 * sigma * xi ^ 2 - 6 * c * lambda1 ^ 2 * sigma - 6 * c * lambda1 * sigma * xi - 6 * c * lambda2 ^ 2 * sigma - 6 * c * lambda2 * sigma * xi + 12 * lambda1 ^ 2 * mu2 * sigma - 24 * lambda1 * lambda2 * mu1 * sigma + 48 * lambda1 * lambda2 * mu2 * sigma + 24 * lambda1 * mu1 * sigma * xi - 36 * lambda1 * mu2 * sigma * xi + 12 * lambda2 ^ 2 * mu2 * sigma + 24 * lambda2 * mu1 * sigma * xi - 36 * lambda2 * mu2 * sigma * xi + 6 * c * lambda1 * sigma + 6 * c * lambda2 * sigma - 24 * lambda1 * mu1 * sigma + 36 * lambda1 * mu2 * sigma - 24 * lambda2 * mu1 * sigma + 36 * lambda2 * mu2 * sigma) * exp((2 * xi)) - (36 * lambda2 * mu2 * sigma * xi) - (36 * lambda1 * mu2 * sigma * xi) + (24 * lambda2 * mu1 * sigma * xi) + (24 * lambda1 * mu1 * sigma * xi) + (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi) + (4 * lambda1 * mu1 * sigma * xi ^ 3) + (4 * lambda2 * mu1 * sigma * xi ^ 3) - (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3) - (2 * c * lambda2 * s2 * sigma * xi ^ 3) + (6 * c * lambda1 ^ 2 * lambda2 * sigma * xi) + (6 * c * lambda1 * lambda2 ^ 2 * sigma * xi) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 3) - (2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 3) - (12 * lambda1 ^ 2 * mu2 * sigma * xi) - (12 * lambda2 ^ 2 * mu2 * sigma * xi) - (2 * c * lambda1 * s2 * sigma * xi ^ 3) + (12 * lambda1 ^ 2 * c ^ 2 * lambda2 ^ 2) + (12 * lambda1 ^ 2 * c * lambda2 ^ 2) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 3) + (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi) - (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi) + (c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2) - (12 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1 * xi ^ 2) + (8 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * xi ^ 2) + (12 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi) + (6 * c * lambda1 ^ 2 * sigma * xi) + (6 * c * lambda2 ^ 2 * sigma * xi) - (6 * c * lambda2 * sigma * xi) + (4 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2) - (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu1) + (24 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * s2 * xi ^ 2) + (2 * c * lambda1 ^ 2 * lambda2 * sigma * xi ^ 2) + (2 * c * lambda1 * lambda2 ^ 2 * sigma * xi ^ 2) - (2 * c * lambda1 * s2 * sigma * xi ^ 2) - (2 * c * lambda2 * s2 * sigma * xi ^ 2) - (6 * c * lambda1 * sigma * xi) - (12 * lambda1 ^ 2 * lambda2 * mu2 * sigma * xi) - (12 * lambda1 * lambda2 ^ 2 * mu2 * sigma * xi) + (2 * c * lambda1 ^ 2 * lambda2 ^ 2 * mu2 * s2 * xi ^ 2) - (2 * c * lambda1 * lambda2 * s2 * sigma * xi ^ 2) + (4 * c * lambda1 ^ 2 * lambda2 ^ 2 * xi ^ 2) + (2 * c * lambda1 ^ 2 * sigma * xi ^ 2) + (2 * c * lambda2 ^ 2 * sigma * xi ^ 2) + (6 * c * lambda1 ^ 2 * lambda2 * sigma) + (6 * c * lambda1 * lambda2 ^ 2 * sigma) - (2 * c * lambda1 * sigma * xi ^ 2) - (2 * c * lambda2 * sigma * xi ^ 2) - (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu1) + (24 * lambda1 ^ 2 * lambda2 ^ 2 * mu2) - (12 * lambda1 ^ 2 * mu2 * sigma) - (12 * lambda2 ^ 2 * mu2 * sigma) + (24 * lambda1 * mu1 * sigma) - (36 * lambda1 * mu2 * sigma) + (24 * lambda2 * mu1 * sigma) - (36 * lambda2 * mu2 * sigma) + (12 * c ^ 2 * lambda1 ^ 2 * lambda2 ^ 2 * xi)) / 0.2e1;
x0=-(1/2)*(alpha^3*c^2*l^2*s2-alpha^3*c*l^3*mu1+3*alpha^3*c^2*l*s2-2*alpha^3*c*l^2*mu1-6*alpha*c*l^3*mu1-12*alpha*c*l^2*mu1+6*alpha*c*l^2*mu2+3*alpha^2*c^2*l*s2-3*alpha^2*l^3*mu1*mu2-3*alpha^2*l^2*mu1*mu2-6*l^3*mu1*mu2-6*l^2*mu1*mu2+9*c^2-alpha^3*l^3*mu1*mu2-alpha^3*l^2*mu1*mu2-6*alpha*l^3*mu1*mu2-6*alpha*l^2*mu1*mu2+3*alpha^2*c^2+9*c^2*l-18*c*mu1+9*alpha*c^2+3*c^2*l^2+alpha^2*c^2*l^2-6*c*l^3*mu1-12*c*l^2*mu1+6*c*l^2*mu2+3*alpha*c^2*l^2+6*alpha^2*c*l*mu2*s2-9*alpha^2*c*l*mu1+6*alpha^2*c*l*mu2+6*alpha^2*c*mu2*s2-3*alpha^3*c*l*mu1+6*alpha^3*c*mu2*s2-18*alpha*c*l*mu1+18*alpha*c*l*mu2-18*alpha*c*mu1+18*alpha*c*mu2+18*mu2*c+3*alpha^2*c^2*l+3*alpha^2*c^2*s2-9*alpha^2*c*mu1+6*alpha^2*c*mu2-18*c*l*mu1+18*c*l*mu2+3*alpha^3*c^2*s2-3*alpha^3*c*mu1+9*alpha*c^2*l+(alpha^3*c^2*l^2*s2-alpha^3*c*l^3*mu1+2*alpha^3*c*l^2*mu2*s2-alpha^3*l^3*mu1*mu2+3*alpha^3*c^2*l*s2-2*alpha^3*c*l^2*mu1+6*alpha^3*c*l*mu2*s2-alpha^3*l^2*mu1*mu2-alpha^2*c^2*l^2*s2+3*alpha^2*c*l^3*mu1-2*alpha^2*c*l^2*mu2*s2+3*alpha^2*l^3*mu1*mu2+3*alpha^3*c^2*s2-3*alpha^3*c*l*mu1+6*alpha^3*c*mu2*s2-alpha^2*c^2*l^2-3*alpha^2*c^2*l*s2+6*alpha^2*c*l^2*mu1-2*alpha^2*c*l^2*mu2-6*alpha^2*c*l*mu2*s2+3*alpha^2*l^2*mu1*mu2-6*alpha*c*l^3*mu1-6*alpha*l^3*mu1*mu2-3*alpha^3*c*mu1-3*alpha^2*c^2*l-3*alpha^2*c^2*s2+9*alpha^2*c*l*mu1-6*alpha^2*c*l*mu2-6*alpha^2*c*mu2*s2+3*alpha*c^2*l^2-12*alpha*c*l^2*mu1+6*alpha*c*l^2*mu2-6*alpha*l^2*mu1*mu2+6*c*l^3*mu1+6*l^3*mu1*mu2-3*alpha^2*c^2+9*alpha^2*c*mu1-6*alpha^2*c*mu2+9*alpha*c^2*l-18*alpha*c*l*mu1+18*alpha*c*l*mu2-3*c^2*l^2+12*c*l^2*mu1-6*c*l^2*mu2+6*l^2*mu1*mu2+9*alpha*c^2-18*alpha*c*mu1+18*alpha*c*mu2-9*c^2*l+18*c*l*mu1-18*c*l*mu2-9*c^2+18*c*mu1-18*c*mu2)*exp(2*alpha)+2*alpha^2*c*l^2*mu2*s2+2*alpha^3*c*l^2*mu2*s2+6*alpha^3*c*l*mu2*s2+alpha^2*c^2*l^2*s2-3*alpha^2*c*l^3*mu1-6*alpha^2*c*l^2*mu1+2*alpha^2*c*l^2*mu2)/(l^2*((alpha^3*c^2*s2+2*alpha^3*c*l*mu1+2*alpha^3*c*mu2*s2+2*alpha^3*l*mu1*mu2+alpha^3*c*mu1+2*alpha^3*mu1*mu2+3*alpha^2*c^2*l-alpha^2*c^2*s2-6*alpha^2*c*l*mu1+9*alpha^2*c*l*mu2-2*alpha^2*c*mu2*s2-6*alpha^2*l*mu1*mu2+6*alpha^2*l*mu2^2+2*alpha^2*c^2-3*alpha^2*c*mu1+7*alpha^2*c*mu2-6*alpha^2*mu1*mu2+6*alpha^2*mu2^2-9*alpha*c^2*l+12*alpha*c*l*mu1-27*alpha*c*l*mu2+12*alpha*l*mu1*mu2-18*alpha*l*mu2^2-6*alpha*c^2+6*alpha*c*mu1-21*alpha*c*mu2+12*alpha*mu1*mu2-18*alpha*mu2^2+9*c^2*l-12*c*l*mu1+27*c*l*mu2-12*l*mu1*mu2+18*l*mu2^2+6*c^2-6*c*mu1+21*c*mu2-12*mu1*mu2+18*mu2^2)*exp(2*alpha)-6*alpha^2*l*mu2^2+6*alpha^2*mu1*mu2+12*l*mu1*mu2-6*alpha^2*mu2^2-18*l*mu2^2+12*mu1*mu2-6*c^2+2*alpha^3*l*mu1*mu2+12*alpha*l*mu1*mu2-2*alpha^2*c^2-9*c^2*l+6*c*mu1-6*alpha*c^2-18*alpha*mu2^2+2*alpha^3*mu1*mu2-18*alpha*l*mu2^2+12*alpha*mu1*mu2+6*alpha^2*l*mu1*mu2+6*alpha^2*c*l*mu1-9*alpha^2*c*l*mu2+2*alpha^2*c*mu2*s2+2*alpha^3*c*l*mu1+2*alpha^3*c*mu2*s2+12*alpha*c*l*mu1-27*alpha*c*l*mu2+6*alpha*c*mu1-21*alpha*c*mu2-21*mu2*c-18*mu2^2-3*alpha^2*c^2*l+alpha^2*c^2*s2+3*alpha^2*c*mu1-7*alpha^2*c*mu2+12*c*l*mu1-27*c*l*mu2+alpha^3*c^2*s2+alpha^3*c*mu1-9*alpha*c^2*l));
val =(3/(4*pi*sigma^2*(rho1-1 -3*x1)))/(1/(2*(2+c)*pi*x0));
function ilt = talbot_inversion(F, t, M)
error('All time values must be positive.');
ilt = zeros(length(t), 1);
(2*M/5) * theta .* (cot_theta + 1i);
gamma(1) = 0.5 * exp(delta(1));
+ 1i*theta.*(1 + cot_theta.^2) ...
f_vals = arrayfun(F, s_vals);
(0.4/t(i)) * sum(real(gamma .* f_vals));
ilt = reshape(ilt, original_size);