求住 谁能帮我用java或c语言 实现一下公式 任意输入数字可出结果~~

(1)L=10log{[Re(P_入/P_透)]²+[Im(P_入/P_透)]²}

(2)Re(P_入/P_透)=1/ 4m (1+cos2kl2){4mcosk(l1-l2)+( 3m -1)cosk(l1+l2)+(m+1)cosk(l1 -3l 2)

(3)Im(P_入/P_透)=1/ 4m (1+cos2kl2){2(m²+1)sink(l1-l2)+( 2m ² -m+1)sink(l1+l2)+(m+1)sink(l1 -3l 2)

式中P_入-------消声器的入射声压，      P_透---------消声器透射声压

m---------扩张比（m=消声器体径向截面积/插入管截面积）

l1---------消声器单室扩张室长度（米）；l2--------插入管长度（米）；

k=2∏f/C;      f--------消声频率（赫）；C----------速度（米）