% w0/d0   phL^2/2D  ph/2(kD)^(1/2)  L(k/D)^(1/4)  delph/2(kD)^(1/2)  Energy Frac

aaa0=[0.001,  12.337,   1.66667 ,     2.7207,     5.625*10^-13,     1.08742*10^-20;

0.1009,  12.3366,   1.66667 ,     2.72065,     5.72681*10^-9,     1.12714*10^-12;

0.2008,  12.3361,   1.66667 ,     2.7206,     2.26812*10^-8,     1.768*10^-11;

0.3007,  12.3357,   1.66667 ,     2.72055,     5.08643*10^-8,     8.89155*10^-11;

0.4006,  12.3352,   1.66667 ,     2.7205,     9.02767*10^-8,     2.80093*10^-10;

0.5005,  12.3348,   1.66667 ,     2.72046,     1.40919*10^-7,     6.82482*10^-10;

0.6004,  12.3344,   1.66667 ,     2.72041,     2.02792*10^-7,     1.41336*10^-9;

0.7003,  12.3339,   1.66667 ,     2.72036,     2.75896*10^-7,     2.61603*10^-9;

0.8002,  12.3335,   1.66667 ,     2.72031,     3.60232*10^-7,     4.4598*10^-9;

0.9001,  12.333,   1.66667 ,     2.72026,     4.558*10^-7,     7.14003*10^-9;

1.,  12.3326,   1.66667 ,     2.72021,     5.62601*10^-7,     1.08781*10^-8];


% i    c/L   phL^2/2D   ph/2(kD)^(1/2)  L(k/D)^(1/4) delph/2(kD)^(1/2)   y(0)/d0   energyfrac  (ph/LD)^2-kL^4/D    Imag/Real

aaa1=[1,  0.0001,  12.3326,     1.66667,     2.72021,     5.62601*10^-7,     1.,   1.08781*10^-8,   97.3393,   0.;

2,  0.0101,  12.3326,     1.66666,     2.72022,     5.63024*10^-7,     1.00038,   1.08945*10^-8,   97.3385,   0.;

3,  0.0201,  12.3324,     1.66661,     2.72024,     5.64273*10^-7,     1.0015,   1.09429*10^-8,   97.3327,   0.;

4,  0.0301,  12.332,     1.66647,     2.72031,     5.66345*10^-7,     1.00336,   1.10234*10^-8,   97.3173,   0.;

5,  0.0401,  12.3313,     1.6662,     2.72045,     5.6924*10^-7,     1.00597,   1.11364*10^-8,   97.2874,   0.;

6,  0.0501,  12.33,     1.66575,     2.72068,     5.72965*10^-7,     1.00934,   1.12826*10^-8,   97.2385,   0.;

7,  0.0601,  12.3282,     1.66509,     2.72102,     5.77529*10^-7,     1.01348,   1.1463*10^-8,   97.166,   0.;

8,  0.0701,  12.3256,     1.66417,     2.72149,     5.82945*10^-7,     1.0184,   1.1679*10^-8,   97.0654,   0.;

9,  0.0801,  12.3223,     1.66295,     2.72212,     5.89233*10^-7,     1.02411,   1.19324*10^-8,   96.9326,   0.;

10,  0.0901,  12.3181,     1.6614,     2.72292,     5.96416*10^-7,     1.03064,   1.2225*10^-8,   96.7633,   0.;

11,  0.1001,  12.3129,     1.65946,     2.72393,     6.0452*10^-7,     1.038,   1.25595*10^-8,   96.5535,   0.;

12,  0.1101,  12.3066,     1.65712,     2.72516,     6.1358*10^-7,     1.04622,   1.29388*10^-8,   96.2994,   0.;

13,  0.1201,  12.2992,     1.65433,     2.72663,     6.23633*10^-7,     1.05532,   1.33663*10^-8,   95.9972,   0.;

14,  0.1301,  12.2905,     1.65106,     2.72837,     6.34722*10^-7,     1.06533,   1.38458*10^-8,   95.6434,   0.;

15,  0.1401,  12.2806,     1.64727,     2.7304,     6.46897*10^-7,     1.07629,   1.43821*10^-8,   95.2345,   0.;

16,  0.1501,  12.2694,     1.64294,     2.73275,     6.60212*10^-7,     1.08823,   1.49802*10^-8,   94.7674,   0.;

17,  0.1601,  12.2568,     1.63804,     2.73543,     6.74729*10^-7,     1.1012,   1.56463*10^-8,   94.2389,   0.;

18,  0.1701,  12.2427,     1.63253,     2.73848,     6.90517*10^-7,     1.11523,   1.63871*10^-8,   93.6459,   0.;

19,  0.1801,  12.2273,     1.62639,     2.7419,     7.07652*10^-7,     1.13037,   1.72104*10^-8,   92.9855,   0.;

20,  0.1901,  12.2104,     1.61961,     2.74574,     7.26217*10^-7,     1.14668,   1.81253*10^-8,   92.2549,   0.;

21,  0.2001,  12.192,     1.61215,     2.75001,     7.46306*10^-7,     1.1642,   1.91419*10^-8,   91.4514,   0.;

22,  0.2101,  12.1721,     1.60399,     2.75474,     7.68021*10^-7,     1.18301,   2.02721*10^-8,   90.5723,   0.;

23,  0.2201,  12.1507,     1.59513,     2.75996,     7.91474*10^-7,     1.20317,   2.15291*10^-8,   89.6147,   0.;

24,  0.2301,  12.1278,     1.58554,     2.76568,     8.1679*10^-7,     1.22474,   2.29283*10^-8,   88.5761,   0.;

25,  0.2401,  12.1034,     1.57521,     2.77195,     8.44105*10^-7,     1.24782,   2.44875*10^-8,   87.4537,   0.;

26,  0.2501,  12.0776,     1.56412,     2.77878,     8.73569*10^-7,     1.27247,   2.62268*10^-8,   86.2446,   0.;

27,  0.2601,  12.0503,     1.55228,     2.78621,     9.05347*10^-7,     1.2988,   2.81697*10^-8,   84.9458,   0.;

28,  0.2701,  12.0216,     1.53966,     2.79427,     9.39622*10^-7,     1.3269,   3.0343*10^-8,   83.5543,   0.;

29,  0.2801,  11.9915,     1.52627,     2.80299,     9.76594*10^-7,     1.35688,   3.27778*10^-8,   82.0667,   0.;

30,  0.2901,  11.9599,     1.51209,     2.81239,     1.01648*10^-6,     1.38886,   3.55101*10^-8,   80.4793,   0.;

31,  0.3001,  11.9271,     1.49713,     2.82252,     1.05953*10^-6,     1.42296,   3.85817*10^-8,   78.7883,   0.;

32,  0.3101,  11.8929,     1.48139,     2.8334,     1.10602*10^-6,     1.45932,   4.20411*10^-8,   76.9891,   0.;

33,  0.3201,  11.8574,     1.46487,     2.84508,     1.15622*10^-6,     1.4981,   4.59447*10^-8,   75.0771,   0.;

34,  0.3301,  11.8206,     1.44757,     2.85759,     1.21049*10^-6,     1.53945,   5.03586*10^-8,   73.0468,   0.;

35,  0.3401,  11.7827,     1.4295,     2.87098,     1.26918*10^-6,     1.58357,   5.53598*10^-8,   70.8923,   0.;

36,  0.3501,  11.7435,     1.41066,     2.88528,     1.33268*10^-6,     1.63063,   6.10389*10^-8,   68.6067,   0.;

37,  0.3601,  11.7031,     1.39106,     2.90054,     1.40147*10^-6,     1.68086,   6.75021*10^-8,   66.1824,   0.;

38,  0.3701,  11.6616,     1.3707,     2.91681,     1.47602*10^-6,     1.73448,   7.48746*10^-8,   63.6108,   0.;

39,  0.3801,  11.619,     1.3496,     2.93414,     1.55689*10^-6,     1.79177,   8.33044*10^-8,   60.8822,   0.;

40,  0.3901,  11.5752,     1.32777,     2.95258,     1.6447*10^-6,     1.85299,   9.29666*10^-8,   57.9855,   0.;

41,  0.4001,  11.5303,     1.30522,     2.9722,     1.74014*10^-6,     1.91846,   1.04069*10^-7,   54.9081,   0.;

42,  0.4101,  11.4843,     1.28196,     2.99306,     1.84398*10^-6,     1.98851,   1.16859*10^-7,   51.6357,   0.;

43,  0.4201,  11.4372,     1.258,     3.01523,     1.95706*10^-6,     2.06354,   1.31632*10^-7,   48.1521,   0.;

44,  0.4301,  11.3889,     1.23335,     3.03877,     2.08037*10^-6,     2.14394,   1.48741*10^-7,   44.4385,   0.;

45,  0.4401,  11.3395,     1.20804,     3.06378,     2.21496*10^-6,     2.23019,   1.6861*10^-7,   40.4737,   0.;

46,  0.4501,  11.2889,     1.18206,     3.09034,     2.36205*10^-6,     2.32279,   1.91748*10^-7,   36.233,   0.;

47,  0.4601,  11.237,     1.15544,     3.11855,     2.52301*10^-6,     2.42231,   2.18771*10^-7,   31.6884,   0.;

48,  0.4701,  11.1838,     1.12818,     3.14852,     2.69936*10^-6,     2.52938,   2.50422*10^-7,   26.8073,   0.;

49,  0.4801,  11.1292,     1.1003,     3.18037,     2.89281*10^-6,     2.64471,   2.87602*10^-7,   21.5519,   0.;

50,  0.4901,  11.0731,     1.0718,     3.21422,     3.10533*10^-6,     2.76907,   3.3141*10^-7,   15.8784,   0.;

51,  0.5001,  11.0152,     1.04271,     3.25024,     3.33908*10^-6,     2.90336,   3.83183*10^-7,   9.73558,   0.;

% i    c/L   phL^2/2D   ph/2(kD)^(1/2)  L(k/D)^(1/4) delph/2(kD)^(1/2)   y(0)/d0   energyfrac  (ph/LD)^2-kL^4/D    Imag/Real

1,  0.51,  10.9561,     1.01331,     3.28819,     3.59386*10^-6,     3.04705,   4.4389*10^-7,   3.133,   0.;

2,  0.52,  10.8943,     0.983031,     3.32902,     3.87758*10^-6,     3.20414,   5.16743*10^-7,   -4.13293,   1.09201*10^-12;

3,  0.53,  10.8302,     0.952162,     3.37259,     4.19095*10^-6,     3.37449,   6.03638*10^-7,   -12.0821,   8.90391*10^-13;

4,  0.54,  10.7634,     0.920704,     3.41912,     4.53753*10^-6,     3.55952,   7.07604*10^-7,   -20.8146,   8.64409*10^-14;

5,  0.55,  10.6934,     0.888654,     3.46891,     4.92134*10^-6,     3.76083,   8.32374*10^-7,   -30.4507,   6.45061*10^-14;

6,  0.56,  10.6199,     0.856003,     3.52227,     5.34693*10^-6,     3.98024,   9.82564*10^-7,   -41.1361,   5.02551*10^-14;

7,  0.57,  10.5421,     0.822737,     3.57959,     5.81939*10^-6,     4.21984,   1.16388*10^-6,   -53.0487,   8.39827*10^-15;

8,  0.58,  10.4593,     0.788835,     3.64131,     6.34445*10^-6,     4.48199,   1.38338*10^-6,   -66.4083,   2.16154*10^-14;

9,  0.59,  10.3704,     0.754269,     3.70796,     6.92847*10^-6,     4.76943,   1.64978*10^-6,   -81.4888,   2.35747*10^-14;

10,  0.6,  10.2743,     0.718996,     3.78018,     7.57848*10^-6,     5.0853,   1.97386*10^-6,   -98.6364,   2.02747*10^-14;

11,  0.61,  10.1692,     0.682963,     3.85873,     8.30215*10^-6,     5.43322,   2.36883*10^-6,   -118.294,   2.91495*10^-14;

12,  0.62,  10.0531,     0.646098,     3.94457,     9.1077*10^-6,     5.81741,   2.85082*10^-6,   -141.039,   2.89742*10^-15;

13,  0.63,  9.92303,     0.608303,     4.03889,     0.0000100036,     6.24277,   3.43928*10^-6,   -167.636,   6.40805*10^-14;

14,  0.64,  9.77523,     0.569451,     4.14319,     0.0000109983,     6.71506,   4.15724*10^-6,   -199.119,   1.45044*10^-13;

15,  0.65,  9.60427,     0.529366,     4.25946,     0.0000120992,     7.24106,   5.03116*10^-6,   -236.926,   5.07975*10^-14;

16,  0.66,  9.40245,     0.487809,     4.39031,     0.0000133115,     7.82883,   6.08983*10^-6,   -283.114,   2.7303*10^-13;

17,  0.67,  9.15827,     0.444446,     4.53939,     0.0000146354,     8.48803,   7.3614*10^-6,   -340.735,   1.04357*10^-13;

18,  0.68,  8.85384,     0.398788,     4.71188,     0.0000160622,     9.23045,   8.86665*10^-6,   -414.532,   6.42459*10^-13;

1,  0.684,  8.70916,     0.379729,     4.78907,     0.0000166568,     9.55387,   9.53528*10^-6,   -450.174,   9.77488*10^-13;

2,  0.6841,  8.70534,     0.379246,     4.79107,     0.0000166718,     9.56216,   9.55247*10^-6,   -451.119,   7.96172*10^-13;

3,  0.6842,  8.70151,     0.378762,     4.79307,     0.0000166868,     9.57047,   9.56968*10^-6,   -452.067,   1.06454*10^-12;

4,  0.6843,  8.69767,     0.378278,     4.79508,     0.0000167018,     9.57878,   9.58691*10^-6,   -453.018,   1.05322*10^-12;

5,  0.6844,  8.69382,     0.377794,     4.79709,     0.0000167168,     9.58711,   9.60417*10^-6,   -453.972,   2.41464*10^-13;

6,  0.6845,  8.68996,     0.37731,     4.7991,     0.0000167319,     9.59545,   9.62145*10^-6,   -454.928,   4.334*10^-13;

7,  0.6846,  8.68608,     0.376825,     4.80112,     0.0000167469,     9.60379,   9.63875*10^-6,   -455.888,   1.3641*10^-12;

8,  0.6847,  8.6822,     0.37634,     4.80314,     0.0000167619,     9.61215,   9.65607*10^-6,   -456.85,   3.51792*10^-13;

9,  0.6848,  8.67831,     0.375854,     4.80516,     0.000016777,     9.62052,   9.67342*10^-6,   -457.815,   8.85455*10^-13;

10,  0.6849,  8.6744,     0.375368,     4.80719,     0.000016792,     9.6289,   9.69078*10^-6,   -458.783,   1.17034*10^-12;

11,  0.685,  8.67048,     0.374882,     4.80922,     0.0000168071,     9.63729,   9.70817*10^-6,   -459.754,   6.14042*10^-13;

12,  0.6851,  8.66656,     0.374395,     4.81125,     0.0000168222,     9.64569,   9.72559*10^-6,   -460.728,   2.95106*10^-13;

13,  0.6852,  8.66262,     0.373908,     4.81329,     0.0000168372,     9.6541,   9.74302*10^-6,   -461.704,   8.45993*10^-13;

14,  0.6853,  8.65867,     0.373421,     4.81533,     0.0000168523,     9.66252,   9.76048*10^-6,   -462.684,   1.57841*10^-12;

15,  0.6854,  8.65471,     0.372933,     4.81738,     0.0000168674,     9.67095,   9.77796*10^-6,   -463.667,   1.58525*10^-13;

16,  0.6855,  8.65074,     0.372445,     4.81943,     0.0000168825,     9.6794,   9.79546*10^-6,   -464.652,   1.42434*10^-13;

17,  0.6856,  8.64676,     0.371957,     4.82148,     0.0000168976,     9.68785,   9.81298*10^-6,   -465.641,   2.95076*10^-13;

18,  0.6857,  8.64276,     0.371468,     4.82354,     0.0000169127,     9.69631,   9.83053*10^-6,   -466.632,   8.54829*10^-13;

19,  0.6858,  8.63876,     0.370979,     4.8256,     0.0000169278,     9.70479,   9.84809*10^-6,   -467.627,   1.78397*10^-12];
