在R中编写一个函数来求解Plank方程

plank_equation <- function(frequency, temperature) {
  h <- 6.62607015e-34  # 普朗克常数
  c <- 2.99792458e8    # 光速
  k <- 1.380649e-23    # 玻尔兹曼常数
  
  energy_density <- (2 * h * frequency^3) / (c^2 * (exp((h * frequency) / (k * temperature)) - 1))
  
  return(energy_density)
}

frequency <- 1e12  # 频率为1 THz
temperature <- 300  # 温度为300K

result <- plank_equation(frequency, temperature)
print(result)