Module[{cw, ca, cvxa, cvza, cvxw, cvzw, cxa, cxw, ua1, uw1, Q1, kc,
et, rn, a1, a2, a3, a4, a5, a6, a7, bt1, bt2, fi1, fi2, fi3, fi4,
fi5, fi6},
cw = (1 - m2w/m1kw)/(m2w/m1kw);
ca = m2a/(m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cvxw = kxw/(rw*m2w); cvzw = kzw/(rw*m2w);
cvxa = (kxa*R1*T1)/(
g*M*ua0*m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cvza = (kza*R1*T1)/(
g*M*ua0*m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cxa = m1ka/(m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cxw = m1kw/m2w;
ua1 = (4*ua0)/((2*k + 1)*\[Pi]);
uw1 = (4*uw0)/((2*k + 1)*\[Pi]);
Q1 = s*(LaplaceTransform[
Which[0 < Mod[t, 3456000] < 1728000, 1.16*10^-4*Mod[t, 3456000],
1728000 <= Mod[t, 3456000] <= 3456000,
200 - 1.16*10^-4*(Mod[t, 3456000] - 1728000)], t, s]) - q0;
a1 = (cvxa*cvxw)/(s*ca);
a2 = (cvxw + cvxa)/ca;
a3 = -s*cw + s/ca;
a4 = (ca*cw - 1)/ca*ua0 + (ca*cxw - cxa)/ca*Q1;
kc = Sqrt[(-a2 - Sqrt[a2*a2 - 4*a1*a3])/(2*a1)];
et = Sqrt[(-a2 + Sqrt[a2*a2 - 4*a1*a3])/(2*a1)];
a5 = -cvxa/(s*ca)*kc*kc - 1/ca;
a6 = (-cvxa*et*et)/(s*ca) - 1/ca;
a7 = a4/(a3*ca) + uw0/s + ua0/(s*ca) + cxa/(s*ca)*Q1;
(-a4*a6 - a3*a7)/(a3*(a5 - a6)*(1 + E^(2 kc*h)))*
E^(kc*z) + (-(a4*a6 + a3*a7)*E^(2*kc*h))/(
a3*(a5 - a6)*(1 + E^(2 kc*h)))*E^(-kc*z) + (a4*a5 + a3*a7)/(
a3*(a5 - a6)*(1 + E^(2 et*h)))*E^(et*z) + (
E^(2*et*h)*(a4*a5 + a3*a7))/(a3*(a5 - a6)*(1 + E^(2 et*h)))*
E^(-et*z) - a4/a3]
time = Table[10^j, {j, 0, 10, 5}];
ft = Durbin[P[1, s], s, time]
et, rn, a1, a2, a3, a4, a5, a6, a7, bt1, bt2, fi1, fi2, fi3, fi4,
fi5, fi6},
cw = (1 - m2w/m1kw)/(m2w/m1kw);
ca = m2a/(m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cvxw = kxw/(rw*m2w); cvzw = kzw/(rw*m2w);
cvxa = (kxa*R1*T1)/(
g*M*ua0*m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cvza = (kza*R1*T1)/(
g*M*ua0*m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cxa = m1ka/(m1ka - m2a - (uatm*n0*(1 - sr0))/(ua0*ua0));
cxw = m1kw/m2w;
ua1 = (4*ua0)/((2*k + 1)*\[Pi]);
uw1 = (4*uw0)/((2*k + 1)*\[Pi]);
Q1 = s*(LaplaceTransform[
Which[0 < Mod[t, 3456000] < 1728000, 1.16*10^-4*Mod[t, 3456000],
1728000 <= Mod[t, 3456000] <= 3456000,
200 - 1.16*10^-4*(Mod[t, 3456000] - 1728000)], t, s]) - q0;
a1 = (cvxa*cvxw)/(s*ca);
a2 = (cvxw + cvxa)/ca;
a3 = -s*cw + s/ca;
a4 = (ca*cw - 1)/ca*ua0 + (ca*cxw - cxa)/ca*Q1;
kc = Sqrt[(-a2 - Sqrt[a2*a2 - 4*a1*a3])/(2*a1)];
et = Sqrt[(-a2 + Sqrt[a2*a2 - 4*a1*a3])/(2*a1)];
a5 = -cvxa/(s*ca)*kc*kc - 1/ca;
a6 = (-cvxa*et*et)/(s*ca) - 1/ca;
a7 = a4/(a3*ca) + uw0/s + ua0/(s*ca) + cxa/(s*ca)*Q1;
(-a4*a6 - a3*a7)/(a3*(a5 - a6)*(1 + E^(2 kc*h)))*
E^(kc*z) + (-(a4*a6 + a3*a7)*E^(2*kc*h))/(
a3*(a5 - a6)*(1 + E^(2 kc*h)))*E^(-kc*z) + (a4*a5 + a3*a7)/(
a3*(a5 - a6)*(1 + E^(2 et*h)))*E^(et*z) + (
E^(2*et*h)*(a4*a5 + a3*a7))/(a3*(a5 - a6)*(1 + E^(2 et*h)))*
E^(-et*z) - a4/a3]
time = Table[10^j, {j, 0, 10, 5}];
ft = Durbin[P[1, s], s, time]