Solve the equation without giving numerical values to the variables:
Reduce[(wh[n, B] == wh1[n, B]) && e0 > 0 && 0 < m && 0 < M && B > 0 &&
e > 0, Reals]
You will see that the only solution in the reals can be n == 0
. Still, the two functions wh
and wh1
are very close in numerical value:
values = {B -> 0.32, e -> 1.6*10^(-19),
e0 -> 8.8*10^(-12),
m -> 0.91*10^(-30),
M -> 1.67*10^(-27)};
Plot[(wh[n, B] - wh1[n, B]) /. values,
{n, 0, 1000000}, PlotRange -> .0001]