> `b := sqrt(2*n+1) `;

                                b := sqrt(2*n+1)
> b^2 / 2 * hbar * w = 1 / 2 * m *w^2;

                                 2                 2
                            1/2 b  hbar w = 1/2 m w
> solve(",b^2);

                                       m w
                                      ----
                                      hbar
