lda0 1[m] "波长，真空"
f0 c_const/lda0 频率
e_a 1 "相对介电常数，真空"
mu_a 1 "相对磁导率，真空"
e_b -1 "相对介电常数，超材料"
mu_b -1 "相对磁导率，超材料"
n_a sqrt(e_a*mu_a) "折射率，真空"
n_b -sqrt(e_b*mu_b) "折射率，超材料"
alpha 30[deg] "入射角"