Physics 360
Winter 2006

Matlab file: Click here to download.

Some examples of output.

Enter the total number of energy units: 6
Enter the number of oscillators in A:   3
Enter the number of oscillators in B:   3

qA OmegaA qB OmegaB Omega
 0    1    6   28     28
 1    3    5   21     63
 2    6    4   15     90
 3   10    3   10    100
 4   15    2    6     90
 5   21    1    3     63
 6   28    0    1     28

Total number of states = 462
The maximum value of Omega is 100 when qA = 3.

 

a

 

Enter the total number of energy units: 100
Enter the number of oscillators in A:   3
Enter the number of oscillators in B:   3

 qA OmegaA  qB OmegaB Omega

  0    1   100  5151   5151
  1    3    99  5050  15150
  2    6    98  4950  29700
  3   10    97  4851  48510
  4   15    96  4753  71295
  5   21    95  4656  97776
  6   28    94  4560  127680
  7   36    93  4465  160740
  8   45    92  4371  196695
  9   55    91  4278  235290
 10   66    90  4186  276276
 11   78    89  4095  319410
 12   91    88  4005  364455
 13  105    87  3916  411180
 14  120    86  3828  459360
 15  136    85  3741  508776
 16  153    84  3655  559215
 17  171    83  3570  610470
 18  190    82  3486  662340
 19  210    81  3403  714630
 20  231    80  3321  767151
 21  253    79  3240  819720
 22  276    78  3160  872160
 23  300    77  3081  924300
 24  325    76  3003  975975
 25  351    75  2926 1027026
 26  378    74  2850 1077300
 27  406    73  2775 1126650
 28  435    72  2701 1174935
 29  465    71  2628 1222020
 30  496    70  2556 1267776
 31  528    69  2485 1312080
 32  561    68  2415 1354815
 33  595    67  2346 1395870
 34  630    66  2278 1435140
 35  666    65  2211 1472526
 36  703    64  2145 1507935
 37  741    63  2080 1541280
 38  780    62  2016 1572480
 39  820    61  1953 1601460
 40  861    60  1891 1628151
 41  903    59  1830 1652490
 42  946    58  1770 1674420
 43  990    57  1711 1693890
 44 1035    56  1653 1710855
 45 1081    55  1596 1725276
 46 1128    54  1540 1737120
 47 1176    53  1485 1746360
 48 1225    52  1431 1752975
 49 1275    51  1378 1756950
 50 1326    50  1326 1758276
 51 1378    49  1275 1756950
 52 1431    48  1225 1752975
 53 1485    47  1176 1746360
 54 1540    46  1128 1737120
 55 1596    45  1081 1725276
 56 1653    44  1035 1710855
 57 1711    43   990 1693890
 58 1770    42   946 1674420
 59 1830    41   903 1652490
 60 1891    40   861 1628151
 61 1953    39   820 1601460
 62 2016    38   780 1572480
 63 2080    37   741 1541280
 64 2145    36   703 1507935
 65 2211    35   666 1472526
 66 2278    34   630 1435140
 67 2346    33   595 1395870
 68 2415    32   561 1354815
 69 2485    31   528 1312080
 70 2556    30   496 1267776
 71 2628    29   465 1222020
 72 2701    28   435 1174935
 73 2775    27   406 1126650
 74 2850    26   378 1077300
 75 2926    25   351 1027026
 76 3003    24   325  975975
 77 3081    23   300  924300
 78 3160    22   276  872160
 79 3240    21   253  819720
 80 3321    20   231  767151
 81 3403    19   210  714630
 82 3486    18   190  662340
 83 3570    17   171  610470
 84 3655    16   153  559215
 85 3741    15   136  508776
 86 3828    14   120  459360
 87 3916    13   105  411180
 88 4005    12    91  364455
 89 4095    11    78  319410
 90 4186    10    66  276276
 91 4278     9    55  235290
 92 4371     8    45  196695
 93 4465     7    36  160740
 94 4560     6    28  127680
 95 4656     5    21   97776
 96 4753     4    15   71295
 97 4851     3    10   48510
 98 4950     2     6   29700
 99 5050     1     3   15150
100 5151     0     1    5151

Total number of states = 96560646
The maximum value of Omega is 1758276 when qA = 50.

 

Enter the total number of energy units: 100
Enter the number of oscillators in A:   40
Enter the number of oscillators in B:   60

 qA        OmegaA  qB        OmegaB         Omega
  0             1 100 2.276717e+044 2.276717e+044
  1            40  99 1.431897e+044 5.727589e+045
  2           820  98 8.972015e+043 7.357052e+046
  3         11480  97 5.600366e+043 6.429220e+047
  4         12341  96 3.482279e+043 4.297480e+048
  5       1086008  95 2.156766e+043 2.342265e+049
  6       8145060  94 1.330473e+043 1.083678e+050
  7      53524680  93 8.174146e+042 4.375185e+050
  8 3.144575e+008  92 5.001287e+042 1.572692e+051
  9 1.677107e+009  91 3.047142e+042 5.110381e+051
 10 8.217823e+009  90 1.848599e+042 1.519146e+052
 11 3.735374e+010  89 1.116604e+042 4.170932e+052
 12 1.587534e+011  88 6.714711e+041 1.065983e+053
 13 6.350136e+011  87 4.019691e+041 2.552558e+053
 14 2.403980e+012  86 2.395295e+041 5.758241e+053
 15 8.654328e+012  85 1.420658e+041 1.229484e+054
 16 2.974925e+013  84 8.385827e+040 2.494721e+054
 17 9.799753e+013  83 4.925940e+040 4.827300e+054
 18 3.103255e+014  82 2.879247e+040 8.935038e+054
 19 9.473095e+014  81 1.674456e+040 1.586228e+055
 20 2.794563e+015  80 9.687922e+039 2.707351e+055
 21 7.984466e+015  79 5.575782e+039 4.451964e+055
 22 2.213875e+016  78 3.191933e+039 7.066540e+055
 23 5.967836e+016  77 1.817305e+039 1.084538e+056
 24 1.566557e+017  76 1.028915e+039 1.611855e+056
 25 4.010386e+017  75 5.792413e+038 2.322981e+056
 26 1.002596e+018  74 3.242022e+038 3.250440e+056
 27 2.450791e+018  73 1.803832e+038 4.420815e+056
 28 5.864393e+018  72 9.975737e+037 5.850164e+056
 29 1.375099e+019  71 5.482848e+037 7.539459e+056
 30 3.162728e+019  70 2.994478e+037 9.470720e+056
 31 7.141644e+019  69 1.624911e+037 1.160453e+057
 32 1.584552e+020  68 8.759284e+036 1.387954e+057
 33 3.457205e+020  67 4.690010e+036 1.621433e+057
 34 7.422822e+020  66 2.493894e+036 1.851174e+057
 35 1.569397e+021  65 1.316776e+036 2.066544e+057
 36 3.269576e+021  64 6.902456e+035 2.256811e+057
 37 6.715887e+021  63 3.591522e+035 2.412025e+057
 38 1.360851e+022  62 1.854638e+035 2.523886e+057
 39 2.721701e+022  61 9.503106e+034 2.586462e+057
 40 5.375360e+022  60 4.830745e+034 2.596700e+057
 41 1.048851e+023  59 2.435670e+034 2.554654e+057
 42 2.022784e+023  58 1.217835e+034 2.463417e+057
 43 3.857402e+023  57 6.037131e+033 2.328764e+057
 44 7.276462e+023  56 2.966521e+033 2.158578e+057
 45 1.358273e+024  55 1.444567e+033 1.962116e+057
 46 2.509852e+024  54 6.969401e+032 1.749217e+057
 47 4.592495e+024  53 3.330510e+032 1.529535e+057
 48 8.323898e+024  52 1.576045e+032 1.311884e+057
 49 1.494904e+025  51 7.383274e+031 1.103729e+057
 50 2.660929e+025  50 3.423154e+031 9.108772e+056
 51 4.695758e+025  49 1.570254e+031 7.373534e+056
 52 8.217576e+025  48 7.124302e+030 5.854449e+056
 53 1.426447e+026  47 3.195949e+030 4.558852e+056
 54 2.456659e+026  46 1.417072e+030 3.481261e+056
 55 4.198653e+026  45 6.208123e+029 2.606576e+056
 56 7.122716e+026  44 2.686207e+029 1.913309e+056
 57 1.199615e+027  43 1.147506e+029 1.376566e+056
 58 2.006253e+027  42 4.837525e+028 9.705299e+055
 59 3.332420e+027  41 2.011644e+028 6.703644e+055
 60 5.498494e+027  40 8.247740e+027 4.535015e+055
 61 9.013924e+027  39 3.332420e+027 3.003818e+055
 62 1.468397e+028  38 1.326167e+027 1.947340e+055
 63 2.377405e+028  37 5.195295e+026 1.235132e+055
 64 3.826136e+028  36 2.002353e+026 7.661276e+054
 65 6.121818e+028  35 7.587864e+025 4.645153e+054
 66 9.739256e+028  34 2.825269e+025 2.751602e+054
 67 1.540838e+029  33 1.032894e+025 1.591522e+054
 68 2.424553e+029  32 3.704946e+024 8.982838e+053
 69 3.794953e+029  31 1.302838e+024 4.944209e+053
 70 5.909284e+029  30 4.487553e+023 2.651823e+053
 71 9.155228e+029  29 1.512658e+023 1.384873e+053
 72 1.411431e+030  28 4.984897e+022 7.035838e+052
 73 2.165483e+030  27 1.604335e+022 3.474160e+052
 74 3.306751e+030  26 5.036865e+021 1.665566e+052
 75 5.026262e+030  25 1.540688e+021 7.743902e+051
 76 7.605528e+030  24 4.585381e+020 3.487424e+051
 77 1.145768e+031  23 1.325893e+020 1.519166e+051
 78 1.718652e+031  22 3.718969e+019 6.391613e+050
 79 2.567100e+031  21 1.010090e+019 2.593003e+050
 80 3.818562e+031  20 2.651487e+018 1.012487e+050
 81 5.657128e+031  19 6.712626e+017 3.797418e+049
 82 8.347714e+031  18 1.635127e+017 1.364957e+049
 83 1.227013e+032  17 3.822374e+016 4.690104e+048
 84 1.796698e+032  16 8.550048e+015 1.536185e+048
 85 2.621065e+032  15 1.824010e+015 4.780850e+047
 86 3.809688e+032  14 3.697318e+014 1.408563e+047
 87 5.517479e+032  13 7.090747e+013 3.912305e+046
 88 7.962726e+032  12 1.280274e+013 1.019447e+046
 89 1.145201e+033  11 2.163843e+012 2.478035e+045
 90 1.641455e+033  10 3.400324e+011 5.581479e+044
 91 2.344935e+033   9 4.928007e+010 1.155586e+044
 92 3.338984e+033   8 6.522362e+009 2.177806e+043
 93 4.739203e+033   7 7.787894e+008 3.690842e+042
 94 6.705469e+033   6 8.259888e+007 5.538642e+041
 95 9.458240e+033   5 7.624512e+006 7.211446e+040
 96 1.330065e+034   4 5.956650e+005 7.922732e+039
 97 1.864833e+034   3 3.782000e+004 7.052800e+038
 98 2.606961e+034   2 1.830000e+003 4.770739e+037
 99 3.633946e+034   1            60 2.180367e+036
100 5.051184e+034   0             1 5.051184e+034

Total number of states = 4.527426e+058
The maximum value of Omega is 2.596700e+057 when qA = 40.