clc; fC1=5000; fC1T=fC1/(2^0.5-1)^0.5; disp('kHz',fC1T/1000,"fC1T=");