%    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.72069,     5.62501*10^-11,     1.08742*10^-18;

0.1009,  12.3325,   1.66667 ,     2.72021,     5.72774*10^-7,     1.12751*10^-10;

0.2008,  12.3281,   1.66667 ,     2.71972,     2.26885*10^-6,     1.76915*10^-9;

0.3007,  12.3237,   1.66667 ,     2.71923,     5.0889*10^-6,     8.9002*10^-9;

0.4006,  12.3193,   1.66667 ,     2.71875,     9.03353*10^-6,     2.80457*10^-8;

0.5005,  12.3149,   1.66667 ,     2.71826,     0.0000141033,     6.83591*10^-8;

0.6004,  12.3105,   1.66667 ,     2.71777,     0.000020299,     1.41612*10^-7;

0.7003,  12.306,   1.66667 ,     2.71728,     0.000027621,     2.62199*10^-7;

0.8002,  12.3016,   1.66667 ,     2.7168,     0.0000360701,     4.47144*10^-7;

0.9001,  12.2972,   1.66667 ,     2.71631,     0.0000456469,     7.16101*10^-7;

1.,  12.2928,   1.66667 ,     2.71582,     0.000056352,     1.09136*10^-6];


%  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.2928,     1.66667,     2.71582,     0.000056352,     1.,   1.09137*10^-6,   96.7121,   0.;

2,  0.0101,  12.2927,     1.66666,     2.71582,     0.0000563943,     1.00038,   1.093*10^-6,   96.711,   0.;

3,  0.0201,  12.2926,     1.66661,     2.71584,     0.0000565191,     1.00149,   1.09785*10^-6,   96.7046,   0.;

4,  0.0301,  12.2921,     1.66647,     2.7159,     0.0000567261,     1.00335,   1.1059*10^-6,   96.6881,   0.;

5,  0.0401,  12.2912,     1.6662,     2.71603,     0.0000570154,     1.00595,   1.11721*10^-6,   96.6569,   0.;

6,  0.0501,  12.2899,     1.66575,     2.71624,     0.0000573876,     1.00931,   1.13185*10^-6,   96.6064,   0.;

7,  0.0601,  12.2879,     1.6651,     2.71656,     0.0000578436,     1.01344,   1.14991*10^-6,   96.5321,   0.;

8,  0.0701,  12.2852,     1.66418,     2.717,     0.0000583848,     1.01834,   1.17153*10^-6,   96.4295,   0.;

9,  0.0801,  12.2816,     1.66297,     2.7176,     0.0000590131,     1.02403,   1.19687*10^-6,   96.2945,   0.;

10,  0.0901,  12.2771,     1.66142,     2.71837,     0.0000597307,     1.03053,   1.22616*10^-6,   96.1229,   0.;

11,  0.1001,  12.2717,     1.6595,     2.71934,     0.0000605404,     1.03787,   1.25963*10^-6,   95.9108,   0.;

12,  0.1101,  12.2651,     1.65717,     2.72052,     0.0000614454,     1.04606,   1.29757*10^-6,   95.6542,   0.;

13,  0.1201,  12.2573,     1.65439,     2.72194,     0.0000624495,     1.05512,   1.34032*10^-6,   95.3497,   0.;

14,  0.1301,  12.2484,     1.65114,     2.72362,     0.0000635571,     1.0651,   1.38829*10^-6,   94.9936,   0.;

15,  0.1401,  12.2381,     1.64737,     2.72559,     0.0000647729,     1.07601,   1.44191*10^-6,   94.5825,   0.;

16,  0.1501,  12.2264,     1.64306,     2.72786,     0.0000661025,     1.08791,   1.50171*10^-6,   94.1134,   0.;

17,  0.1601,  12.2134,     1.63818,     2.73046,     0.0000675519,     1.10082,   1.56829*10^-6,   93.583,   0.;

18,  0.1701,  12.1989,     1.63271,     2.73342,     0.000069128,     1.11478,   1.64233*10^-6,   92.9883,   0.;

19,  0.1801,  12.1829,     1.62661,     2.73675,     0.0000708382,     1.12985,   1.72459*10^-6,   92.3267,   0.;

20,  0.1901,  12.1655,     1.61986,     2.74047,     0.0000726909,     1.14608,   1.81598*10^-6,   91.5952,   0.;

21,  0.2001,  12.1465,     1.61244,     2.74463,     0.0000746952,     1.16352,   1.91751*10^-6,   90.7911,   0.;

22,  0.2101,  12.126,     1.60434,     2.74923,     0.0000768613,     1.18223,   2.03033*10^-6,   89.912,   0.;

23,  0.2201,  12.1039,     1.59553,     2.7543,     0.0000792003,     1.20227,   2.15578*10^-6,   88.9551,   0.;

24,  0.2301,  12.0803,     1.586,     2.75987,     0.0000817244,     1.22372,   2.29538*10^-6,   87.9178,   0.;

25,  0.2401,  12.0552,     1.57573,     2.76597,     0.0000844471,     1.24665,   2.45087*10^-6,   86.7975,   0.;

26,  0.2501,  12.0286,     1.56472,     2.77262,     0.0000873831,     1.27115,   2.62426*10^-6,   85.5914,   0.;

27,  0.2601,  12.0005,     1.55295,     2.77985,     0.0000905488,     1.2973,   2.81784*10^-6,   84.2969,   0.;

28,  0.2701,  11.9709,     1.54042,     2.78768,     0.0000939619,     1.3252,   3.03428*10^-6,   82.9108,   0.;

29,  0.2801,  11.9398,     1.52712,     2.79616,     0.0000976422,     1.35496,   3.27663*10^-6,   81.43,   0.;

30,  0.2901,  11.9073,     1.51305,     2.8053,     0.000101611,     1.38668,   3.54843*10^-6,   79.8512,   0.;

31,  0.3001,  11.8734,     1.49821,     2.81514,     0.000105893,     1.4205,   3.85379*10^-6,   78.1707,   0.;

32,  0.3101,  11.838,     1.4826,     2.82572,     0.000110514,     1.45654,   4.19746*10^-6,   76.3843,   0.;

33,  0.3201,  11.8014,     1.46621,     2.83706,     0.000115503,     1.49495,   4.58497*10^-6,   74.4876,   0.;

34,  0.3301,  11.7634,     1.44906,     2.8492,     0.000120891,     1.53589,   5.02277*10^-6,   72.4757,   0.;

35,  0.3401,  11.724,     1.43114,     2.86218,     0.000126716,     1.57952,   5.51839*10^-6,   70.343,   0.;

36,  0.3501,  11.6834,     1.41247,     2.87604,     0.000133014,     1.62605,   6.08062*10^-6,   68.0833,   0.;

37,  0.3601,  11.6416,     1.39306,     2.89082,     0.000139831,     1.67567,   6.71981*10^-6,   65.6897,   0.;

38,  0.3701,  11.5985,     1.37291,     2.90656,     0.000147213,     1.7286,   7.44808*10^-6,   63.1542,   0.;

39,  0.3801,  11.5541,     1.35203,     2.92332,     0.000155215,     1.78509,   8.27974*10^-6,   60.4682,   0.;

40,  0.3901,  11.5086,     1.33044,     2.94113,     0.000163894,     1.84541,   9.23167*10^-6,   57.6215,   0.;

41,  0.4001,  11.4618,     1.30815,     2.96005,     0.000173319,     1.90984,   0.0000103239,   54.6029,   0.;

42,  0.4101,  11.4138,     1.28517,     2.98014,     0.00018356,     1.97871,   0.00001158,   51.3997,   0.;

43,  0.4201,  11.3646,     1.26151,     3.00145,     0.000194701,     2.05236,   0.0000130282,   47.9974,   0.;

44,  0.4301,  11.3141,     1.2372,     3.02406,     0.000206831,     2.13119,   0.0000147022,   44.3795,   0.;

45,  0.4401,  11.2624,     1.21224,     3.04803,     0.000220053,     2.21562,   0.000016642,   40.5273,   0.;

46,  0.4501,  11.2093,     1.18666,     3.07344,     0.00023448,     2.30612,   0.0000188957,   36.4194,   0.;

47,  0.4601,  11.1548,     1.16046,     3.10038,     0.000250239,     2.4032,   0.0000215209,   32.0316,   0.;

48,  0.4701,  11.0988,     1.13367,     3.12893,     0.000267471,     2.50745,   0.000024587,   27.3359,   0.;

49,  0.4801,  11.0413,     1.10629,     3.15919,     0.000286336,     2.61948,   0.0000281775,   22.3004,   0.;

50,  0.4901,  10.9821,     1.07835,     3.19127,     0.000307011,     2.74,   0.0000323936,   16.8883,   0.;

51,  0.5001,  10.9211,     1.04985,     3.2253,     0.000329695,     2.8698,   0.0000373574,   11.057,   0.;

52,  0.5101,  10.8581,     1.02081,     3.2614,     0.000354612,     3.00973,   0.0000432175,   4.75719,   0.;

53,  0.5201,  10.7928,     0.991238,     3.29973,     0.000382012,     3.16075,   0.000050154,   -2.06834,   1.42832*10^-13;

54,  0.5301,  10.7251,     0.961151,     3.34046,     0.000412174,     3.32394,   0.0000583865,   -9.48665,   2.09613*10^-13;

55,  0.5401,  10.6547,     0.930554,     3.38376,     0.00044541,     3.5005,   0.0000681823,   -17.5764,   1.01828*10^-13;

56,  0.5501,  10.5811,     0.899452,     3.42986,     0.000482069,     3.69175,   0.0000798676,   -26.4307,   1.33292*10^-14;

57,  0.5601,  10.504,     0.86785,     3.479,     0.00052254,     3.89918,   0.0000938405,   -36.1599,   1.09909*10^-13;

58,  0.5701,  10.4227,     0.835748,     3.53145,     0.000567251,     4.12448,   0.000110587,   -46.8961,   2.67225*10^-14;

59,  0.5801,  10.3367,     0.803141,     3.58753,     0.000616679,     4.36951,   0.000130698,   -58.7988,   2.49653*10^-14;

60,  0.5901,  10.2452,     0.77002,     3.64762,     0.000671342,     4.63638,   0.000154896,   -72.0616,   1.97256*10^-14;

61,  0.6001,  10.1471,     0.73637,     3.71213,     0.000731805,     4.92743,   0.000184053,   -86.9223,   2.77147*10^-14;

62,  0.6101,  10.0413,     0.702167,     3.7816,     0.000798673,     5.24534,   0.000219225,   -103.676,   6.5097*10^-14;

63,  0.6201,  9.92624,     0.667377,     3.85662,     0.000872581,     5.59308,   0.000261675,   -122.691,   7.20961*10^-14;

64,  0.6301,  9.79996,     0.631954,     3.93794,     0.000954175,     5.97402,   0.000312902,   -144.44,   4.14579*10^-14;

65,  0.6401,  9.65997,     0.595834,     4.02648,     0.00104409,     6.39195,   0.000374651,   -169.531,   8.38684*10^-15;

66,  0.6501,  9.50297,     0.558932,     4.12335,     0.00114289,     6.85116,   0.00044891,   -198.762,   7.48798*10^-14;

67,  0.6601,  9.32455,     0.52113,     4.23,     0.001251,     7.3565,   0.000537855,   -233.21,   1.91401*10^-14;

68,  0.6701,  9.11866,     0.482271,     4.3483,     0.00136859,     7.91345,   0.00064372,   -274.352,   1.27142*10^-13;

69,  0.6801,  8.87669,     0.442135,     4.48072,     0.00149538,     8.52825,   0.000768517,   -324.286,   9.19091*10^-14;

70,  0.6901,  8.58608,     0.400411,     4.63068,     0.00163035,     9.20803,   0.00091351,   -386.09,   2.15747*10^-13;

1,  0.695,  8.42045,     0.379257,     4.71195,     0.00169886,     9.56729,   0.000991902,   -422.046,   4.36027*10^-13;

2,  0.6951,  8.41688,     0.37882,     4.71367,     0.00170027,     9.57481,   0.00099355,   -422.826,   5.72871*10^-13;

3,  0.6952,  8.4133,     0.378382,     4.71539,     0.00170169,     9.58234,   0.0009952,   -423.608,   7.75602*10^-14;

4,  0.6953,  8.40971,     0.377945,     4.71712,     0.0017031,     9.58988,   0.000996852,   -424.392,   2.63082*10^-13;

5,  0.6954,  8.40612,     0.377507,     4.71884,     0.00170451,     9.59743,   0.000998506,   -425.178,   7.59438*10^-13;

6,  0.6955,  8.40252,     0.377068,     4.72057,     0.00170592,     9.60498,   0.00100016,   -425.966,   8.0704*10^-13;

7,  0.6956,  8.39891,     0.37663,     4.7223,     0.00170733,     9.61254,   0.00100182,   -426.755,   9.69901*10^-13;

8,  0.6957,  8.39529,     0.376191,     4.72404,     0.00170875,     9.62011,   0.00100348,   -427.547,   5.18795*10^-13;

9,  0.6958,  8.39167,     0.375753,     4.72578,     0.00171016,     9.62769,   0.00100514,   -428.341,   8.37502*10^-13;

10,  0.6959,  8.38803,     0.375313,     4.72752,     0.00171158,     9.63527,   0.0010068,   -429.137,   4.55876*10^-13;

11,  0.696,  8.38439,     0.374874,     4.72926,     0.00171299,     9.64287,   0.00100847,   -429.935,   8.65354*10^-13;

12,  0.6961,  8.38074,     0.374434,     4.731,     0.00171441,     9.65047,   0.00101013,   -430.735,   1.2497*10^-12;

13,  0.6962,  8.37708,     0.373995,     4.73275,     0.00171582,     9.65807,   0.0010118,   -431.536,   2.42978*10^-12;

14,  0.6963,  8.37341,     0.373555,     4.7345,     0.00171724,     9.66569,   0.00101347,   -432.34,   6.61344*10^-13;

15,  0.6964,  8.36974,     0.373114,     4.73625,     0.00171865,     9.67331,   0.00101515,   -433.146,   1.00961*10^-12;

16,  0.6965,  8.36605,     0.372674,     4.73801,     0.00172007,     9.68095,   0.00101682,   -433.955,   6.32612*10^-13;

17,  0.6966,  8.36236,     0.372233,     4.73977,     0.00172149,     9.68859,   0.00101849,   -434.765,   8.89444*10^-14;

18,  0.6967,  8.35866,     0.371792,     4.74153,     0.0017229,     9.69623,   0.00102017,   -435.577,   5.55895*10^-13;

19,  0.6968,  8.35495,     0.371351,     4.74329,     0.00172432,     9.70389,   0.00102185,   -436.391,   3.54613*10^-13];