// Scilab Code Ex15.1: Page-323 (2010)
k = 1.38e-023;    // Boltzmann constant, J/K
h = 6.626e-034;    // Planck's constant, Js
f_D = 64e+011;    // Debye frequency for Al, Hz
theta_D = h*f_D/k;    // Debye temperature, K
printf("\nThe Debye temperature of aluminium = %5.1f K", theta_D);

// Result
// The Debye temperature of aluminium = 307.3 K