clc; //e.g 13.6 L1=40*10**-6; L2=80*10**-6; M=11.3*10**-6; k=M/sqrt(L1*L2); disp(k);