On Mon, Aug 3, 2009 at 7:56 AM, Barton Willis<willisb at unk.edu> wrote:
> (%o3) [e=((sqrt(3)*%pi^(4/3)*%i-%pi^(4/3))*hbar^(4/3)*mu^(1/3)*(2*n+1)^
> (4/3))/(2*2^(2/3)*beta(1/4,3/2)^(4/3)*m^(2/3)),e=-((sqrt(3)*%pi^(4/3)*%i
> +%pi^(4/3))*hbar^(4/3)*mu^(1/3)*(2*n+1)^(4/3))/(2*2^(2/3)*beta(1/4,3/2)^
> (4/3)*m^(2/3)),e=(%pi^(4/3)*hbar^(4/3)*mu^(1/3)*(2*n+1)^(4/3))/(2^
> (2/3)*beta(1/4,3/2)^(4/3)*m^(2/3))]
FWIW the value of hbar is in a share package physical_constants
(as %h_bar). physical_constants makes use of the ezunits units stuff.
?? physical_constants
=> info
load (physical_constants);
propvars (physical_constant);
=> [%c, %mu_0, %e_0, %Z_0, %G, %h, %h_bar, %m_P, %T_P, %l_P, %t_P, %%e,
%Phi_0, %G_0, %K_J, %R_K, %mu_B, %mu_N, %alpha, %R_inf, %a_0, %E_h,
%ratio_h_me, %m_e, %N_A, %m_u, %F, %R, %%k, %V_m, %n_0, %ratio_S0_R, %sigma,
%c_1, %c_1L, %c_2, %b, %b_prime]
for x in propvars (physical_constant) do disp (x = get (x, 'description));
=> each constant & its description
constvalue (%h_bar);
=> 41412931/(125000000000000000000000000000000000000000*%pi) ` s*J
float (%);
=> 1.054571628251774E-34 ` s J
Be advised I fixed some bugs in this stuff recently (incorrect values, I think).
HTH in some way.
Robert Dodier