Ka = abs(simps(S1 * exp(1J * f * Tm02), f)) / m[0]