diff --git a/data.txt b/data.txt
index 27c3ead..3b95842 100644
--- a/data.txt
+++ b/data.txt
@@ -1,1237 +1,1237 @@
 wl	n
-0.40	1.41491
-0.41	1.41403
-0.42	1.41300
-0.43	1.41179
-0.44	1.41113
-0.45	1.41033
-0.46	1.40943
-0.47	1.40866
-0.48	1.40802
-0.49	1.40733
-0.50	1.40675
-0.51	1.40614
-0.52	1.40557
-0.53	1.40513
-0.54	1.40430
-0.55	1.40420
-0.56	1.40375
-0.57	1.40338
-0.58	1.40292
-0.59	1.40259
-0.60	1.40202
-0.61	1.40158
-0.62	1.40139
-0.63	1.40121
-0.64	1.40059
-0.65	1.40047
-0.66	1.40022
-0.67	1.40006
-0.68	1.39984
-0.69	1.39958
-0.70	1.39907
-0.71	1.39905
-0.72	1.39871
-0.73	1.39817
-0.74	1.39814
-0.75	1.39765
-0.76	1.39750
-0.77	1.39727
-0.78	1.39660
-0.79	1.39645
-0.80	1.39601
-0.81	1.39617
-0.82	1.39684
-0.83	1.39701
-0.84	1.39629
-0.85	1.39687
-0.86	1.39671
-0.87	1.39652
-0.88	1.39624
-0.89	1.39629
-0.90	1.39616
-0.91	1.39610
-0.92	1.39583
-0.93	1.39584
-0.94	1.39554
-0.95	1.39546
-0.96	1.39552
-0.97	1.39537
-0.98	1.39513
-0.99	1.39506
-1.00	1.39525
-1.01	1.39523
-1.02	1.39544
-1.03	1.39121
-1.04	1.39127
-1.05	1.39181
-1.06	1.39102
-1.07	1.39193
-1.08	1.39217
-1.09	1.39220
-1.10	1.39140
-1.11	1.39202
-1.12	1.39164
-1.13	1.39084
-1.14	1.39105
-1.15	1.39087
-1.16	1.39135
-1.17	1.39134
-1.18	1.39199
-1.19	1.39145
-1.20	1.39106
-1.21	1.38948
-1.22	1.39053
-1.23	1.38988
-1.24	1.38929
-1.25	1.38975
-1.26	1.38999
-1.27	1.39049
-1.28	1.39009
-1.29	1.38910
-1.30	1.38878
-1.31	1.38850
-1.32	1.38844
-1.33	1.38895
-1.34	1.38870
-1.35	1.38897
-1.36	1.38914
-1.37	1.39000
-1.38	1.38972
-1.39	1.38923
-1.40	1.38938
-1.41	1.38836
-1.42	1.38932
-1.43	1.38933
-1.44	1.38929
-1.45	1.38840
-1.46	1.38840
-1.47	1.38895
-1.48	1.38968
-1.49	1.38914
-1.50	1.38886
-1.51	1.38844
-1.52	1.38850
-1.53	1.38804
-1.54	1.38863
-1.55	1.38922
-1.56	1.38841
-1.57	1.38815
-1.58	1.38787
-1.59	1.38755
-1.60	1.38748
-1.61	1.38731
-1.62	1.38811
-1.63	1.38725
-1.64	1.38609
-1.65	1.38759
-1.66	1.38751
-1.67	1.38797
-1.68	1.38773
-1.69	1.38719
-1.70	1.38694
-1.71	1.38755
-1.72	1.38689
-1.73	1.38756
-1.74	1.38736
-1.75	1.38736
-1.76	1.38611
-1.77	1.38703
-1.78	1.38651
-1.79	1.38645
-1.80	1.38644
-1.81	1.38680
-1.82	1.38704
-1.83	1.38639
-1.84	1.38730
-1.85	1.38712
-1.86	1.38634
-1.87	1.38571
-1.88	1.38651
-1.89	1.38631
-1.90	1.38600
-1.91	1.38610
-1.92	1.38585
-1.93	1.38623
-1.94	1.38625
-1.95	1.38571
-1.96	1.38546
-1.97	1.38623
-1.98	1.38586
-1.99	1.38653
-2.00	1.38648
-2.0097	1.39924
-2.0190	1.39793
-2.0285	1.39715
-2.0381	1.39551
-2.0478	1.39481
-2.0575	1.39529
-2.0673	1.39562
-2.0773	1.39677
-2.0873	1.39707
-2.0975	1.39523
-2.1077	1.39490
-2.1180	1.39475
-2.1285	1.39453
-2.1390	1.39614
-2.1496	1.39796
-2.1568	1.39915
-2.1676	1.39923
-2.1785	1.39721
-2.1896	1.39679
-2.1970	1.39814
-2.2082	1.40006
-2.2196	1.40093
-2.2272	1.40083
-2.2387	1.39913
-2.2465	1.39783
-2.2582	1.39753
-2.2661	1.39782
-2.2781	1.39884
-2.2861	1.39983
-2.2983	1.39921
-2.3065	1.39840
-2.3188	1.39812
-2.3272	1.39817
-2.3398	1.39834
-2.3482	1.39873
-2.3568	1.39875
-2.3697	1.39769
-2.3784	1.39710
-2.3872	1.39638
-2.3960	1.39638
-2.4093	1.39737
-2.4183	1.39788
-2.4274	1.39844
-2.4365	1.39870
-2.4457	1.39872
-2.4596	1.39708
-2.4690	1.39680
-2.4784	1.39698
-2.4880	1.39703
-2.4975	1.39646
-2.5072	1.39605
-2.5169	1.39625
-2.5268	1.39632
-2.5366	1.39528
-2.5466	1.39446
-2.5567	1.39468
-2.5668	1.39547
-2.5770	1.39589
-2.5873	1.39548
-2.5976	1.39486
-2.6081	1.39458
-2.6186	1.39407
-2.6293	1.39401
-2.6400	1.39501
-2.6454	1.39538
-2.6562	1.39534
-2.6671	1.39507
-2.6782	1.39542
-2.6893	1.39518
-2.6949	1.39492
-2.7061	1.39514
-2.7175	1.39539
-2.7289	1.39509
-2.7347	1.39508
-2.7462	1.39535
-2.7579	1.39543
-2.7697	1.39596
-2.7756	1.39632
-2.7876	1.39618
-2.7996	1.39533
-2.8057	1.39520
-2.8179	1.39570
-2.8240	1.39604
-2.8364	1.39633
-2.8489	1.39575
-2.8551	1.39494
-2.8678	1.39319
-2.8741	1.39313
-2.8869	1.39458
-2.8998	1.39513
-2.9063	1.39478
-2.9194	1.39497
-2.9260	1.39484
-2.9393	1.39403
-2.9460	1.39363
-2.9594	1.39380
-2.9662	1.39377
-2.9798	1.39379
-2.9867	1.39366
-2.9936	1.39342
-3.0075	1.39390
-3.0145	1.39431
-3.0286	1.39431
-3.0357	1.39427
-3.0499	1.39371
-3.0571	1.39288
-3.0644	1.39204
-3.0789	1.39114
-3.0863	1.39119
-3.0936	1.39144
-3.1085	1.39225
-3.1159	1.39210
-3.1234	1.39168
-3.1386	1.39114
-3.1462	1.39112
-3.1538	1.39088
-3.1693	1.39129
-3.1770	1.39169
-3.1848	1.39173
-3.1927	1.39124
-3.2085	1.39036
-3.2164	1.38993
-3.2244	1.38921
-3.2325	1.38876
-3.2487	1.38817
-3.2569	1.38737
-3.2651	1.38650
-3.2733	1.38627
-3.2899	1.38658
-3.2983	1.38591
-3.3067	1.38448
-3.3152	1.38272
-3.3237	1.38100
-3.3322	1.37863
-3.3494	1.37097
-3.3581	1.36967
-3.3668	1.37687
-3.3756	1.39316
-3.3844	1.40866
-3.3933	1.41427
-3.4022	1.41184
-3.4111	1.40699
-3.4292	1.40079
-3.4383	1.40068
-3.4474	1.40129
-3.4566	1.40155
-3.4658	1.40126
-3.4751	1.40054
-3.4845	1.39949
-3.4939	1.39857
-3.5033	1.39804
-3.5128	1.39751
-3.5224	1.39708
-3.5320	1.39706
-3.5416	1.39677
-3.5513	1.39639
-3.5611	1.39644
-3.5709	1.39645
-3.5807	1.39613
-3.5907	1.39573
-3.6006	1.39526
-3.6107	1.39488
-3.6207	1.39425
-3.6309	1.39387
-3.6411	1.39371
-3.6513	1.39336
-3.6617	1.39299
-3.6720	1.39275
-3.6825	1.39279
-3.6930	1.39285
-3.7035	1.39292
-3.7141	1.39338
-3.7248	1.39353
-3.7355	1.39284
-3.7463	1.39221
-3.7572	1.39204
-3.7681	1.39178
-3.7791	1.39137
-3.7901	1.39117
-3.8013	1.39096
-3.8124	1.39052
-3.8237	1.39012
-3.8350	1.39015
-3.8464	1.39046
-3.8578	1.39035
-3.8693	1.38984
-3.8809	1.38949
-3.8926	1.38922
-3.9043	1.38909
-3.9161	1.38911
-3.9280	1.38899
-3.9399	1.38885
-3.9519	1.38895
-3.9640	1.38901
-3.9762	1.38868
-3.9884	1.38831
-4.0007	1.38835
-4.0131	1.38836
-4.0255	1.38825
-4.0381	1.38827
-4.0507	1.38831
-4.0634	1.38821
-4.0762	1.38796
-4.0890	1.38747
-4.1020	1.38699
-4.1150	1.38686
-4.1281	1.38674
-4.1413	1.38644
-4.1546	1.38645
-4.1679	1.38652
-4.1814	1.38611
-4.1949	1.38578
-4.2085	1.38612
-4.2222	1.38707
-4.2360	1.38743
-4.2499	1.38623
-4.2639	1.38489
-4.2780	1.38445
-4.2921	1.38487
-4.3064	1.38482
-4.3208	1.38393
-4.3352	1.38330
-4.3498	1.38299
-4.3644	1.38257
-4.3791	1.38213
-4.3940	1.38189
-4.4089	1.38216
-4.4240	1.38250
-4.4391	1.38217
-4.4544	1.38162
-4.4697	1.38145
-4.4852	1.38125
-4.5008	1.38079
-4.5165	1.38026
-4.5323	1.37972
-4.5482	1.37925
-4.5642	1.37880
-4.5803	1.37837
-4.5965	1.37807
-4.6129	1.37826
-4.6294	1.37878
-4.6460	1.37978
-4.6627	1.38045
-4.6795	1.38033
-4.6965	1.37971
-4.7136	1.37893
-4.7308	1.37851
-4.7481	1.37861
-4.7655	1.37872
-4.7831	1.37836
-4.8008	1.37783
-4.8187	1.37745
-4.8367	1.37699
-4.8548	1.37667
-4.8730	1.37645
-4.8914	1.37643
-4.9100	1.37627
-4.9286	1.37551
-4.9474	1.37478
-4.9664	1.37452
-4.9855	1.37413
-5.0047	1.37342
-5.0241	1.37301
-5.0437	1.37286
-5.0634	1.37230
-5.0832	1.37140
-5.1033	1.37101
-5.1234	1.37112
-5.1438	1.37151
-5.1643	1.37161
-5.1849	1.37116
-5.2057	1.37061
-5.2267	1.37010
-5.2479	1.36974
-5.2692	1.36945
-5.2907	1.36930
-5.3124	1.36857
-5.3343	1.36819
-5.3563	1.36785
-5.3785	1.36698
-5.4009	1.36623
-5.4235	1.36576
-5.4463	1.36563
-5.4693	1.36493
-5.4925	1.36400
-5.5159	1.36376
-5.5394	1.36366
-5.5632	1.36332
-5.5872	1.36268
-5.6114	1.36188
-5.6358	1.36124
-5.6604	1.36063
-5.6852	1.36017
-5.7103	1.35952
-5.7355	1.35851
-5.7610	1.35788
-5.7867	1.35756
-5.8127	1.35700
-5.8389	1.35636
-5.8653	1.35602
-5.8919	1.35546
-5.9188	1.35444
-5.9460	1.35327
-5.9734	1.35208
-6.0011	1.35141
-6.0290	1.35101
-6.0571	1.34985
-6.0856	1.34856
-6.1143	1.34770
-6.1433	1.34668
-6.1725	1.34572
-6.2020	1.34448
-6.2319	1.34381
-6.2620	1.34327
-6.2924	1.34257
-6.3231	1.34188
-6.3541	1.34096
-6.3854	1.33948
-6.4170	1.33844
-6.4489	1.33698
-6.4811	1.33518
-6.5137	1.33420
-6.5466	1.33287
-6.5798	1.33102
-6.6134	1.32961
-6.6473	1.32711
-6.6816	1.32364
-6.7162	1.32168
-6.7512	1.32053
-6.7865	1.31839
-6.8223	1.31494
-6.8583	1.31223
-6.8948	1.31191
-6.9317	1.31155
-6.9690	1.30835
-7.0066	1.30306
-7.0447	1.30038
-7.0832	1.30332
-7.1221	1.30790
-7.1615	1.30939
-7.2013	1.30643
-7.2415	1.30087
-7.2822	1.29479
-7.3233	1.28942
-7.3649	1.28363
-7.4070	1.27697
-7.4496	1.27066
-7.4926	1.26412
-7.5362	1.25650
-7.5803	1.24758
-7.6249	1.23651
-7.6700	1.22417
-7.7156	1.20747
-7.7618	1.18420
-7.8086	1.13872
-7.8559	1.06920
-7.9038	1.13086
-7.9523	1.31619
-8.0014	1.42903
-8.0511	1.38024
-8.1014	1.31925
-8.1524	1.28177
-8.2040	1.25491
-8.2562	1.23390
-8.3091	1.21820
-8.3628	1.20668
-8.4171	1.19571
-8.4721	1.18132
-8.5278	1.16258
-8.5843	1.14039
-8.6415	1.11398
-8.6995	1.07840
-8.7583	1.03705
-8.8179	0.99077
-8.8783	0.94243
-8.9395	0.90850
-9.0016	0.91694
-9.0645	0.99296
-9.1284	1.12582
-9.1931	1.27070
-9.2588	1.38675
-9.3254	1.46460
-9.3930	1.51165
-9.4615	1.53749
-9.5311	1.55207
-9.6017	1.56879
-9.6733	1.60177
-9.7461	1.66994
-9.8199	1.77468
-9.8949	1.85128
-9.9710	1.83876
-10.048	1.77344
-10.127	1.70999
-10.207	1.66025
-10.288	1.61895
-10.370	1.58348
-10.453	1.55340
-10.538	1.52658
-10.625	1.50134
-10.713	1.47663
-10.802	1.45099
-10.893	1.42981
-10.985	1.41565
-11.079	1.39727
-11.174	1.36630
-11.272	1.32958
-11.370	1.30103
-11.471	1.30308
-11.573	1.32585
-11.678	1.33199
-11.784	1.33215
-11.892	1.33290
-12.002	1.28572
-12.114	1.23134
-12.229	1.25325
-12.345	1.37130
-12.464	1.60462
-12.585	1.86355
-12.708	1.96263
-12.834	1.89414
-12.962	1.79653
-13.093	1.73460
-13.227	1.71230
-13.363	1.69554
-13.502	1.66481
-13.645	1.63533
-13.790	1.61279
-13.938	1.59113
-14.089	1.57651
-14.244	1.57656
-14.403	1.58157
-14.564	1.58549
-14.730	1.58457
-14.899	1.57902
-15.072	1.57898
-15.250	1.57744
-15.431	1.56572
-15.617	1.55054
-15.808	1.53957
-16.003	1.53249
-16.203	1.52491
-16.408	1.51758
-16.618	1.51612
-16.834	1.51225
-17.056	1.50642
-17.283	1.50336
-17.517	1.50038
-17.757	1.49764
-18.003	1.50001
-18.257	1.50080
-18.518	1.49350
-18.786	1.49025
-19.062	1.49254
-19.347	1.49056
-19.640	1.48533
-19.942	1.48051
+0.40	1.42109
+0.41	1.41984
+0.42	1.41859
+0.43	1.41780
+0.44	1.41642
+0.45	1.41581
+0.46	1.41499
+0.47	1.41395
+0.48	1.41342
+0.49	1.41243
+0.50	1.41189
+0.51	1.41128
+0.52	1.41051
+0.53	1.41000
+0.54	1.40954
+0.55	1.40929
+0.56	1.40851
+0.57	1.40819
+0.58	1.40776
+0.59	1.40739
+0.60	1.40715
+0.61	1.40653
+0.62	1.40620
+0.63	1.40553
+0.64	1.40527
+0.65	1.40507
+0.66	1.40484
+0.67	1.40451
+0.68	1.40453
+0.69	1.40389
+0.70	1.40397
+0.71	1.40335
+0.72	1.40324
+0.73	1.40313
+0.74	1.40307
+0.75	1.40234
+0.76	1.40235
+0.77	1.40206
+0.78	1.40201
+0.79	1.40158
+0.80	1.40124
+0.81	1.40126
+0.82	1.40150
+0.83	1.40150
+0.84	1.40137
+0.85	1.40153
+0.86	1.40137
+0.87	1.40120
+0.88	1.40091
+0.89	1.40086
+0.90	1.40077
+0.91	1.40061
+0.92	1.40047
+0.93	1.40040
+0.94	1.40025
+0.95	1.40015
+0.96	1.40015
+0.97	1.39997
+0.98	1.39977
+0.99	1.39970
+1.00	1.39976
+1.01	1.39982
+1.02	1.39993
+1.03	1.39734
+1.04	1.39708
+1.05	1.39739
+1.06	1.39685
+1.07	1.39677
+1.08	1.39696
+1.09	1.39675
+1.10	1.39673
+1.11	1.39660
+1.12	1.39654
+1.13	1.39622
+1.14	1.39568
+1.15	1.39599
+1.16	1.39631
+1.17	1.39619
+1.18	1.39656
+1.19	1.39597
+1.20	1.39567
+1.21	1.39517
+1.22	1.39569
+1.23	1.39558
+1.24	1.39494
+1.25	1.39481
+1.26	1.39515
+1.27	1.39555
+1.28	1.39490
+1.29	1.39445
+1.30	1.39483
+1.31	1.39424
+1.32	1.39390
+1.33	1.39453
+1.34	1.39340
+1.35	1.39398
+1.36	1.39455
+1.37	1.39453
+1.38	1.39425
+1.39	1.39359
+1.40	1.39387
+1.41	1.39361
+1.42	1.39392
+1.43	1.39373
+1.44	1.39366
+1.45	1.39347
+1.46	1.39318
+1.47	1.39373
+1.48	1.39400
+1.49	1.39380
+1.50	1.39423
+1.51	1.39279
+1.52	1.39283
+1.53	1.39321
+1.54	1.39284
+1.55	1.39328
+1.56	1.39309
+1.57	1.39291
+1.58	1.39231
+1.59	1.39276
+1.60	1.39251
+1.61	1.39299
+1.62	1.39257
+1.63	1.39185
+1.64	1.39271
+1.65	1.39170
+1.66	1.39185
+1.67	1.39190
+1.68	1.39262
+1.69	1.39197
+1.70	1.39116
+1.71	1.39032
+1.72	1.39027
+1.73	1.39079
+1.74	1.39191
+1.75	1.39224
+1.76	1.39094
+1.77	1.39111
+1.78	1.39197
+1.79	1.39189
+1.80	1.39084
+1.81	1.39131
+1.82	1.39026
+1.83	1.39172
+1.84	1.39211
+1.85	1.39078
+1.86	1.39137
+1.87	1.39159
+1.88	1.39150
+1.89	1.39093
+1.90	1.39098
+1.91	1.39019
+1.92	1.39004
+1.93	1.39071
+1.94	1.39044
+1.95	1.39051
+1.96	1.39081
+1.97	1.39043
+1.98	1.39098
+1.99	1.38973
+2.00	1.39092
+2.0097	1.40028
+2.0190	1.39960
+2.0285	1.39855
+2.0381	1.39720
+2.0478	1.39612
+2.0575	1.39699
+2.0673	1.39903
+2.0773	1.39881
+2.0873	1.39687
+2.0975	1.39696
+2.1077	1.39914
+2.1180	1.39967
+2.1285	1.39987
+2.1390	1.40233
+2.1496	1.40370
+2.1568	1.40326
+2.1676	1.40133
+2.1785	1.39864
+2.1896	1.39930
+2.1970	1.40095
+2.2082	1.40247
+2.2196	1.40440
+2.2272	1.40536
+2.2387	1.40374
+2.2465	1.40253
+2.2582	1.40081
+2.2661	1.40030
+2.2781	1.40186
+2.2861	1.40307
+2.2983	1.40251
+2.3065	1.40157
+2.3188	1.40147
+2.3272	1.40085
+2.3398	1.40004
+2.3482	1.40030
+2.3568	1.40021
+2.3697	1.40048
+2.3784	1.40061
+2.3872	1.40076
+2.3960	1.40089
+2.4093	1.40133
+2.4183	1.40200
+2.4274	1.40179
+2.4365	1.40081
+2.4457	1.40033
+2.4596	1.40147
+2.4690	1.40204
+2.4784	1.40209
+2.4880	1.40105
+2.4975	1.40027
+2.5072	1.40027
+2.5169	1.39981
+2.5268	1.39979
+2.5366	1.40073
+2.5466	1.40055
+2.5567	1.39877
+2.5668	1.39652
+2.5770	1.39631
+2.5873	1.39760
+2.5976	1.39832
+2.6081	1.39925
+2.6186	1.40041
+2.6293	1.39945
+2.6400	1.39704
+2.6454	1.39654
+2.6562	1.39786
+2.6671	1.39947
+2.6782	1.39874
+2.6893	1.39730
+2.6949	1.39750
+2.7061	1.39903
+2.7175	1.39970
+2.7289	1.39986
+2.7347	1.40003
+2.7462	1.40070
+2.7579	1.39972
+2.7697	1.39755
+2.7756	1.39727
+2.7876	1.39815
+2.7996	1.39933
+2.8057	1.39987
+2.8179	1.39966
+2.8240	1.39890
+2.8364	1.39710
+2.8489	1.39627
+2.8551	1.39635
+2.8678	1.39733
+2.8741	1.39768
+2.8869	1.39797
+2.8998	1.39849
+2.9063	1.39860
+2.9194	1.39795
+2.9260	1.39761
+2.9393	1.39743
+2.9460	1.39741
+2.9594	1.39716
+2.9662	1.39669
+2.9798	1.39594
+2.9867	1.39628
+2.9936	1.39706
+3.0075	1.39749
+3.0145	1.39719
+3.0286	1.39643
+3.0357	1.39647
+3.0499	1.39718
+3.0571	1.39727
+3.0644	1.39740
+3.0789	1.39742
+3.0863	1.39713
+3.0936	1.39680
+3.1085	1.39584
+3.1159	1.39496
+3.1234	1.39379
+3.1386	1.39324
+3.1462	1.39379
+3.1538	1.39429
+3.1693	1.39476
+3.1770	1.39433
+3.1848	1.39372
+3.1927	1.39301
+3.2085	1.39247
+3.2164	1.39263
+3.2244	1.39267
+3.2325	1.39253
+3.2487	1.39146
+3.2569	1.39121
+3.2651	1.39100
+3.2733	1.39072
+3.2899	1.39062
+3.2983	1.38988
+3.3067	1.38851
+3.3152	1.38681
+3.3237	1.38496
+3.3322	1.38255
+3.3494	1.37358
+3.3581	1.37134
+3.3668	1.37853
+3.3756	1.39535
+3.3844	1.41081
+3.3933	1.41582
+3.4022	1.41294
+3.4111	1.40856
+3.4292	1.40376
+3.4383	1.40340
+3.4474	1.40403
+3.4566	1.40436
+3.4658	1.40409
+3.4751	1.40326
+3.4845	1.40221
+3.4939	1.40125
+3.5033	1.40063
+3.5128	1.40039
+3.5224	1.40016
+3.5320	1.39951
+3.5416	1.39853
+3.5513	1.39792
+3.5611	1.39841
+3.5709	1.39924
+3.5807	1.39928
+3.5907	1.39851
+3.6006	1.39810
+3.6107	1.39884
+3.6207	1.39913
+3.6309	1.39830
+3.6411	1.39753
+3.6513	1.39735
+3.6617	1.39690
+3.6720	1.39586
+3.6825	1.39530
+3.6930	1.39548
+3.7035	1.39588
+3.7141	1.39610
+3.7248	1.39566
+3.7355	1.39495
+3.7463	1.39450
+3.7572	1.39453
+3.7681	1.39495
+3.7791	1.39499
+3.7901	1.39450
+3.8013	1.39397
+3.8124	1.39377
+3.8237	1.39349
+3.8350	1.39339
+3.8464	1.39371
+3.8578	1.39355
+3.8693	1.39276
+3.8809	1.39219
+3.8926	1.39184
+3.9043	1.39149
+3.9161	1.39113
+3.9280	1.39115
+3.9399	1.39175
+3.9519	1.39239
+3.9640	1.39221
+3.9762	1.39149
+3.9884	1.39120
+4.0007	1.39128
+4.0131	1.39093
+4.0255	1.39025
+4.0381	1.38993
+4.0507	1.38997
+4.0634	1.39025
+4.0762	1.39022
+4.0890	1.38974
+4.1020	1.38930
+4.1150	1.38929
+4.1281	1.38936
+4.1413	1.38912
+4.1546	1.38879
+4.1679	1.38875
+4.1814	1.38874
+4.1949	1.38886
+4.2085	1.38916
+4.2222	1.38994
+4.2360	1.38929
+4.2499	1.38714
+4.2639	1.38749
+4.2780	1.38893
+4.2921	1.38864
+4.3064	1.38701
+4.3208	1.38723
+4.3352	1.38775
+4.3498	1.38746
+4.3644	1.38661
+4.3791	1.38593
+4.3940	1.38561
+4.4089	1.38551
+4.4240	1.38542
+4.4391	1.38535
+4.4544	1.38523
+4.4697	1.38444
+4.4852	1.38377
+4.5008	1.38368
+4.5165	1.38353
+4.5323	1.38336
+4.5482	1.38332
+4.5642	1.38294
+4.5803	1.38233
+4.5965	1.38243
+4.6129	1.38292
+4.6294	1.38280
+4.6460	1.38261
+4.6627	1.38250
+4.6795	1.38226
+4.6965	1.38173
+4.7136	1.38121
+4.7308	1.38039
+4.7481	1.37976
+4.7655	1.37976
+4.7831	1.37999
+4.8008	1.37979
+4.8187	1.37925
+4.8367	1.37867
+4.8548	1.37801
+4.8730	1.37747
+4.8914	1.37730
+4.9100	1.37710
+4.9286	1.37647
+4.9474	1.37595
+4.9664	1.37552
+4.9855	1.37505
+5.0047	1.37520
+5.0241	1.37538
+5.0437	1.37499
+5.0634	1.37463
+5.0832	1.37450
+5.1033	1.37445
+5.1234	1.37467
+5.1438	1.37432
+5.1643	1.37384
+5.1849	1.37399
+5.2057	1.37398
+5.2267	1.37307
+5.2479	1.37199
+5.2692	1.37158
+5.2907	1.37117
+5.3124	1.37023
+5.3343	1.36935
+5.3563	1.36952
+5.3785	1.36938
+5.4009	1.36887
+5.4235	1.36838
+5.4463	1.36772
+5.4693	1.36709
+5.4925	1.36624
+5.5159	1.36549
+5.5394	1.36500
+5.5632	1.36410
+5.5872	1.36347
+5.6114	1.36311
+5.6358	1.36284
+5.6604	1.36302
+5.6852	1.36229
+5.7103	1.36127
+5.7355	1.36110
+5.7610	1.36045
+5.7867	1.35957
+5.8127	1.35862
+5.8389	1.35789
+5.8653	1.35815
+5.8919	1.35778
+5.9188	1.35573
+5.9460	1.35404
+5.9734	1.35385
+6.0011	1.35358
+6.0290	1.35284
+6.0571	1.35203
+6.0856	1.35100
+6.1143	1.35000
+6.1433	1.34880
+6.1725	1.34723
+6.2020	1.34606
+6.2319	1.34605
+6.2620	1.34595
+6.2924	1.34520
+6.3231	1.34372
+6.3541	1.34273
+6.3854	1.34234
+6.4170	1.34110
+6.4489	1.33891
+6.4811	1.33700
+6.5137	1.33512
+6.5466	1.33371
+6.5798	1.33250
+6.6134	1.33083
+6.6473	1.32824
+6.6816	1.32662
+6.7162	1.32489
+6.7512	1.32235
+6.7865	1.31995
+6.8223	1.31759
+6.8583	1.31451
+6.8948	1.31220
+6.9317	1.31140
+6.9690	1.30858
+7.0066	1.30449
+7.0447	1.30322
+7.0832	1.30571
+7.1221	1.30963
+7.1615	1.31189
+7.2013	1.30840
+7.2415	1.30095
+7.2822	1.29452
+7.3233	1.28974
+7.3649	1.28464
+7.4070	1.27893
+7.4496	1.27215
+7.4926	1.26451
+7.5362	1.25588
+7.5803	1.24676
+7.6249	1.23572
+7.6700	1.22274
+7.7156	1.20633
+7.7618	1.18428
+7.8086	1.13732
+7.8559	1.06690
+7.9038	1.13149
+7.9523	1.31562
+8.0014	1.42755
+8.0511	1.37931
+8.1014	1.31645
+8.1524	1.27762
+8.2040	1.24968
+8.2562	1.22830
+8.3091	1.21345
+8.3628	1.20275
+8.4171	1.19147
+8.4721	1.17666
+8.5278	1.15959
+8.5843	1.14061
+8.6415	1.11362
+8.6995	1.07813
+8.7583	1.03538
+8.8179	0.98636
+8.8783	0.93665
+8.9395	0.90189
+9.0016	0.91161
+9.0645	0.99110
+9.1284	1.12930
+9.1931	1.28201
+9.2588	1.40749
+9.3254	1.49129
+9.3930	1.54005
+9.4615	1.56597
+9.5311	1.57798
+9.6017	1.58843
+9.6733	1.61404
+9.7461	1.67584
+9.8199	1.77617
+9.8949	1.85190
+9.9710	1.84298
+10.048	1.78062
+10.127	1.71731
+10.207	1.66568
+10.288	1.62340
+10.370	1.58894
+10.453	1.55986
+10.538	1.53370
+10.625	1.50932
+10.713	1.48704
+10.802	1.46405
+10.893	1.43983
+10.985	1.41698
+11.079	1.39386
+11.174	1.36267
+11.272	1.32503
+11.370	1.29663
+11.471	1.30210
+11.573	1.32967
+11.678	1.33658
+11.784	1.33938
+11.892	1.34524
+12.002	1.29807
+12.114	1.24035
+12.229	1.26473
+12.345	1.38610
+12.464	1.61522
+12.585	1.86651
+12.708	1.96260
+12.834	1.89309
+12.962	1.79395
+13.093	1.73241
+13.227	1.71030
+13.363	1.69424
+13.502	1.66425
+13.645	1.63582
+13.790	1.61199
+13.938	1.59044
+14.089	1.57874
+14.244	1.57910
+14.403	1.58302
+14.564	1.58508
+14.730	1.58381
+14.899	1.57887
+15.072	1.57711
+15.250	1.57350
+15.431	1.56094
+15.617	1.54848
+15.808	1.54164
+16.003	1.53726
+16.203	1.52800
+16.408	1.51659
+16.618	1.51254
+16.834	1.51107
+17.056	1.50150
+17.283	1.49402
+17.517	1.49619
+17.757	1.50036
+18.003	1.49790
+18.257	1.49605
+18.518	1.49709
+18.786	1.49362
+19.062	1.48785
+19.347	1.48798
+19.640	1.48985
+19.942	1.49006
 
 wl	k
-0.40	1.40E-06
-0.41	1.38E-06
-0.42	1.38E-06
-0.43	1.40E-06
-0.44	1.40E-06
-0.45	1.41E-06
-0.46	1.42E-06
-0.47	1.43E-06
-0.48	1.43E-06
-0.49	1.46E-06
-0.50	1.46E-06
-0.51	1.47E-06
-0.52	1.48E-06
-0.53	1.48E-06
-0.54	1.49E-06
-0.55	1.51E-06
-0.56	1.50E-06
-0.57	1.53E-06
-0.58	1.55E-06
-0.59	1.54E-06
-0.60	1.57E-06
-0.61	1.57E-06
-0.62	1.57E-06
-0.63	1.58E-06
-0.64	1.58E-06
-0.65	1.59E-06
-0.66	1.58E-06
-0.67	1.59E-06
-0.68	1.60E-06
-0.69	1.60E-06
-0.70	1.61E-06
-0.71	1.62E-06
-0.72	1.63E-06
-0.73	1.64E-06
-0.74	1.69E-06
-0.75	1.66E-06
-0.76	1.65E-06
-0.77	1.67E-06
-0.78	1.67E-06
-0.79	1.67E-06
-0.80	1.69E-06
-0.81	1.70E-06
-0.82	1.72E-06
-0.83	1.73E-06
-0.84	1.62E-06
-0.85	1.76E-06
-0.86	1.70E-06
-0.87	2.22E-06
-0.88	2.34E-06
-0.89	2.39E-06
-0.90	2.32E-06
-0.91	2.31E-06
-0.92	2.79E-06
-0.93	2.60E-06
-0.94	2.58E-06
-0.95	2.32E-06
-0.96	2.78E-06
-0.97	2.63E-06
-0.98	2.45E-06
-0.99	2.56E-06
-1.00	2.88E-06
-1.01	2.26E-06
-1.02	2.61E-06
-1.03	2.54E-06
-1.04	2.44E-06
-1.05	2.38E-06
-1.06	2.47E-06
-1.07	2.52E-06
-1.08	2.57E-06
-1.09	2.76E-06
-1.10	2.70E-06
-1.11	2.40E-06
-1.12	2.66E-06
-1.13	2.85E-06
-1.14	3.42E-06
-1.15	4.81E-06
-1.16	4.84E-06
-1.17	5.34E-06
-1.18	1.10E-05
-1.19	1.46E-05
-1.20	4.75E-06
-1.21	3.87E-06
-1.22	3.70E-06
-1.23	3.00E-06
-1.24	2.89E-06
-1.25	3.00E-06
-1.26	3.04E-06
-1.27	3.16E-06
-1.28	3.19E-06
-1.29	3.11E-06
-1.30	3.14E-06
-1.31	3.11E-06
-1.32	3.12E-06
-1.33	3.25E-06
-1.34	3.54E-06
-1.35	4.13E-06
-1.36	5.43E-06
-1.37	9.45E-06
-1.38	9.79E-06
-1.39	1.14E-05
-1.40	1.42E-05
-1.41	1.24E-05
-1.42	7.87E-06
-1.43	6.99E-06
-1.44	6.05E-06
-1.45	5.11E-06
-1.46	5.12E-06
-1.47	5.10E-06
-1.48	3.03E-06
-1.49	5.34E-06
-1.50	5.84E-06
-1.51	6.90E-06
-1.52	7.80E-06
-1.53	7.85E-06
-1.54	7.98E-06
-1.55	7.83E-06
-1.56	7.28E-06
-1.57	5.95E-06
-1.58	5.23E-06
-1.59	5.01E-06
-1.60	5.04E-06
-1.61	5.17E-06
-1.62	5.53E-06
-1.63	6.14E-06
-1.64	7.67E-06
-1.65	9.77E-06
-1.66	1.02E-05
-1.67	1.08E-05
-1.68	1.81E-05
-1.69	7.25E-05
-1.70	5.97E-05
-1.71	3.22E-05
-1.72	2.16E-05
-1.73	3.96E-05
-1.74	7.03E-05
-1.75	7.41E-05
-1.76	3.46E-05
-1.77	2.28E-05
-1.78	2.78E-05
-1.79	2.36E-05
-1.80	1.26E-05
-1.81	9.26E-06
-1.82	8.87E-06
-1.83	1.39E-05
-1.84	1.97E-05
-1.85	2.29E-05
-1.86	1.35E-05
-1.87	8.20E-06
-1.88	7.43E-06
-1.89	7.87E-06
-1.90	8.45E-06
-1.91	8.83E-06
-1.92	9.57E-06
-1.93	1.01E-05
-1.94	1.08E-05
-1.95	1.11E-05
-1.96	1.25E-05
-1.97	1.45E-05
-1.98	1.42E-05
-1.99	1.40E-05
-2.00	1.31E-05
-2.0097	1.20E-05
-2.0190	1.23E-05
-2.0285	1.03E-05
-2.0381	1.01E-05
-2.0478	9.51E-07
-2.0575	8.81E-06
-2.0673	8.39E-06
-2.0773	8.30E-06
-2.0873	8.51E-06
-2.0975	8.99E-06
-2.1077	9.64E-06
-2.1180	1.05E-05
-2.1285	1.20E-05
-2.1390	1.42E-05
-2.1496	1.73E-05
-2.1568	1.96E-05
-2.1676	2.38E-05
-2.1785	2.98E-05
-2.1896	3.83E-05
-2.1970	4.46E-05
-2.2082	5.46E-05
-2.2196	6.48E-05
-2.2272	7.02E-05
-2.2387	7.63E-05
-2.2465	8.03E-05
-2.2582	8.48E-05
-2.2661	8.59E-05
-2.2781	8.61E-05
-2.2861	8.63E-05
-2.2983	8.77E-05
-2.3065	8.90E-05
-2.3188	8.91E-05
-2.3272	8.67E-05
-2.3398	8.38E-05
-2.3482	8.49E-05
-2.3568	8.79E-05
-2.3697	9.29E-05
-2.3784	9.52E-05
-2.3872	9.55E-05
-2.3960	9.51E-05
-2.4093	9.40E-05
-2.4183	9.05E-05
-2.4274	8.79E-05
-2.4365	9.00E-05
-2.4457	9.57E-05
-2.4596	9.98E-05
-2.4690	9.93E-05
-2.4784	9.86E-05
-2.4880	9.80E-05
-2.4975	9.74E-05
-2.5072	9.51E-05
-2.5169	9.13E-05
-2.5268	8.86E-05
-2.5366	8.74E-05
-2.5466	8.69E-05
-2.5567	8.80E-05
-2.5668	9.20E-05
-2.5770	9.70E-05
-2.5873	1.00E-04
-2.5976	1.03E-04
-2.6081	1.06E-04
-2.6186	1.08E-04
-2.6293	1.08E-04
-2.6400	1.08E-04
-2.6454	1.08E-04
-2.6562	1.09E-04
-2.6671	1.09E-04
-2.6782	1.08E-04
-2.6893	1.09E-04
-2.6949	1.10E-04
-2.7061	1.11E-04
-2.7175	1.12E-04
-2.7289	1.14E-04
-2.7347	1.14E-04
-2.7462	1.14E-04
-2.7579	1.13E-04
-2.7697	1.15E-04
-2.7756	1.16E-04
-2.7876	1.16E-04
-2.7996	1.17E-04
-2.8057	1.18E-04
-2.8179	1.19E-04
-2.8240	1.19E-04
-2.8364	1.20E-04
-2.8489	1.21E-04
-2.8551	1.21E-04
-2.8678	1.22E-04
-2.8741	1.21E-04
-2.8869	1.21E-04
-2.8998	1.21E-04
-2.9063	1.21E-04
-2.9194	1.22E-04
-2.9260	1.22E-04
-2.9393	1.22E-04
-2.9460	1.22E-04
-2.9594	1.24E-04
-2.9662	1.24E-04
-2.9798	1.23E-04
-2.9867	1.24E-04
-2.9936	1.26E-04
-3.0075	1.26E-04
-3.0145	1.25E-04
-3.0286	1.25E-04
-3.0357	1.27E-04
-3.0499	1.30E-04
-3.0571	1.31E-04
-3.0644	1.32E-04
-3.0789	1.33E-04
-3.0863	1.33E-04
-3.0936	1.33E-04
-3.1085	1.35E-04
-3.1159	1.37E-04
-3.1234	1.39E-04
-3.1386	1.40E-04
-3.1462	1.40E-04
-3.1538	1.40E-04
-3.1693	1.39E-04
-3.1770	1.37E-04
-3.1848	1.37E-04
-3.1927	1.38E-04
-3.2085	1.42E-04
-3.2164	1.45E-04
-3.2244	1.46E-04
-3.2325	1.45E-04
-3.2487	1.46E-04
-3.2569	1.47E-04
-3.2651	1.49E-04
-3.2733	1.51E-04
-3.2899	1.51E-04
-3.2983	1.48E-04
-3.3067	1.47E-04
-3.3152	1.48E-04
-3.3237	1.50E-04
-3.3322	1.00E-03
-3.3494	8.07E-03
-3.3581	1.86E-02
-3.3668	3.19E-02
-3.3756	3.83E-02
+0.40	9.29E-07
+0.41	9.08E-07
+0.42	9.01E-07
+0.43	9.07E-07
+0.44	9.05E-07
+0.45	8.99E-07
+0.46	9.02E-07
+0.47	9.02E-07
+0.48	9.08E-07
+0.49	9.17E-07
+0.50	9.06E-07
+0.51	9.21E-07
+0.52	9.12E-07
+0.53	9.04E-07
+0.54	9.09E-07
+0.55	9.16E-07
+0.56	9.01E-07
+0.57	9.12E-07
+0.58	9.14E-07
+0.59	9.02E-07
+0.60	9.32E-07
+0.61	9.04E-07
+0.62	9.10E-07
+0.63	8.99E-07
+0.64	9.04E-07
+0.65	8.95E-07
+0.66	8.72E-07
+0.67	8.66E-07
+0.68	8.71E-07
+0.69	8.48E-07
+0.70	8.66E-07
+0.71	8.57E-07
+0.72	8.61E-07
+0.73	8.66E-07
+0.74	9.06E-07
+0.75	8.60E-07
+0.76	8.40E-07
+0.77	8.43E-07
+0.78	8.32E-07
+0.79	8.13E-07
+0.80	8.06E-07
+0.81	8.32E-07
+0.82	6.80E-07
+0.83	7.14E-07
+0.84	7.12E-07
+0.85	8.60E-07
+0.86	7.98E-07
+0.87	1.32E-06
+0.88	1.42E-06
+0.89	1.40E-06
+0.90	1.43E-06
+0.91	1.53E-06
+0.92	1.47E-06
+0.93	1.76E-06
+0.94	1.52E-06
+0.95	1.40E-06
+0.96	1.75E-06
+0.97	1.79E-06
+0.98	1.53E-06
+0.99	1.37E-06
+1.00	1.80E-06
+1.01	1.55E-06
+1.02	1.78E-06
+1.03	1.69E-06
+1.04	1.55E-06
+1.05	1.45E-06
+1.06	1.51E-06
+1.07	1.60E-06
+1.08	1.70E-06
+1.09	1.98E-06
+1.10	1.77E-06
+1.11	1.58E-06
+1.12	1.72E-06
+1.13	1.98E-06
+1.14	2.62E-06
+1.15	4.44E-06
+1.16	4.43E-06
+1.17	5.13E-06
+1.18	1.22E-05
+1.19	1.68E-05
+1.20	4.17E-06
+1.21	3.01E-06
+1.22	2.83E-06
+1.23	1.95E-06
+1.24	1.89E-06
+1.25	1.95E-06
+1.26	1.96E-06
+1.27	2.30E-06
+1.28	2.11E-06
+1.29	2.07E-06
+1.30	2.11E-06
+1.31	2.02E-06
+1.32	2.03E-06
+1.33	2.16E-06
+1.34	2.55E-06
+1.35	3.20E-06
+1.36	4.93E-06
+1.37	9.94E-06
+1.38	1.01E-05
+1.39	1.22E-05
+1.40	1.59E-05
+1.41	1.35E-05
+1.42	7.84E-06
+1.43	6.69E-06
+1.44	5.54E-06
+1.45	4.38E-06
+1.46	4.41E-06
+1.47	4.43E-06
+1.48	2.54E-06
+1.49	4.64E-06
+1.50	5.18E-06
+1.51	6.58E-06
+1.52	7.61E-06
+1.53	7.71E-06
+1.54	7.77E-06
+1.55	7.85E-06
+1.56	6.97E-06
+1.57	5.25E-06
+1.58	4.25E-06
+1.59	3.94E-06
+1.60	3.99E-06
+1.61	4.17E-06
+1.62	4.68E-06
+1.63	5.52E-06
+1.64	7.63E-06
+1.65	1.16E-05
+1.66	1.09E-05
+1.67	1.13E-05
+1.68	2.04E-05
+1.69	8.84E-05
+1.70	7.16E-05
+1.71	3.78E-05
+1.72	2.48E-05
+1.73	4.66E-05
+1.74	8.54E-05
+1.75	9.00E-05
+1.76	4.01E-05
+1.77	2.54E-05
+1.78	2.54E-05
+1.79	2.63E-05
+1.80	1.30E-05
+1.81	8.95E-06
+1.82	8.47E-06
+1.83	1.50E-05
+1.84	2.21E-05
+1.85	2.61E-05
+1.86	1.39E-05
+1.87	7.61E-06
+1.88	6.69E-06
+1.89	7.04E-06
+1.90	7.68E-06
+1.91	8.38E-06
+1.92	9.07E-06
+1.93	9.60E-06
+1.94	1.03E-05
+1.95	1.02E-05
+1.96	1.22E-05
+1.97	1.48E-05
+1.98	1.46E-05
+1.99	1.45E-05
+2.00	1.36E-05
+2.0097	1.57E-05
+2.0190	1.41E-05
+2.0285	1.12E-05
+2.0381	1.09E-05
+2.0478	1.03E-05
+2.0575	9.48E-06
+2.0673	8.99E-06
+2.0773	8.88E-06
+2.0873	9.15E-06
+2.0975	9.70E-06
+2.1077	1.04E-05
+2.1180	1.14E-05
+2.1285	1.32E-05
+2.1390	1.63E-05
+2.1496	2.04E-05
+2.1568	2.28E-05
+2.1676	2.64E-05
+2.1785	3.25E-05
+2.1896	4.14E-05
+2.1970	4.82E-05
+2.2082	5.87E-05
+2.2196	6.87E-05
+2.2272	7.38E-05
+2.2387	7.99E-05
+2.2465	8.41E-05
+2.2582	8.90E-05
+2.2661	9.09E-05
+2.2781	9.16E-05
+2.2861	9.17E-05
+2.2983	9.27E-05
+2.3065	9.36E-05
+2.3188	9.28E-05
+2.3272	9.02E-05
+2.3398	8.73E-05
+2.3482	8.88E-05
+2.3568	9.24E-05
+2.3697	9.74E-05
+2.3784	9.91E-05
+2.3872	9.95E-05
+2.3960	9.95E-05
+2.4093	9.74E-05
+2.4183	9.34E-05
+2.4274	9.11E-05
+2.4365	9.33E-05
+2.4457	9.86E-05
+2.4596	1.03E-04
+2.4690	1.04E-04
+2.4784	1.05E-04
+2.4880	1.05E-04
+2.4975	1.03E-04
+2.5072	1.00E-04
+2.5169	9.67E-05
+2.5268	9.33E-05
+2.5366	9.10E-05
+2.5466	9.01E-05
+2.5567	9.13E-05
+2.5668	9.45E-05
+2.5770	9.85E-05
+2.5873	1.03E-04
+2.5976	1.08E-04
+2.6081	1.10E-04
+2.6186	1.12E-04
+2.6293	1.13E-04
+2.6400	1.13E-04
+2.6454	1.13E-04
+2.6562	1.13E-04
+2.6671	1.13E-04
+2.6782	1.14E-04
+2.6893	1.15E-04
+2.6949	1.16E-04
+2.7061	1.17E-04
+2.7175	1.19E-04
+2.7289	1.20E-04
+2.7347	1.20E-04
+2.7462	1.20E-04
+2.7579	1.19E-04
+2.7697	1.19E-04
+2.7756	1.20E-04
+2.7876	1.22E-04
+2.7996	1.23E-04
+2.8057	1.24E-04
+2.8179	1.25E-04
+2.8240	1.25E-04
+2.8364	1.25E-04
+2.8489	1.25E-04
+2.8551	1.25E-04
+2.8678	1.27E-04
+2.8741	1.27E-04
+2.8869	1.28E-04
+2.8998	1.28E-04
+2.9063	1.29E-04
+2.9194	1.28E-04
+2.9260	1.28E-04
+2.9393	1.28E-04
+2.9460	1.29E-04
+2.9594	1.29E-04
+2.9662	1.30E-04
+2.9798	1.31E-04
+2.9867	1.32E-04
+2.9936	1.33E-04
+3.0075	1.32E-04
+3.0145	1.32E-04
+3.0286	1.33E-04
+3.0357	1.34E-04
+3.0499	1.36E-04
+3.0571	1.37E-04
+3.0644	1.37E-04
+3.0789	1.39E-04
+3.0863	1.40E-04
+3.0936	1.40E-04
+3.1085	1.42E-04
+3.1159	1.44E-04
+3.1234	1.44E-04
+3.1386	1.43E-04
+3.1462	1.44E-04
+3.1538	1.46E-04
+3.1693	1.47E-04
+3.1770	1.46E-04
+3.1848	1.47E-04
+3.1927	1.47E-04
+3.2085	1.49E-04
+3.2164	1.50E-04
+3.2244	1.50E-04
+3.2325	1.50E-04
+3.2487	1.52E-04
+3.2569	1.54E-04
+3.2651	1.55E-04
+3.2733	1.57E-04
+3.2899	1.58E-04
+3.2983	1.57E-04
+3.3067	1.57E-04
+3.3152	1.57E-04
+3.3237	1.57E-04
+3.3322	1.50E-03
+3.3494	8.33E-03
+3.3581	1.87E-02
+3.3668	3.17E-02
+3.3756	3.79E-02
 3.3844	3.13E-02
-3.3933	1.83E-02
-3.4022	8.36E-03
-3.4111	3.55E-03
-3.4292	3.33E-03
-3.4383	4.63E-03
-3.4474	4.56E-03
-3.4566	3.00E-03
-3.4658	1.45E-03
-3.4751	1.05E-03
-3.4845	1.63E-04
-3.4939	1.63E-04
-3.5033	1.64E-04
-3.5128	1.65E-04
-3.5224	1.67E-04
-3.5320	1.66E-04
-3.5416	1.65E-04
-3.5513	1.67E-04
-3.5611	1.70E-04
-3.5709	1.69E-04
-3.5807	1.70E-04
-3.5907	1.71E-04
-3.6006	1.70E-04
-3.6107	1.69E-04
-3.6207	1.70E-04
-3.6309	1.71E-04
-3.6411	1.73E-04
-3.6513	1.76E-04
-3.6617	1.75E-04
-3.6720	1.74E-04
-3.6825	1.75E-04
-3.6930	1.78E-04
-3.7035	1.80E-04
-3.7141	1.81E-04
-3.7248	1.82E-04
-3.7355	1.81E-04
-3.7463	1.78E-04
-3.7572	1.79E-04
-3.7681	1.82E-04
-3.7791	1.70E-04
-3.7901	1.43E-04
-3.8013	1.19E-04
-3.8124	1.02E-04
-3.8237	8.86E-05
-3.8350	7.94E-05
-3.8464	7.28E-05
-3.8578	6.83E-05
-3.8693	6.53E-05
-3.8809	6.36E-05
-3.8926	6.30E-05
-3.9043	6.36E-05
-3.9161	6.56E-05
-3.9280	6.92E-05
-3.9399	7.42E-05
-3.9519	8.17E-05
-3.9640	9.42E-05
-3.9762	1.15E-04
-3.9884	1.46E-04
-4.0007	1.78E-04
-4.0131	1.65E-04
-4.0255	1.29E-04
-4.0381	1.02E-04
-4.0507	8.45E-05
-4.0634	7.40E-05
-4.0762	6.74E-05
-4.0890	6.32E-05
-4.1020	6.05E-05
-4.1150	5.90E-05
-4.1281	5.87E-05
-4.1413	5.93E-05
-4.1546	6.07E-05
-4.1679	6.28E-05
-4.1814	6.53E-05
-4.1949	6.77E-05
-4.2085	6.99E-05
-4.2222	7.22E-05
-4.2360	7.61E-05
-4.2499	8.06E-05
-4.2639	8.41E-05
-4.2780	8.87E-05
-4.2921	9.50E-05
-4.3064	1.03E-04
-4.3208	1.11E-04
-4.3352	1.20E-04
-4.3498	1.31E-04
-4.3644	1.45E-04
-4.3791	1.59E-04
-4.3940	1.73E-04
-4.4089	1.87E-04
-4.4240	1.98E-04
-4.4391	2.09E-04
-4.4544	2.19E-04
-4.4697	2.26E-04
-4.4852	2.26E-04
-4.5008	2.27E-04
-4.5165	2.24E-04
-4.5323	2.25E-04
-4.5482	2.33E-04
-4.5642	2.37E-04
-4.5803	2.36E-04
-4.5965	2.35E-04
-4.6129	2.33E-04
-4.6294	2.32E-04
-4.6460	2.34E-04
-4.6627	2.36E-04
-4.6795	2.34E-04
-4.6965	2.35E-04
-4.7136	2.40E-04
-4.7308	2.38E-04
-4.7481	2.39E-04
-4.7655	2.42E-04
-4.7831	2.45E-04
-4.8008	2.49E-04
-4.8187	2.49E-04
-4.8367	2.45E-04
-4.8548	2.45E-04
-4.8730	2.49E-04
-4.8914	2.47E-04
-4.9100	2.50E-04
-4.9286	2.53E-04
-4.9474	2.58E-04
-4.9664	2.61E-04
-4.9855	2.63E-04
-5.0047	2.58E-04
-5.0241	2.55E-04
-5.0437	2.64E-04
-5.0634	2.60E-04
-5.0832	2.58E-04
-5.1033	2.62E-04
-5.1234	2.63E-04
-5.1438	2.65E-04
-5.1643	2.60E-04
-5.1849	2.60E-04
-5.2057	2.67E-04
-5.2267	2.80E-04
-5.2479	2.72E-04
-5.2692	2.69E-04
-5.2907	2.76E-04
-5.3124	2.74E-04
-5.3343	2.71E-04
-5.3563	2.67E-04
-5.3785	2.82E-04
-5.4009	2.80E-04
-5.4235	2.81E-04
-5.4463	2.77E-04
-5.4693	2.77E-04
-5.4925	2.76E-04
-5.5159	2.82E-04
-5.5394	2.84E-04
-5.5632	2.79E-04
-5.5872	2.81E-04
-5.6114	2.77E-04
-5.6358	2.79E-04
-5.6604	2.91E-04
-5.6852	2.90E-04
-5.7103	2.79E-04
-5.7355	2.77E-04
-5.7610	2.82E-04
-5.7867	2.95E-04
-5.8127	2.91E-04
-5.8389	2.82E-04
-5.8653	2.75E-04
-5.8919	2.85E-04
-5.9188	2.94E-04
-5.9460	3.07E-04
-5.9734	3.03E-04
-6.0011	2.94E-04
-6.0290	5.02E-04
-6.0571	3.99E-04
-6.0856	2.38E-04
-6.1143	4.84E-04
-6.1433	1.09E-03
-6.1725	1.21E-03
-6.2020	9.51E-04
-6.2319	1.20E-03
-6.2620	1.10E-03
-6.2924	5.76E-04
-6.3231	3.39E-04
-6.3541	3.13E-04
-6.3854	3.94E-04
-6.4170	4.60E-04
-6.4489	7.13E-04
-6.4811	7.17E-04
-6.5137	4.91E-04
-6.5466	5.19E-04
-6.5798	2.60E-04
-6.6134	5.84E-04
-6.6473	7.76E-04
-6.6816	6.03E-04
-6.7162	5.89E-04
-6.7512	6.78E-04
-6.7865	1.42E-03
-6.8223	2.57E-03
-6.8583	4.35E-03
-6.8948	6.58E-03
-6.9317	7.04E-03
-6.9690	6.48E-03
-7.0066	8.61E-03
-7.0447	1.37E-02
-7.0832	1.80E-02
-7.1221	1.67E-02
-7.1615	1.21E-02
-7.2013	7.40E-03
-7.2415	4.83E-03
-7.2822	3.47E-03
-7.3233	3.07E-03
-7.3649	2.87E-03
-7.4070	3.11E-03
-7.4496	3.12E-03
-7.4926	2.59E-03
-7.5362	2.52E-03
-7.5803	3.96E-03
-7.6249	6.05E-03
-7.6700	7.77E-03
-7.7156	9.38E-03
-7.7618	1.47E-02
-7.8086	2.93E-02
-7.8559	1.71E-01
-7.9038	3.31E-01
-7.9523	3.32E-01
-8.0014	1.74E-01
-8.0511	5.40E-02
-8.1014	2.73E-02
-8.1524	2.43E-02
-8.2040	2.83E-02
-8.2562	3.56E-02
-8.3091	4.49E-02
-8.3628	5.22E-02
-8.4171	5.81E-02
-8.4721	6.32E-02
-8.5278	6.79E-02
-8.5843	7.41E-02
-8.6415	8.41E-02
-8.6995	1.02E-01
-8.7583	1.31E-01
-8.8179	1.75E-01
-8.8783	2.40E-01
-8.9395	3.37E-01
-9.0016	4.69E-01
-9.0645	6.02E-01
-9.1284	6.88E-01
-9.1931	7.05E-01
-9.2588	6.71E-01
-9.3254	6.16E-01
-9.3930	5.62E-01
-9.4615	5.19E-01
-9.5311	4.91E-01
-9.6017	4.83E-01
-9.6733	4.91E-01
-9.7461	4.96E-01
-9.8199	4.53E-01
+3.3933	1.86E-02
+3.4022	9.32E-03
+3.4111	5.23E-03
+3.4292	4.49E-03
+3.4383	5.43E-03
+3.4474	5.70E-03
+3.4566	4.68E-03
+3.4658	3.05E-03
+3.4751	1.64E-03
+3.4845	1.70E-04
+3.4939	1.71E-04
+3.5033	1.72E-04
+3.5128	1.73E-04
+3.5224	1.74E-04
+3.5320	1.74E-04
+3.5416	1.75E-04
+3.5513	1.77E-04
+3.5611	1.77E-04
+3.5709	1.76E-04
+3.5807	1.78E-04
+3.5907	1.80E-04
+3.6006	1.79E-04
+3.6107	1.78E-04
+3.6207	1.78E-04
+3.6309	1.81E-04
+3.6411	1.84E-04
+3.6513	1.85E-04
+3.6617	1.86E-04
+3.6720	1.84E-04
+3.6825	1.83E-04
+3.6930	1.84E-04
+3.7035	1.84E-04
+3.7141	1.83E-04
+3.7248	1.86E-04
+3.7355	1.88E-04
+3.7463	1.89E-04
+3.7572	1.90E-04
+3.7681	1.90E-04
+3.7791	1.76E-04
+3.7901	1.50E-04
+3.8013	1.26E-04
+3.8124	1.08E-04
+3.8237	9.43E-05
+3.8350	8.48E-05
+3.8464	7.81E-05
+3.8578	7.33E-05
+3.8693	7.01E-05
+3.8809	6.80E-05
+3.8926	6.71E-05
+3.9043	6.74E-05
+3.9161	6.93E-05
+3.9280	7.27E-05
+3.9399	7.78E-05
+3.9519	8.58E-05
+3.9640	9.92E-05
+3.9762	1.21E-04
+3.9884	1.53E-04
+4.0007	1.87E-04
+4.0131	1.74E-04
+4.0255	1.36E-04
+4.0381	1.08E-04
+4.0507	8.99E-05
+4.0634	7.87E-05
+4.0762	7.17E-05
+4.0890	6.73E-05
+4.1020	6.45E-05
+4.1150	6.30E-05
+4.1281	6.26E-05
+4.1413	6.32E-05
+4.1546	6.45E-05
+4.1679	6.64E-05
+4.1814	6.87E-05
+4.1949	7.11E-05
+4.2085	7.34E-05
+4.2222	7.60E-05
+4.2360	8.01E-05
+4.2499	8.44E-05
+4.2639	8.78E-05
+4.2780	9.23E-05
+4.2921	9.86E-05
+4.3064	1.06E-04
+4.3208	1.14E-04
+4.3352	1.24E-04
+4.3498	1.34E-04
+4.3644	1.47E-04
+4.3791	1.62E-04
+4.3940	1.75E-04
+4.4089	1.86E-04
+4.4240	1.93E-04
+4.4391	1.99E-04
+4.4544	2.05E-04
+4.4697	2.13E-04
+4.4852	2.22E-04
+4.5008	2.30E-04
+4.5165	2.35E-04
+4.5323	2.39E-04
+4.5482	2.42E-04
+4.5642	2.44E-04
+4.5803	2.44E-04
+4.5965	2.43E-04
+4.6129	2.40E-04
+4.6294	2.41E-04
+4.6460	2.46E-04
+4.6627	2.50E-04
+4.6795	2.48E-04
+4.6965	2.47E-04
+4.7136	2.47E-04
+4.7308	2.49E-04
+4.7481	2.50E-04
+4.7655	2.54E-04
+4.7831	2.57E-04
+4.8008	2.59E-04
+4.8187	2.60E-04
+4.8367	2.59E-04
+4.8548	2.60E-04
+4.8730	2.61E-04
+4.8914	2.59E-04
+4.9100	2.62E-04
+4.9286	2.64E-04
+4.9474	2.65E-04
+4.9664	2.68E-04
+4.9855	2.68E-04
+5.0047	2.67E-04
+5.0241	2.68E-04
+5.0437	2.73E-04
+5.0634	2.70E-04
+5.0832	2.74E-04
+5.1033	2.75E-04
+5.1234	2.76E-04
+5.1438	2.80E-04
+5.1643	2.74E-04
+5.1849	2.72E-04
+5.2057	2.79E-04
+5.2267	2.83E-04
+5.2479	2.81E-04
+5.2692	2.81E-04
+5.2907	2.87E-04
+5.3124	2.93E-04
+5.3343	2.90E-04
+5.3563	2.80E-04
+5.3785	2.87E-04
+5.4009	2.93E-04
+5.4235	2.96E-04
+5.4463	2.93E-04
+5.4693	2.94E-04
+5.4925	2.94E-04
+5.5159	2.97E-04
+5.5394	2.97E-04
+5.5632	2.92E-04
+5.5872	2.93E-04
+5.6114	2.97E-04
+5.6358	2.97E-04
+5.6604	3.02E-04
+5.6852	3.00E-04
+5.7103	2.99E-04
+5.7355	2.96E-04
+5.7610	2.96E-04
+5.7867	2.99E-04
+5.8127	2.98E-04
+5.8389	2.96E-04
+5.8653	2.96E-04
+5.8919	2.97E-04
+5.9188	2.98E-04
+5.9460	3.04E-04
+5.9734	3.09E-04
+6.0011	3.05E-04
+6.0290	1.60E-03
+6.0571	1.55E-03
+6.0856	2.08E-03
+6.1143	1.90E-03
+6.1433	1.35E-03
+6.1725	1.30E-03
+6.2020	2.06E-03
+6.2319	2.28E-03
+6.2620	1.58E-03
+6.2924	1.60E-03
+6.3231	1.86E-03
+6.3541	2.10E-03
+6.3854	1.76E-03
+6.4170	1.44E-03
+6.4489	1.52E-03
+6.4811	1.54E-03
+6.5137	1.40E-03
+6.5466	1.85E-03
+6.5798	1.88E-03
+6.6134	1.75E-03
+6.6473	1.60E-03
+6.6816	1.59E-03
+6.7162	1.63E-03
+6.7512	2.29E-03
+6.7865	3.02E-03
+6.8223	3.16E-03
+6.8583	4.75E-03
+6.8948	8.02E-03
+6.9317	8.69E-03
+6.9690	7.82E-03
+7.0066	9.06E-03
+7.0447	1.39E-02
+7.0832	1.84E-02
+7.1221	1.81E-02
+7.1615	1.40E-02
+7.2013	8.61E-03
+7.2415	4.65E-03
+7.2822	3.90E-03
+7.3233	3.87E-03
+7.3649	3.74E-03
+7.4070	3.56E-03
+7.4496	2.47E-03
+7.4926	2.53E-03
+7.5362	3.80E-03
+7.5803	5.70E-03
+7.6249	6.82E-03
+7.6700	7.81E-03
+7.7156	9.63E-03
+7.7618	1.48E-02
+7.8086	3.02E-02
+7.8559	1.79E-01
+7.9038	3.34E-01
+7.9523	3.31E-01
+8.0014	1.75E-01
+8.0511	5.71E-02
+8.1014	3.05E-02
+8.1524	2.87E-02
+8.2040	3.47E-02
+8.2562	4.22E-02
+8.3091	5.07E-02
+8.3628	5.97E-02
+8.4171	6.70E-02
+8.4721	7.03E-02
+8.5278	7.41E-02
+8.5843	8.07E-02
+8.6415	8.99E-02
+8.6995	1.09E-01
+8.7583	1.39E-01
+8.8179	1.82E-01
+8.8783	2.49E-01
+8.9395	3.50E-01
+9.0016	4.87E-01
+9.0645	6.25E-01
+9.1284	7.16E-01
+9.1931	7.36E-01
+9.2588	7.00E-01
+9.3254	6.37E-01
+9.3930	5.75E-01
+9.4615	5.24E-01
+9.5311	4.88E-01
+9.6017	4.74E-01
+9.6733	4.79E-01
+9.7461	4.87E-01
+9.8199	4.49E-01
 9.8949	3.34E-01
-9.9710	1.94E-01
-10.048	1.06E-01
-10.127	6.35E-02
-10.207	4.21E-02
-10.288	3.12E-02
-10.370	2.52E-02
-10.453	2.08E-02
-10.538	1.83E-02
-10.625	1.81E-02
-10.713	1.82E-02
-10.802	2.14E-02
-10.893	3.13E-02
-10.985	4.01E-02
-11.079	4.13E-02
-11.174	4.66E-02
-11.272	6.75E-02
-11.370	1.12E-01
-11.471	1.67E-01
-11.573	1.93E-01
-11.678	2.02E-01
-11.784	2.22E-01
-11.892	2.27E-01
-12.002	2.44E-01
-12.114	3.43E-01
-12.229	4.89E-01
-12.345	6.40E-01
-12.464	7.21E-01
-12.585	6.09E-01
-12.708	3.61E-01
-12.834	1.74E-01
-12.962	1.02E-01
-13.093	8.96E-02
-13.227	8.48E-02
-13.363	6.27E-02
-13.502	4.22E-02
-13.645	3.49E-02
-13.790	3.34E-02
-13.938	3.87E-02
-14.089	5.18E-02
-14.244	6.28E-02
-14.403	6.42E-02
-14.564	6.06E-02
-14.730	5.30E-02
-14.899	4.68E-02
-15.072	4.00E-02
-15.250	2.85E-02
-15.431	1.97E-02
-15.617	1.94E-02
-15.808	2.16E-02
-16.003	2.07E-02
-16.203	1.86E-02
-16.408	2.12E-02
-16.618	2.41E-02
-16.834	2.32E-02
-17.056	2.41E-02
-17.283	3.07E-02
-17.517	3.56E-02
-17.757	3.87E-02
-18.003	3.99E-02
-18.257	4.02E-02
-18.518	4.26E-02
-18.786	4.38E-02
-19.062	4.41E-02
-19.347	4.57E-02
-19.640	4.47E-02
-19.942	4.20E-02
\ No newline at end of file
+9.9710	1.97E-01
+10.048	1.09E-01
+10.127	6.59E-02
+10.207	4.40E-02
+10.288	3.20E-02
+10.370	2.55E-02
+10.453	2.21E-02
+10.538	2.00E-02
+10.625	1.82E-02
+10.713	1.78E-02
+10.802	1.99E-02
+10.893	2.47E-02
+10.985	2.98E-02
+11.079	3.55E-02
+11.174	4.52E-02
+11.272	6.50E-02
+11.370	1.13E-01
+11.471	1.75E-01
+11.573	2.04E-01
+11.678	2.12E-01
+11.784	2.32E-01
+11.892	2.31E-01
+12.002	2.43E-01
+12.114	3.45E-01
+12.229	4.91E-01
+12.345	6.36E-01
+12.464	7.11E-01
+12.585	6.02E-01
+12.708	3.58E-01
+12.834	1.70E-01
+12.962	9.71E-02
+13.093	8.80E-02
+13.227	8.51E-02
+13.363	6.13E-02
+13.502	4.10E-02
+13.645	3.61E-02
+13.790	3.56E-02
+13.938	3.95E-02
+14.089	5.24E-02
+14.244	6.53E-02
+14.403	6.98E-02
+14.564	6.55E-02
+14.730	5.43E-02
+14.899	4.63E-02
+15.072	4.27E-02
+15.250	3.37E-02
+15.431	2.33E-02
+15.617	2.18E-02
+15.808	2.24E-02
+16.003	1.97E-02
+16.203	1.92E-02
+16.408	2.21E-02
+16.618	2.57E-02
+16.834	2.96E-02
+17.056	3.42E-02
+17.283	3.82E-02
+17.517	4.15E-02
+17.757	4.37E-02
+18.003	4.80E-02
+18.257	5.46E-02
+18.518	5.65E-02
+18.786	5.49E-02
+19.062	5.22E-02
+19.347	4.96E-02
+19.640	4.96E-02
+19.942	5.50E-02
\ No newline at end of file
diff --git a/data2.txt b/data2.txt
new file mode 100644
index 0000000..64e89e7
--- /dev/null
+++ b/data2.txt
@@ -0,0 +1,1533 @@
+wl	n
+2.5000	1.386
+2.5031	1.386
+2.5063	1.386
+2.5094	1.386
+2.5126	1.386
+2.5157	1.386
+2.5189	1.386
+2.5221	1.386
+2.5253	1.386
+2.5284	1.386
+2.5316	1.386
+2.5349	1.385
+2.5381	1.385
+2.5413	1.385
+2.5445	1.386
+2.5478	1.386
+2.5510	1.386
+2.5543	1.386
+2.5575	1.385
+2.5608	1.385
+2.5641	1.385
+2.5674	1.386
+2.5707	1.385
+2.5740	1.385
+2.5773	1.385
+2.5806	1.385
+2.5840	1.385
+2.5873	1.385
+2.5907	1.385
+2.5940	1.385
+2.5974	1.385
+2.6008	1.385
+2.6042	1.385
+2.6076	1.384
+2.6110	1.384
+2.6144	1.384
+2.6178	1.384
+2.6212	1.385
+2.6247	1.384
+2.6281	1.384
+2.6316	1.384
+2.6350	1.384
+2.6385	1.384
+2.6420	1.384
+2.6455	1.384
+2.6490	1.384
+2.6525	1.384
+2.6560	1.384
+2.6596	1.384
+2.6631	1.383
+2.6667	1.383
+2.6702	1.384
+2.6738	1.383
+2.6774	1.384
+2.6810	1.384
+2.6846	1.384
+2.6882	1.384
+2.6918	1.384
+2.6954	1.384
+2.6991	1.384
+2.7027	1.384
+2.7064	1.384
+2.7100	1.384
+2.7137	1.384
+2.7174	1.384
+2.7211	1.385
+2.7248	1.384
+2.7285	1.384
+2.7322	1.384
+2.7360	1.384
+2.7397	1.384
+2.7435	1.384
+2.7473	1.385
+2.7510	1.384
+2.7548	1.384
+2.7586	1.384
+2.7624	1.384
+2.7663	1.385
+2.7701	1.385
+2.7739	1.385
+2.7778	1.384
+2.7816	1.384
+2.7855	1.384
+2.7894	1.384
+2.7933	1.384
+2.7972	1.384
+2.8011	1.383
+2.8050	1.384
+2.8090	1.384
+2.8129	1.384
+2.8169	1.383
+2.8209	1.384
+2.8249	1.384
+2.8289	1.383
+2.8329	1.383
+2.8369	1.383
+2.8409	1.382
+2.8450	1.383
+2.8490	1.382
+2.8531	1.383
+2.8571	1.383
+2.8612	1.383
+2.8653	1.383
+2.8694	1.383
+2.8736	1.383
+2.8777	1.382
+2.8818	1.382
+2.8860	1.383
+2.8902	1.383
+2.8944	1.382
+2.8986	1.382
+2.9028	1.382
+2.9070	1.382
+2.9112	1.382
+2.9155	1.382
+2.9197	1.382
+2.9240	1.382
+2.9283	1.382
+2.9326	1.382
+2.9369	1.382
+2.9412	1.382
+2.9455	1.382
+2.9499	1.382
+2.9542	1.382
+2.9586	1.382
+2.9630	1.381
+2.9674	1.381
+2.9718	1.381
+2.9762	1.381
+2.9806	1.381
+2.9851	1.381
+2.9895	1.381
+2.9940	1.381
+2.9985	1.381
+3.0030	1.381
+3.0075	1.381
+3.0120	1.381
+3.0166	1.380
+3.0211	1.380
+3.0257	1.380
+3.0303	1.380
+3.0349	1.380
+3.0395	1.380
+3.0441	1.380
+3.0488	1.380
+3.0534	1.380
+3.0581	1.379
+3.0628	1.380
+3.0675	1.379
+3.0722	1.379
+3.0769	1.379
+3.0817	1.379
+3.0864	1.379
+3.0912	1.379
+3.0960	1.379
+3.1008	1.378
+3.1056	1.378
+3.1104	1.378
+3.1153	1.378
+3.1201	1.378
+3.1250	1.378
+3.1299	1.378
+3.1348	1.378
+3.1397	1.378
+3.1447	1.378
+3.1496	1.378
+3.1546	1.378
+3.1596	1.378
+3.1646	1.377
+3.1696	1.377
+3.1746	1.377
+3.1797	1.377
+3.1847	1.376
+3.1898	1.377
+3.1949	1.376
+3.2000	1.376
+3.2051	1.376
+3.2103	1.377
+3.2154	1.376
+3.2206	1.376
+3.2258	1.376
+3.2310	1.376
+3.2362	1.376
+3.2415	1.375
+3.2468	1.375
+3.2520	1.375
+3.2573	1.375
+3.2626	1.374
+3.2680	1.374
+3.2733	1.374
+3.2787	1.373
+3.2841	1.373
+3.2895	1.373
+3.2949	1.372
+3.3003	1.371
+3.3058	1.371
+3.3113	1.370
+3.3167	1.369
+3.3223	1.368
+3.3278	1.367
+3.3333	1.365
+3.3389	1.363
+3.3445	1.360
+3.3501	1.354
+3.3557	1.350
+3.3613	1.348
+3.3670	1.354
+3.3727	1.372
+3.3784	1.394
+3.3841	1.406
+3.3898	1.407
+3.3956	1.402
+3.4014	1.398
+3.4072	1.394
+3.4130	1.391
+3.4188	1.389
+3.4247	1.387
+3.4305	1.385
+3.4364	1.386
+3.4423	1.386
+3.4483	1.388
+3.4542	1.388
+3.4602	1.388
+3.4662	1.388
+3.4722	1.386
+3.4783	1.386
+3.4843	1.385
+3.4904	1.385
+3.4965	1.384
+3.5026	1.384
+3.5088	1.383
+3.5149	1.383
+3.5211	1.383
+3.5273	1.383
+3.5336	1.383
+3.5398	1.382
+3.5461	1.382
+3.5524	1.381
+3.5587	1.382
+3.5651	1.382
+3.5714	1.381
+3.5778	1.382
+3.5842	1.382
+3.5907	1.382
+3.5971	1.381
+3.6036	1.381
+3.6101	1.381
+3.6166	1.380
+3.6232	1.380
+3.6298	1.380
+3.6364	1.380
+3.6430	1.380
+3.6496	1.379
+3.6563	1.379
+3.6630	1.379
+3.6697	1.379
+3.6765	1.379
+3.6832	1.379
+3.6900	1.379
+3.6969	1.378
+3.7037	1.378
+3.7106	1.378
+3.7175	1.378
+3.7244	1.378
+3.7313	1.378
+3.7383	1.378
+3.7453	1.378
+3.7523	1.378
+3.7594	1.378
+3.7665	1.377
+3.7736	1.377
+3.7807	1.377
+3.7879	1.377
+3.7951	1.377
+3.8023	1.376
+3.8095	1.376
+3.8168	1.376
+3.8241	1.376
+3.8314	1.376
+3.8388	1.376
+3.8462	1.376
+3.8536	1.376
+3.8610	1.375
+3.8685	1.375
+3.8760	1.375
+3.8835	1.375
+3.8911	1.375
+3.8986	1.375
+3.9063	1.375
+3.9139	1.375
+3.9216	1.375
+3.9293	1.375
+3.9370	1.375
+3.9448	1.375
+3.9526	1.375
+3.9604	1.375
+3.9683	1.374
+3.9761	1.374
+3.9841	1.374
+3.9920	1.374
+4.0000	1.374
+4.0080	1.374
+4.0161	1.374
+4.0241	1.374
+4.0323	1.373
+4.0404	1.373
+4.0486	1.373
+4.0568	1.373
+4.0650	1.373
+4.0733	1.373
+4.0816	1.372
+4.0900	1.372
+4.0984	1.372
+4.1068	1.372
+4.1152	1.372
+4.1237	1.372
+4.1322	1.372
+4.1408	1.372
+4.1494	1.372
+4.1580	1.371
+4.1667	1.371
+4.1754	1.371
+4.1841	1.370
+4.1929	1.370
+4.2017	1.370
+4.2105	1.370
+4.2194	1.370
+4.2283	1.370
+4.2373	1.370
+4.2463	1.369
+4.2553	1.369
+4.2644	1.369
+4.2735	1.368
+4.2827	1.368
+4.2918	1.368
+4.3011	1.368
+4.3103	1.368
+4.3197	1.368
+4.3290	1.368
+4.3384	1.368
+4.3478	1.367
+4.3573	1.367
+4.3668	1.367
+4.3856	1.367
+4.3764	1.367
+4.3956	1.367
+4.4053	1.367
+4.4150	1.367
+4.4248	1.366
+4.4446	1.366
+4.4444	1.367
+4.4543	1.366
+4.4643	1.366
+4.4743	1.366
+4.4843	1.366
+4.4944	1.366
+4.5045	1.365
+4.5147	1.365
+4.5249	1.365
+4.5351	1.365
+4.5455	1.365
+4.5558	1.365
+4.5662	1.365
+4.5767	1.364
+4.5872	1.364
+4.5977	1.364
+4.6083	1.364
+4.6189	1.364
+4.6296	1.364
+4.6404	1.364
+4.6512	1.364
+4.6620	1.363
+4.6729	1.363
+4.6838	1.363
+4.6948	1.363
+4.7059	1.363
+4.7170	1.363
+4.7281	1.362
+4.7393	1.362
+4.7506	1.362
+4.7619	1.362
+4.7733	1.362
+4.7847	1.361
+4.7962	1.361
+4.8077	1.361
+4.8193	1.360
+4.8309	1.360
+4.8426	1.360
+4.8544	1.360
+4.8662	1.360
+4.8780	1.360
+4.8900	1.360
+4.9020	1.360
+4.9140	1.359
+4.9261	1.359
+4.9383	1.359
+4.9505	1.358
+4.9628	1.358
+4.9751	1.358
+4.9875	1.358
+5.0000	1.358
+5.0125	1.358
+5.0251	1.358
+5.0378	1.358
+5.0505	1.358
+5.0633	1.358
+5.0761	1.358
+5.0891	1.357
+5.1020	1.357
+5.1151	1.357
+5.1282	1.357
+5.1414	1.357
+5.1546	1.358
+5.1680	1.357
+5.1813	1.357
+5.1948	1.357
+5.2083	1.357
+5.2219	1.356
+5.2356	1.356
+5.2493	1.356
+5.2632	1.355
+5.2770	1.355
+5.2910	1.354
+5.3050	1.354
+5.3191	1.354
+5.3333	1.354
+5.3476	1.354
+5.3619	1.354
+5.3763	1.354
+5.4054	1.352
+5.4201	1.352
+5.4201	1.352
+5.4348	1.351
+5.4496	1.351
+5.4645	1.350
+5.4795	1.350
+5.4945	1.349
+5.5096	1.349
+5.5249	1.349
+5.5402	1.349
+5.5556	1.349
+5.5710	1.348
+5.5866	1.348
+5.6022	1.346
+5.6180	1.348
+5.6338	1.347
+5.6497	1.347
+5.6657	1.347
+5.6818	1.346
+5.6980	1.346
+5.7143	1.345
+5.7307	1.345
+5.7471	1.344
+5.7637	1.344
+5.7803	1.344
+5.7971	1.343
+5.8140	1.343
+5.8309	1.343
+5.8480	1.342
+5.8651	1.342
+5.882	1.342
+5.8997	1.341
+5.9172	1.341
+5.9347	1.340
+5.9524	1.340
+5.9701	1.339
+5.9880	1.339
+6.0060	1.338
+6.0241	1.338
+6.0423	1.337
+6.0606	1.337
+6.0790	1.336
+6.0976	1.336
+6.1162	1.336
+6.1350	1.335
+6.1538	1.334
+6.1728	1.334
+6.1920	1.333
+6.2112	1.333
+6.2305	1.332
+6.2500	1.332
+6.2696	1.331
+6.2893	1.331
+6.3091	1.330
+6.3291	1.329
+6.3492	1.329
+6.3694	1.328
+6.3898	1.328
+6.4103	1.328
+6.4309	1.328
+6.4516	1.326
+6.4725	1.324
+6.4935	1.324
+6.5147	1.324
+6.5359	1.322
+6.5574	1.322
+6.5789	1.322
+6.6007	1.320
+6.6225	1.318
+6.6445	1.318
+6.6667	1.316
+6.6890	1.315
+6.7114	1.315
+6.7340	1.313
+6.7568	1.311
+6.7797	1.310
+6.8027	1.309
+6.8259	1.306
+6.8493	1.304
+6.8729	1.303
+6.8966	1.303
+6.9204	1.303
+6.9444	1.302
+6.9686	1.300
+6.9930	1.297
+7.0175	1.293
+7.0423	1.292
+7.0671	1.295
+7.0922	1.299
+7.1174	1.300
+7.1429	1.302
+7.1685	1.304
+7.1942	1.302
+7.2202	1.297
+7.2464	1.293
+7.2727	1.290
+7.2993	1.287
+7.3260	1.284
+7.3529	1.281
+7.3801	1.278
+7.4074	1.275
+7.4349	1.271
+7.4627	1.267
+7.4906	1.262
+7.5188	1.258
+7.5472	1.252
+7.5758	1.247
+7.6046	1.240
+7.6336	1.232
+7.6628	1.223
+7.6923	1.212
+7.7220	1.200
+7.7519	1.183
+7.7821	1.161
+7.8225	1.125
+7.8431	1.063
+7.8740	1.009
+7.9051	1.127
+7.9365	1.443
+7.9681	1.551
+8.0000	1.474
+8.0321	1.409
+8.0645	1.367
+8.0972	1.339
+8.1301	1.319
+8.1633	1.303
+8.1967	1.289
+8.2305	1.276
+8.2645	1.263
+8.2988	1.251
+8.3333	1.239
+8.3682	1.228
+8.4034	1.217
+8.4388	1.203
+8.4746	1.191
+8.5106	1.176
+8.5470	1.163
+8.5837	1.150
+8.6207	1.133
+8.6580	1.115
+8.6957	1.095
+8.7336	1.072
+8.7719	1.048
+8.8106	1.021
+8.8496	0.991
+8.8889	0.961
+8.9286	0.938
+8.9686	0.932
+9.0090	0.956
+9.0498	1.011
+9.0909	1.092
+9.1324	1.182
+9.1743	1.261
+9.2166	1.315
+9.2593	1.346
+9.3023	1.361
+9.3458	1.366
+9.3897	1.367
+9.4340	1.366
+9.4787	1.367
+9.5238	1.372
+9.5694	1.385
+9.6154	1.408
+9.6618	1.445
+9.7087	1.504
+9.7561	1.592
+9.8039	1.708
+9.8522	1.822
+9.9010	1.877
+9.9502	1.860
+10.0000	1.807
+10.0503	1.756
+10.1010	1.710
+10.1523	1.672
+10.2041	1.639
+10.2564	1.612
+10.3093	1.587
+10.3627	1.565
+10.4167	1.546
+10.4712	1.529
+10.5263	1.512
+10.5820	1.497
+10.6383	1.481
+10.6952	1.466
+10.7527	1.452
+10.8108	1.439
+10.8696	1.425
+10.9290	1.410
+10.9890	1.396
+11.0497	1.381
+11.1111	1.366
+11.1732	1.350
+11.2360	1.333
+11.2994	1.318
+11.3636	1.309
+11.4286	1.307
+11.4943	1.307
+11.5607	1.309
+11.6279	1.305
+11.6959	1.287
+11.7647	1.258
+11.8343	1.230
+11.9048	1.196
+11.9760	1.142
+12.0482	1.075
+12.1212	1.041
+12.1951	1.078
+12.2699	1.143
+12.3457	1.241
+12.4224	1.442
+12.5000	1.744
+12.5786	2.005
+12.6582	2.095
+12.7389	2.043
+12.8205	1.951
+12.9032	1.866
+12.9870	1.799
+13.0719	1.748
+13.1579	1.710
+13.2450	1.685
+13.3333	1.666
+13.4228	1.646
+13.5135	1.629
+13.6054	1.612
+13.6986	1.596
+13.7931	1.579
+13.8889	1.564
+13.9860	1.549
+14.0845	1.540
+14.1844	1.546
+14.2857	1.553
+14.3885	1.556
+14.4928	1.558
+14.5985	1.562
+14.7059	1.565
+14.8148	1.562
+14.9254	1.558
+15.0376	1.562
+15.1515	1.572
+15.2672	1.568
+15.3846	1.554
+15.5039	1.545
+15.6250	1.538
+15.7480	1.535
+15.8730	1.535
+16.0000	1.530
+16.1290	1.527
+16.2602	1.523
+16.3934	1.520
+16.5289	1.517
+16.6667	1.514
+16.8067	1.510
+16.9492	1.508
+17.0940	1.504
+17.2414	1.501
+17.3913	1.498
+17.5439	1.495
+17.6991	1.492
+17.8571	1.489
+18.0180	1.486
+18.1818	1.485
+18.3486	1.482
+18.5185	1.480
+18.6916	1.478
+18.8679	1.476
+19.0476	1.473
+19.2308	1.471
+19.4175	1.470
+19.6078	1.468
+19.8020	1.466
+20.0000	1.465
+20.2020	1.464
+20.4082	1.464
+20.6186	1.462
+20.8333	1.460
+21.0526	1.455
+21.2766	1.450
+21.5054	1.444
+21.7391	1.437
+21.9780	1.431
+22.2222	1.423
+22.4719	1.411
+22.7273	1.400
+22.9885	1.387
+23.2558	1.370
+23.5294	1.348
+23.8095	1.323
+24.0964	1.307
+24.3902	1.311
+24.6914	1.336
+25.0000	1.379
+25.3165	1.433
+25.6410	1.485
+25.9740	1.530
+26.3158	1.567
+26.6667	1.589
+27.0270	1.597
+27.3973	1.596
+27.7778	1.590
+28.1690	1.582
+28.5714	1.576
+28.9855	1.566
+29.4118	1.563
+29.8507	1.560
+30.3030	1.555
+30.7692	1.550
+31.2500	1.546
+31.7460	1.543
+32.2581	1.540
+32.7869	1.537
+33.3333	1.536
+33.8983	1.536
+34.4828	1.536
+35.0877	1.539
+35.7143	1.544
+36.3636	1.549
+37.0370	1.551
+37.7358	1.550
+38.4615	1.554
+39.2157	1.555
+40.0000	1.555
+40.8163	1.553
+41.6667	1.551
+42.5532	1.546
+43.4783	1.542
+44.4444	1.539
+45.4545	1.535
+46.5116	1.532
+47.6190	1.532
+48.7805	1.531
+50.0000	1.529
+51.2821	1.530
+52.6316	1.530
+54.0541	1.531
+55.5556	1.533
+
+wl	k
+2.5000	0.000
+2.5031	0.000
+2.5063	0.000
+2.5094	0.000
+2.5126	0.000
+2.5157	0.000
+2.5189	0.000
+2.5221	0.000
+2.5253	0.000
+2.5284	0.000
+2.5316	0.000
+2.5349	0.000
+2.5381	0.000
+2.5413	0.000
+2.5445	0.000
+2.5478	0.000
+2.5510	0.000
+2.5543	0.000
+2.5575	0.000
+2.5608	0.000
+2.5641	0.000
+2.5674	0.000
+2.5707	0.000
+2.5740	0.000
+2.5773	0.000
+2.5806	0.000
+2.5840	0.000
+2.5873	0.000
+2.5907	0.000
+2.5940	0.000
+2.5974	0.000
+2.6008	0.000
+2.6042	0.000
+2.6076	0.000
+2.6110	0.000
+2.6144	0.000
+2.6178	0.000
+2.6212	0.000
+2.6247	0.000
+2.6281	0.000
+2.6316	0.000
+2.6350	0.000
+2.6385	0.000
+2.6420	0.000
+2.6455	0.000
+2.6490	0.000
+2.6525	0.000
+2.6560	0.000
+2.6596	0.000
+2.6631	0.000
+2.6667	0.001
+2.6702	0.000
+2.6738	0.001
+2.6774	0.001
+2.6810	0.001
+2.6846	0.000
+2.6882	0.001
+2.6918	0.001
+2.6954	0.001
+2.6991	0.001
+2.7027	0.001
+2.7064	0.001
+2.7100	0.001
+2.7137	0.001
+2.7174	0.001
+2.7211	0.001
+2.7248	0.000
+2.7285	0.001
+2.7322	0.001
+2.7360	0.001
+2.7397	0.001
+2.7435	0.001
+2.7473	0.000
+2.7510	0.001
+2.7548	0.001
+2.7586	0.001
+2.7624	0.001
+2.7663	0.001
+2.7701	0.000
+2.7739	0.000
+2.7778	0.000
+2.7816	0.000
+2.7855	0.000
+2.7894	0.000
+2.7933	0.000
+2.7972	0.000
+2.8011	0.000
+2.8050	0.000
+2.8090	0.000
+2.8129	0.000
+2.8169	0.000
+2.8209	0.000
+2.8249	0.000
+2.8289	0.000
+2.8329	0.000
+2.8369	0.000
+2.8409	0.000
+2.8450	0.000
+2.8490	0.000
+2.8531	0.000
+2.8571	0.000
+2.8612	0.000
+2.8653	0.000
+2.8694	0.000
+2.8736	0.000
+2.8777	0.000
+2.8818	0.000
+2.8860	0.000
+2.8902	0.000
+2.8944	0.000
+2.8986	0.000
+2.9028	0.000
+2.9070	0.000
+2.9112	0.000
+2.9155	0.000
+2.9197	0.000
+2.9240	0.000
+2.9283	0.000
+2.9326	0.000
+2.9369	0.000
+2.9412	0.000
+2.9455	0.000
+2.9499	0.000
+2.9542	0.000
+2.9586	0.000
+2.9630	0.000
+2.9674	0.000
+2.9718	0.000
+2.9762	0.000
+2.9806	0.000
+2.9851	0.000
+2.9895	0.000
+2.9940	0.000
+2.9985	0.000
+3.0030	0.000
+3.0075	0.000
+3.0120	0.000
+3.0166	0.000
+3.0211	0.000
+3.0257	0.000
+3.0303	0.000
+3.0349	0.000
+3.0395	0.000
+3.0441	0.000
+3.0488	0.000
+3.0534	0.000
+3.0581	0.000
+3.0628	0.000
+3.0675	0.000
+3.0722	0.000
+3.0769	0.000
+3.0817	0.000
+3.0864	0.000
+3.0912	0.000
+3.0960	0.000
+3.1008	0.000
+3.1056	0.000
+3.1104	0.000
+3.1153	0.000
+3.1201	0.000
+3.1250	0.000
+3.1299	0.000
+3.1348	0.000
+3.1397	0.000
+3.1447	0.000
+3.1496	0.000
+3.1546	0.000
+3.1596	0.000
+3.1646	0.000
+3.1696	0.000
+3.1746	0.000
+3.1797	0.000
+3.1847	0.000
+3.1898	0.001
+3.1949	0.000
+3.2000	0.001
+3.2051	0.001
+3.2103	0.001
+3.2154	0.001
+3.2206	0.001
+3.2258	0.001
+3.2310	0.001
+3.2362	0.001
+3.2415	0.000
+3.2468	0.001
+3.2520	0.000
+3.2573	0.000
+3.2626	0.001
+3.2680	0.001
+3.2733	0.001
+3.2787	0.001
+3.2841	0.001
+3.2895	0.001
+3.2949	0.001
+3.3003	0.001
+3.3058	0.001
+3.3113	0.001
+3.3167	0.002
+3.3223	0.000
+3.3278	0.002
+3.3333	0.003
+3.3389	0.004
+3.3445	0.005
+3.3501	0.008
+3.3557	0.015
+3.3613	0.028
+3.3670	0.043
+3.3727	0.053
+3.3784	0.050
+3.3841	0.035
+3.3898	0.019
+3.3956	0.011
+3.4014	0.007
+3.4072	0.005
+3.4130	0.004
+3.4188	0.004
+3.4247	0.004
+3.4305	0.005
+3.4364	0.006
+3.4423	0.007
+3.4483	0.006
+3.4542	0.004
+3.4602	0.004
+3.4662	0.003
+3.4722	0.002
+3.4783	0.002
+3.4843	0.001
+3.4904	0.000
+3.4965	0.000
+3.5026	0.001
+3.5088	0.001
+3.5149	0.001
+3.5211	0.001
+3.5273	0.001
+3.5336	0.001
+3.5398	0.000
+3.5461	0.001
+3.5524	0.001
+3.5587	0.001
+3.5651	0.000
+3.5714	0.001
+3.5778	0.001
+3.5842	0.001
+3.5907	0.000
+3.5971	0.000
+3.6036	0.000
+3.6101	0.000
+3.6166	0.000
+3.6232	0.000
+3.6298	0.000
+3.6364	0.000
+3.6430	0.000
+3.6496	0.000
+3.6563	0.000
+3.6630	0.000
+3.6697	0.000
+3.6765	0.000
+3.6832	0.000
+3.6900	0.000
+3.6969	0.000
+3.7037	0.000
+3.7106	0.000
+3.7175	0.000
+3.7244	0.000
+3.7313	0.000
+3.7383	0.000
+3.7453	0.000
+3.7523	0.000
+3.7594	0.000
+3.7665	0.000
+3.7736	0.000
+3.7807	0.000
+3.7879	0.000
+3.7951	0.000
+3.8023	0.000
+3.8095	0.000
+3.8168	0.000
+3.8241	0.000
+3.8314	0.000
+3.8388	0.000
+3.8462	0.000
+3.8536	0.000
+3.8610	0.000
+3.8685	0.000
+3.8760	0.000
+3.8835	0.000
+3.8911	0.000
+3.8986	0.000
+3.9063	0.000
+3.9139	0.000
+3.9216	0.000
+3.9293	0.000
+3.9370	0.000
+3.9448	0.000
+3.9526	0.000
+3.9604	0.000
+3.9683	0.000
+3.9761	0.000
+3.9841	0.000
+3.9920	0.000
+4.0000	0.000
+4.0080	0.000
+4.0161	0.000
+4.0241	0.000
+4.0323	0.000
+4.0404	0.000
+4.0486	0.000
+4.0568	0.000
+4.0650	0.000
+4.0733	0.000
+4.0816	0.000
+4.0900	0.000
+4.0984	0.000
+4.1068	0.000
+4.1152	0.000
+4.1237	0.000
+4.1322	0.000
+4.1408	0.000
+4.1494	0.000
+4.1580	0.000
+4.1667	0.000
+4.1754	0.000
+4.1841	0.000
+4.1929	0.000
+4.2017	0.000
+4.2105	0.000
+4.2194	0.000
+4.2283	0.000
+4.2373	0.000
+4.2463	0.000
+4.2553	0.000
+4.2644	0.000
+4.2735	0.000
+4.2827	0.000
+4.2918	0.000
+4.3011	0.000
+4.3103	0.000
+4.3197	0.000
+4.3290	0.000
+4.3384	0.000
+4.3478	0.000
+4.3573	0.000
+4.3668	0.000
+4.3856	0.000
+4.3764	0.000
+4.3956	0.000
+4.4053	0.000
+4.4150	0.000
+4.4248	0.000
+4.4446	0.001
+4.4444	0.000
+4.4543	0.001
+4.4643	0.001
+4.4743	0.001
+4.4843	0.001
+4.4944	0.001
+4.5045	0.001
+4.5147	0.001
+4.5249	0.001
+4.5351	0.001
+4.5455	0.001
+4.5558	0.001
+4.5662	0.001
+4.5767	0.001
+4.5872	0.001
+4.5977	0.001
+4.6083	0.001
+4.6189	0.001
+4.6296	0.001
+4.6404	0.001
+4.6512	0.001
+4.6620	0.001
+4.6729	0.001
+4.6838	0.001
+4.6948	0.001
+4.7059	0.001
+4.7170	0.001
+4.7281	0.002
+4.7393	0.002
+4.7506	0.002
+4.7619	0.001
+4.7733	0.001
+4.7847	0.002
+4.7962	0.002
+4.8077	0.002
+4.8193	0.002
+4.8309	0.002
+4.8426	0.002
+4.8544	0.002
+4.8662	0.002
+4.8780	0.002
+4.8900	0.002
+4.9020	0.002
+4.9140	0.002
+4.9261	0.002
+4.9383	0.002
+4.9505	0.002
+4.9628	0.003
+4.9751	0.003
+4.9875	0.004
+5.0000	0.004
+5.0125	0.004
+5.0251	0.004
+5.0378	0.004
+5.0505	0.004
+5.0633	0.004
+5.0761	0.004
+5.0891	0.004
+5.1020	0.004
+5.1151	0.004
+5.1282	0.004
+5.1414	0.004
+5.1546	0.004
+5.1680	0.004
+5.1813	0.004
+5.1948	0.003
+5.2083	0.003
+5.2219	0.003
+5.2356	0.003
+5.2493	0.003
+5.2632	0.003
+5.2770	0.003
+5.2910	0.003
+5.3050	0.003
+5.3191	0.003
+5.3333	0.003
+5.3476	0.003
+5.3619	0.003
+5.3763	0.002
+5.4054	0.002
+5.4201	0.002
+5.4201	0.002
+5.4348	0.002
+5.4496	0.003
+5.4645	0.003
+5.4795	0.003
+5.4945	0.003
+5.5096	0.003
+5.5249	0.004
+5.5402	0.003
+5.5556	0.004
+5.5710	0.004
+5.5866	0.004
+5.6022	0.004
+5.6180	0.004
+5.6338	0.004
+5.6497	0.004
+5.6657	0.004
+5.6818	0.004
+5.6980	0.004
+5.7143	0.004
+5.7307	0.004
+5.7471	0.005
+5.7637	0.005
+5.7803	0.005
+5.7971	0.005
+5.8140	0.005
+5.8309	0.005
+5.8480	0.005
+5.8651	0.005
+5.882	0.005
+5.8997	0.005
+5.9172	0.006
+5.9347	0.005
+5.9524	0.005
+5.9701	0.005
+5.9880	0.006
+6.0060	0.006
+6.0241	0.006
+6.0423	0.006
+6.0606	0.006
+6.0790	0.006
+6.0976	0.006
+6.1162	0.006
+6.1350	0.006
+6.1538	0.006
+6.1728	0.007
+6.1920	0.007
+6.2112	0.007
+6.2305	0.007
+6.2500	0.007
+6.2696	0.007
+6.2893	0.007
+6.3091	0.007
+6.3291	0.008
+6.3492	0.008
+6.3694	0.007
+6.3898	0.008
+6.4103	0.008
+6.4309	0.007
+6.4516	0.006
+6.4725	0.007
+6.4935	0.008
+6.5147	0.007
+6.5359	0.008
+6.5574	0.009
+6.5789	0.008
+6.6007	0.006
+6.6225	0.007
+6.6445	0.007
+6.6667	0.007
+6.6890	0.008
+6.7114	0.008
+6.7340	0.007
+6.7568	0.008
+6.7797	0.009
+6.8027	0.008
+6.8259	0.008
+6.8493	0.011
+6.8729	0.012
+6.8966	0.014
+6.9204	0.014
+6.9444	0.014
+6.9686	0.013
+6.9930	0.014
+7.0175	0.016
+7.0423	0.022
+7.0671	0.026
+7.0922	0.026
+7.1174	0.025
+7.1429	0.023
+7.1685	0.020
+7.1942	0.015
+7.2202	0.012
+7.2464	0.012
+7.2727	0.012
+7.2993	0.012
+7.3260	0.012
+7.3529	0.012
+7.3801	0.012
+7.4074	0.012
+7.4349	0.012
+7.4627	0.012
+7.4906	0.012
+7.5188	0.013
+7.5472	0.014
+7.5758	0.015
+7.6046	0.015
+7.6336	0.016
+7.6628	0.018
+7.6923	0.020
+7.7220	0.025
+7.7519	0.030
+7.7821	0.037
+7.8225	0.051
+7.8431	0.092
+7.8740	0.232
+7.9051	0.448
+7.9365	0.450
+7.9681	0.201
+8.0000	0.070
+8.0321	0.032
+8.0645	0.020
+8.0972	0.016
+8.1301	0.014
+8.1633	0.013
+8.1967	0.013
+8.2305	0.012
+8.2645	0.013
+8.2988	0.014
+8.3333	0.016
+8.3682	0.018
+8.4034	0.019
+8.4388	0.022
+8.4746	0.026
+8.5106	0.030
+8.5470	0.037
+8.5837	0.042
+8.6207	0.048
+8.6580	0.057
+8.6957	0.067
+8.7336	0.081
+8.7719	0.098
+8.8106	0.121
+8.8496	0.154
+8.8889	0.199
+8.9286	0.262
+8.9686	0.341
+9.0090	0.428
+9.0498	0.507
+9.0909	0.561
+9.1324	0.578
+9.1743	0.564
+9.2166	0.531
+9.2593	0.496
+9.3023	0.467
+9.3458	0.448
+9.3897	0.438
+9.4340	0.438
+9.4787	0.446
+9.5238	0.462
+9.5694	0.486
+9.6154	0.514
+9.6618	0.547
+9.7087	0.580
+9.7561	0.600
+9.8039	0.582
+9.8522	0.499
+9.9010	0.363
+9.9502	0.236
+10.0000	0.154
+10.0503	0.106
+10.1010	0.077
+10.1523	0.058
+10.2041	0.045
+10.2564	0.037
+10.3093	0.030
+10.3627	0.026
+10.4167	0.024
+10.4712	0.020
+10.5263	0.019
+10.5820	0.017
+10.6383	0.017
+10.6952	0.017
+10.7527	0.017
+10.8108	0.019
+10.8696	0.020
+10.9290	0.022
+10.9890	0.025
+11.0497	0.029
+11.1111	0.034
+11.1732	0.043
+11.2360	0.053
+11.2994	0.069
+11.3636	0.089
+11.4286	0.106
+11.4943	0.118
+11.5607	0.124
+11.6279	0.121
+11.6959	0.116
+11.7647	0.123
+11.8343	0.139
+11.9048	0.156
+11.9760	0.186
+12.0482	0.258
+12.1212	0.391
+12.1951	0.532
+12.2699	0.650
+12.3457	0.790
+12.4224	0.923
+12.5000	0.930
+12.5786	0.745
+12.6582	0.475
+12.7389	0.264
+12.8205	0.146
+12.9032	0.088
+12.9870	0.060
+13.0719	0.048
+13.1579	0.044
+13.2450	0.043
+13.3333	0.038
+13.4228	0.033
+13.5135	0.030
+13.6054	0.028
+13.6986	0.027
+13.7931	0.028
+13.8889	0.032
+13.9860	0.040
+14.0845	0.058
+14.1844	0.071
+14.2857	0.073
+14.3885	0.071
+14.4928	0.072
+14.5985	0.070
+14.7059	0.063
+14.8148	0.057
+14.9254	0.057
+15.0376	0.059
+15.1515	0.050
+15.2672	0.032
+15.3846	0.025
+15.5039	0.025
+15.6250	0.025
+15.7480	0.028
+15.8730	0.026
+16.0000	0.024
+16.1290	0.023
+16.2602	0.022
+16.3934	0.022
+16.5289	0.021
+16.6667	0.020
+16.8067	0.019
+16.9492	0.019
+17.0940	0.018
+17.2414	0.018
+17.3913	0.018
+17.5439	0.019
+17.6991	0.019
+17.8571	0.019
+18.0180	0.020
+18.1818	0.021
+18.3486	0.021
+18.5185	0.022
+18.6916	0.022
+18.8679	0.023
+19.0476	0.024
+19.2308	0.025
+19.4175	0.026
+19.6078	0.026
+19.8020	0.028
+20.0000	0.028
+20.2020	0.028
+20.4082	0.028
+20.6186	0.025
+20.8333	0.023
+21.0526	0.021
+21.2766	0.020
+21.5054	0.019
+21.7391	0.019
+21.9780	0.019
+22.2222	0.019
+22.4719	0.020
+22.7273	0.024
+22.9885	0.027
+23.2558	0.034
+23.5294	0.048
+23.8095	0.073
+24.0964	0.116
+24.3902	0.169
+24.6914	0.219
+25.0000	0.257
+25.3165	0.275
+25.6410	0.273
+25.9740	0.258
+26.3158	0.231
+26.6667	0.198
+27.0270	0.164
+27.3973	0.136
+27.7778	0.116
+28.1690	0.101
+28.5714	0.091
+28.9855	0.084
+29.4118	0.082
+29.8507	0.073
+30.3030	0.069
+30.7692	0.067
+31.2500	0.065
+31.7460	0.064
+32.2581	0.063
+32.7869	0.065
+33.3333	0.066
+33.8983	0.067
+34.4828	0.069
+35.0877	0.070
+35.7143	0.071
+36.3636	0.067
+37.0370	0.063
+37.7358	0.060
+38.4615	0.056
+39.2157	0.051
+40.0000	0.046
+40.8163	0.040
+41.6667	0.035
+42.5532	0.031
+43.4783	0.030
+44.4444	0.029
+45.4545	0.028
+46.5116	0.030
+47.6190	0.032
+48.7805	0.031
+50.0000	0.033
+51.2821	0.034
+52.6316	0.035
+54.0541	0.036
+55.5556	0.034
\ No newline at end of file
diff --git a/org/chatgpt2/q1_2.py b/org/chatgpt2/q1_2.py
index 933472e..c05c928 100644
--- a/org/chatgpt2/q1_2.py
+++ b/org/chatgpt2/q1_2.py
@@ -2,7 +2,21 @@ import numpy as np
 import matplotlib.pyplot as plt
 from scipy.interpolate import CubicSpline
 
+def make_strictly_increasing(wl, n, k):
+    # 去除完全相同的点
+    unique_wl, indices = np.unique(wl, return_index=True)
+    if len(unique_wl) != len(wl):
+        print(f"Removed {len(wl) - len(unique_wl)} duplicate wavelength points")
+        wl, n, k = wl[indices], n[indices], k[indices]
 
+    # 确保严格递增
+    is_increasing = np.diff(wl) > 0
+    if not all(is_increasing):
+        # 移除不满足严格递增的点
+        valid_indices = np.concatenate([[True], is_increasing])
+        wl, n, k = wl[valid_indices], n[valid_indices], k[valid_indices]
+
+    return wl, n, k
 # -----------------------------
 # 1. 从data.txt读取分块格式数据（先wl+n，再wl+k）
 # -----------------------------
@@ -62,6 +76,8 @@ def read_split_data(file_path):
 # 读取数据
 wl_all, n_all, k_all = read_split_data('/Users/spasolreisa/IdeaProjects/asiaMath/data.txt')
 
+wl_all, n_all, k_all = make_strictly_increasing(wl_all, n_all, k_all)
+#
 # 三次样条插值（覆盖全波段，保证计算精度）
 cs_n = CubicSpline(wl_all, n_all)  # 折射率n的插值函数
 cs_k = CubicSpline(wl_all, k_all)  # 消光系数k的插值函数
@@ -163,7 +179,7 @@ plt.xlim(wl_min, wl_max)
 
 # 保存图片（高分辨率）
 plt.tight_layout()
-plt.savefig('PDMS_emissivity_spectrum.png', dpi=300, bbox_inches='tight')
+# plt.savefig('PDMS_emissivity_spectrum.png', dpi=300, bbox_inches='tight')
 plt.show()
 
 # -----------------------------
diff --git a/org/chatgpt2/q1_3.py b/org/chatgpt2/q1_3.py
new file mode 100644
index 0000000..57182e9
--- /dev/null
+++ b/org/chatgpt2/q1_3.py
@@ -0,0 +1,347 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy.interpolate import CubicSpline
+from scipy.integrate import simpson
+
+
+# -----------------------------
+# 1. 材料级：PDMS光学性能计算（你的模型核心逻辑，修复变量定义错误）
+# -----------------------------
+def make_strictly_increasing(wl, n, k):
+    # 去除重复点并确保波长严格递增
+    unique_wl, indices = np.unique(wl, return_index=True)
+    if len(unique_wl) != len(wl):
+        print(f"Removed {len(wl) - len(unique_wl)} duplicate wavelength points")
+        wl, n, k = wl[indices], n[indices], k[indices]
+    # 确保严格递增
+    is_increasing = np.diff(wl) > 0
+    if not all(is_increasing):
+        valid_indices = np.concatenate([[True], is_increasing])
+        wl, n, k = wl[valid_indices], n[valid_indices], k[valid_indices]
+    return wl, n, k
+
+
+def read_split_data(file_path):
+    with open(file_path, 'r', encoding='utf-8') as f:
+        lines = [line.strip() for line in f if line.strip() and not line.startswith('#')]
+    split_idx = None
+    for i, line in enumerate(lines):
+        if line == 'wl	k':
+            split_idx = i
+            break
+    if split_idx is None:
+        raise ValueError("未找到'wl	k'表头，请检查数据格式！")
+    n_lines = lines[1:split_idx]
+    wl_n, n_list = [], []
+    for line in n_lines:
+        parts = line.split()
+        if len(parts) != 2:
+            continue  # 跳过格式错误的行
+        wl, n_val = parts
+        wl_n.append(float(wl)), n_list.append(float(n_val))
+    k_lines = lines[split_idx + 1:]
+    wl_k, k_list = [], []
+    for line in k_lines:
+        parts = line.split()
+        if len(parts) != 2:
+            continue  # 跳过格式错误的行
+        wl, k_val = parts
+        wl_k.append(float(wl)), k_list.append(float(k_val))
+    # 转换为numpy数组
+    wl_n, n_list = np.array(wl_n), np.array(n_list)
+    wl_k, k_list = np.array(wl_k), np.array(k_list)
+    # 确保n和k的波长完全一致
+    assert np.allclose(wl_n, wl_k), "n和k的波长列表不一致！"
+    # 排序
+    sorted_idx = np.argsort(wl_n)
+    return wl_n[sorted_idx], n_list[sorted_idx], k_list[sorted_idx]
+
+
+def fresnel_reflectance(n1, k1, n2, k2):
+    m1, m2 = n1 + 1j * k1, n2 + 1j * k2
+    return np.abs((m1 - m2) / (m1 + m2)) ** 2
+
+
+def thin_film_optical_properties(n_film, k_film, d, wl):
+    """修复denominator未定义的错误，完整计算R_total和T_total"""
+    R12 = fresnel_reflectance(1.0, 0.0, n_film, k_film)  # 空气→薄膜
+    R23 = fresnel_reflectance(n_film, k_film, 1.0, 0.0)  # 薄膜→空气
+    delta = 2 * np.pi * n_film * d / wl  # 干涉相位差
+    alpha_abs = 4 * np.pi * k_film * d / wl  # 吸收衰减系数
+
+    # 计算分母（关键修复：补充denominator的定义）
+    denominator = 1 + R12 * R23 * np.exp(-alpha_abs) + 2 * np.sqrt(R12 * R23 * np.exp(-alpha_abs)) * np.cos(2 * delta)
+
+    # 总反射率
+    numerator_R = R12 + R23 * np.exp(-alpha_abs) + 2 * np.sqrt(R12 * R23 * np.exp(-alpha_abs)) * np.cos(2 * delta)
+    R_total = numerator_R / denominator
+
+    # 总透射率（修复后）
+    T_total = (1 - R12) * (1 - R23) * np.exp(-alpha_abs) / denominator
+
+    # 吸收率=发射率（热平衡下）
+    alpha_total = 1 - R_total - T_total
+    return alpha_total, R_total, T_total  # alpha_total即ε（发射率）
+
+
+# -----------------------------
+# 2. 数据读取与预处理
+# -----------------------------
+# 替换为你的data.txt实际路径（确保正确）
+DATA_PATH = '/Users/spasolreisa/IdeaProjects/asiaMath/data.txt'
+wl_all, n_all, k_all = read_split_data(DATA_PATH)
+wl_all, n_all, k_all = make_strictly_increasing(wl_all, n_all, k_all)
+print(f"数据读取成功：波长范围 {wl_all.min():.2f}–{wl_all.max():.2f} μm，共{len(wl_all)}个数据点")
+
+# 三次样条插值
+cs_n = CubicSpline(wl_all, n_all)
+cs_k = CubicSpline(wl_all, k_all)
+
+# 定义PDMS厚度和计算波长范围
+thicknesses = [0.5, 1.0, 1.5, 2.0]
+wl_fine = np.linspace(wl_all.min(), wl_all.max(), 500)  # 细化解析度
+
+# -----------------------------
+# 3. 预计算材料级关键参数（ε和α）
+# -----------------------------
+# 存储平均发射率（8-13μm，黑体辐射加权）和平均太阳吸收率（0.3-2.5μm，太阳光谱加权）
+avg_eps_dict = {}  # 平均发射率 ε_avg
+avg_alpha_dict = {}  # 平均太阳吸收率 α_avg
+
+
+# 定义权重光谱（辐射冷却+太阳吸收关键波段）
+def planck_spectrum(wl, T):
+    """普朗克黑体光谱（8-13μm波段权重）"""
+    wl_m = wl * 1e-6  # 转换为米
+    c1 = 3.7418e8  # 第一辐射常数 (W·μm⁴/m²)
+    c2 = 14388  # 第二辐射常数 (μm·K)
+    return c1 / (wl_m ** 5 * (np.exp(c2 / (wl * T)) - 1))
+
+
+def solar_spectrum_am15(wl):
+    """AM1.5太阳光谱（0.3-2.5μm波段权重）"""
+    spectrum = np.zeros_like(wl)
+    mask = (wl >= 0.3) & (wl <= 2.5)
+    wl_masked = wl[mask]
+    # 经验拟合AM1.5标准光谱
+    spectrum[mask] = np.where(
+        wl_masked < 0.5, 800 + 400 * wl_masked,
+        np.where(wl_masked < 1.0, 1000 - 200 * (wl_masked - 0.5),
+                 np.where(wl_masked < 1.5, 900 - 100 * (wl_masked - 1.0),
+                          750 - 200 * (wl_masked - 1.5)))
+    )
+    return spectrum
+
+
+# 计算各厚度的平均ε和α
+for d in thicknesses:
+    print(f"\n正在计算厚度 {d} μm 的光学性能...")
+
+    # -----------------------------
+    # 计算平均发射率 ε_avg（8-13μm，黑体辐射加权）
+    # -----------------------------
+    if wl_all.min() <= 13 and wl_all.max() >= 8:
+        wl_rad = np.linspace(8, 13, 300)  # 辐射冷却核心波段
+        n_rad = cs_n(wl_rad)
+        k_rad = cs_k(wl_rad)
+        eps_rad, _, _ = thin_film_optical_properties(n_rad, k_rad, d, wl_rad)
+        planck_weight = planck_spectrum(wl_rad, T=298)  # 25℃黑体光谱权重
+        # 加权平均
+        eps_avg = simpson(eps_rad * planck_weight, wl_rad) / simpson(planck_weight, wl_rad)
+    else:
+        print(f"警告：数据未覆盖8-13μm波段，使用全波段平均发射率替代")
+        n_film = cs_n(wl_fine)
+        k_film = cs_k(wl_fine)
+        eps_full, _, _ = thin_film_optical_properties(n_film, k_film, d, wl_fine)
+        eps_avg = np.mean(eps_full)
+    avg_eps_dict[d] = eps_avg
+
+    # -----------------------------
+    # 计算平均太阳吸收率 α_avg（0.3-2.5μm，太阳光谱加权）
+    # -----------------------------
+    if wl_all.min() <= 2.5 and wl_all.max() >= 0.3:
+        wl_solar = np.linspace(0.3, 2.5, 300)  # 太阳光谱波段
+        n_solar = cs_n(wl_solar)
+        k_solar = cs_k(wl_solar)
+        alpha_solar, _, _ = thin_film_optical_properties(n_solar, k_solar, d, wl_solar)
+        solar_weight = solar_spectrum_am15(wl_solar)  # AM1.5太阳光谱权重
+        # 加权平均
+        alpha_avg = simpson(alpha_solar * solar_weight, wl_solar) / simpson(solar_weight, wl_solar)
+    else:
+        print(f"警告：数据未覆盖0.3-2.5μm太阳波段，使用PDMS典型值α=0.08")
+        alpha_avg = 0.08  # PDMS在太阳波段的典型吸收率（低吸收）
+    avg_alpha_dict[d] = alpha_avg
+
+# -----------------------------
+# 4. 系统级：净冷却功率计算（解答思路核心逻辑）
+# -----------------------------
+# 系统物理参数（可根据实际场景调整）
+sigma = 5.67e-8  # 斯特藩-玻尔兹曼常数 (W/m²·K⁴)
+G_sun_list = [500, 700, 900, 1100]  # 不同太阳辐照强度（对应多云到晴天）
+T_amb_list = np.linspace(293, 318, 6)  # 环境温度（20-45℃，转换为开尔文）
+v_wind = 1.5  # 风速 (m/s)
+h_conv = 5.6 + 3.1 * v_wind  # 对流换热系数（经验公式，W/m²·K）
+
+
+def net_cooling_power(eps, alpha, T_s, T_amb, G_sun, h_conv, sigma):
+    """净冷却功率公式：P_net = 辐射散热 - 太阳吸收热 - 对流换热损失"""
+    # 辐射散热（向宇宙太空）
+    rad散热 = eps * sigma * T_s ** 4
+    # 太阳吸收热（从太阳光获取的热量）
+    solar吸热 = alpha * G_sun
+    # 对流换热损失（向环境散热/吸热）
+    conv损失 = h_conv * (T_s - T_amb)
+    # 环境辐射吸收（从环境获取的辐射热）
+    amb_rad吸热 = eps * sigma * T_amb ** 4
+    # 净冷却功率（正值表示主动冷却，负值表示吸热）
+    return rad散热 - solar吸热 - conv损失 - amb_rad吸热
+
+
+def solve_surface_temperature(eps, alpha, T_amb, G_sun, h_conv, sigma):
+    """迭代求解PDMS薄膜表面温度T_s（净冷却功率=0时的热平衡温度）"""
+    T_s_guess = T_amb - 5  # 初始猜测（比环境低5℃）
+    tol = 1e-3  # 收敛精度
+    max_iter = 100  # 最大迭代次数
+    for _ in range(max_iter):
+        P_net = net_cooling_power(eps, alpha, T_s_guess, T_amb, G_sun, h_conv, sigma)
+        # 数值微分求导（牛顿迭代法，确保收敛）
+        dP_dT = (net_cooling_power(eps, alpha, T_s_guess + 1e-4, T_amb, G_sun, h_conv, sigma) -
+                 net_cooling_power(eps, alpha, T_s_guess - 1e-4, T_amb, G_sun, h_conv, sigma)) / (2e-4)
+        if abs(dP_dT) < 1e-6:
+            break  # 避免除以零
+        # 更新猜测值
+        T_s_new = T_s_guess - P_net / dP_dT
+        # 限制温度范围（物理合理值）
+        T_s_new = max(250, min(T_amb + 5, T_s_new))
+        # 检查收敛
+        if abs(T_s_new - T_s_guess) < tol:
+            return T_s_new
+        T_s_guess = T_s_new
+    return T_s_guess  # 若未收敛，返回最后一次猜测值
+
+
+# -----------------------------
+# 5. 全链条分析：材料性能→系统冷却性能
+# -----------------------------
+# 存储各厚度的系统级结果
+system_results = {}
+for d in thicknesses:
+    eps = avg_eps_dict[d]
+    alpha = avg_alpha_dict[d]
+    # 初始化结果矩阵（太阳辐照×环境温度）
+    P_net_matrix = np.zeros((len(G_sun_list), len(T_amb_list)))
+    T_s_matrix = np.zeros((len(G_sun_list), len(T_amb_list)))
+    # 遍历所有太阳辐照和环境温度组合
+    for i, G_sun in enumerate(G_sun_list):
+        for j, T_amb in enumerate(T_amb_list):
+            # 求解表面温度
+            T_s = solve_surface_temperature(eps, alpha, T_amb, G_sun, h_conv, sigma)
+            T_s_matrix[i, j] = T_s
+            # 计算净冷却功率
+            P_net = net_cooling_power(eps, alpha, T_s, T_amb, G_sun, h_conv, sigma)
+            P_net_matrix[i, j] = P_net
+    # 存储结果
+    system_results[d] = {
+        "eps_avg": eps,
+        "alpha_avg": alpha,
+        "P_net": P_net_matrix,
+        "T_s": T_s_matrix
+    }
+
+# -----------------------------
+# 6. 结果可视化（全链条分析图表）
+# -----------------------------
+plt.rcParams['font.sans-serif'] = ['Arial']  # 统一字体
+plt.rcParams['axes.unicode_minus'] = False  # 支持负号
+
+# 图1：材料级性能（平均ε和α随厚度变化）
+fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(12, 5))
+thicknesses_arr = np.array(thicknesses)
+
+# 平均发射率
+ax1.bar(thicknesses_arr - 0.08, [system_results[d]["eps_avg"] for d in thicknesses],
+        width=0.15, label='Avg Emissivity (8-13μm)', color='darkred', alpha=0.8)
+ax1.set_xlabel('PDMS Thickness (μm)', fontsize=12)
+ax1.set_ylabel('Emissivity', fontsize=12)
+ax1.set_title('Average Emissivity (Radiative Cooling Window)', fontsize=14, fontweight='bold')
+ax1.grid(True, alpha=0.3)
+ax1.set_ylim(0, 1.05)
+
+# 平均太阳吸收率
+ax2.bar(thicknesses_arr - 0.08, [system_results[d]["alpha_avg"] for d in thicknesses],
+        width=0.15, label='Avg Solar Absorptivity (0.3-2.5μm)', color='darkblue', alpha=0.8)
+ax2.set_xlabel('PDMS Thickness (μm)', fontsize=12)
+ax2.set_ylabel('Absorptivity', fontsize=12)
+ax2.set_title('Average Solar Absorptivity', fontsize=14, fontweight='bold')
+ax2.grid(True, alpha=0.3)
+ax2.set_ylim(0, 0.2)  # 限制范围，更清晰
+
+plt.tight_layout()
+plt.savefig('material_performance.png', dpi=300, bbox_inches='tight')
+
+# 图2：系统级性能（净冷却功率随环境温度变化，选最优厚度）
+# 最优厚度：ε/α比值最大（平衡高发射和低吸收）
+optimal_d = max(thicknesses, key=lambda x: system_results[x]["eps_avg"] / (system_results[x]["alpha_avg"] + 0.01))
+print(
+    f"\n最优厚度：{optimal_d} μm（ε={system_results[optimal_d]['eps_avg']:.4f}, α={system_results[optimal_d]['alpha_avg']:.4f}）")
+
+fig, ax = plt.subplots(figsize=(12, 6))
+T_amb_c = T_amb_list - 273.15  # 转换为摄氏度
+colors = ['red', 'orange', 'green', 'blue']
+
+for i, G_sun in enumerate(G_sun_list):
+    P_net = system_results[optimal_d]["P_net"][i, :]
+    ax.plot(T_amb_c, P_net, marker='o', markersize=6, linewidth=2,
+            color=colors[i], label=f'Solar Irradiance = {G_sun} W/m²')
+
+ax.set_xlabel('Ambient Temperature (°C)', fontsize=12)
+ax.set_ylabel('Net Cooling Power (W/m²)', fontsize=12)
+ax.set_title(f'Net Cooling Power vs Ambient Temperature (PDMS Thickness = {optimal_d} μm)',
+             fontsize=14, fontweight='bold')
+ax.grid(True, alpha=0.3)
+ax.legend(fontsize=11)
+# 添加零线（区分冷却/吸热）
+ax.axhline(y=0, color='black', linestyle='--', alpha=0.5, label='Zero Cooling Power')
+plt.tight_layout()
+plt.savefig('net_cooling_power.png', dpi=300, bbox_inches='tight')
+
+# 图3：表面温度随环境温度变化
+fig, ax = plt.subplots(figsize=(12, 6))
+for i, G_sun in enumerate(G_sun_list):
+    T_s = system_results[optimal_d]["T_s"][i, :] - 273.15  # 转换为摄氏度
+    ax.plot(T_amb_c, T_s, marker='s', markersize=6, linewidth=2,
+            color=colors[i], label=f'Solar Irradiance = {G_sun} W/m²')
+
+ax.set_xlabel('Ambient Temperature (°C)', fontsize=12)
+ax.set_ylabel('PDMS Surface Temperature (°C)', fontsize=12)
+ax.set_title(f'Surface Temperature vs Ambient Temperature (PDMS Thickness = {optimal_d} μm)',
+             fontsize=14, fontweight='bold')
+ax.grid(True, alpha=0.3)
+ax.legend(fontsize=11)
+# 添加环境温度参考线（y=x）
+ax.plot(T_amb_c, T_amb_c, color='black', linestyle='--', alpha=0.5, label='Ambient Temperature')
+plt.tight_layout()
+plt.savefig('surface_temperature.png', dpi=300, bbox_inches='tight')
+
+plt.show()
+
+# -----------------------------
+# 7. 关键结果输出（量化分析）
+# -----------------------------
+print("\n" + "=" * 60)
+print("材料-系统全链条关键结果")
+print("=" * 60)
+for d in thicknesses:
+    print(f"\n厚度 {d} μm:")
+    print(f"  - 平均发射率（8-13μm）: {system_results[d]['eps_avg']:.4f}")
+    print(f"  - 平均太阳吸收率（0.3-2.5μm）: {system_results[d]['alpha_avg']:.4f}")
+    print(f"  - 最优工况净冷却功率（T_amb=30℃, G_sun=900 W/m²）: {system_results[d]['P_net'][2, 2]:.2f} W/m²")
+    print(f"  - 对应表面温度: {system_results[d]['T_s'][2, 2] - 273.15:.2f} ℃")
+
+print("\n" + "=" * 60)
+print("结论：PDMS薄膜的最优厚度为 {} μm，在典型工况下（30℃环境、900 W/m²太阳辐照）".format(optimal_d))
+print("可实现 {:.2f} W/m² 的净冷却功率，表面温度比环境低 {:.2f} ℃".format(
+    system_results[optimal_d]['P_net'][2, 2],
+    30 - (system_results[optimal_d]['T_s'][2, 2] - 273.15)
+))
+print("=" * 60)
\ No newline at end of file
diff --git a/org/chatgpt2/q1_anaylis.py b/org/chatgpt2/q1_anaylis.py
new file mode 100644
index 0000000..e5f6d61
--- /dev/null
+++ b/org/chatgpt2/q1_anaylis.py
@@ -0,0 +1,329 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy.interpolate import CubicSpline
+from scipy.stats import pearsonr
+from sklearn.metrics import mean_squared_error
+import warnings
+
+from sklearn.metrics.pairwise import cosine_similarity
+
+from org.chatgpt2.q1_2 import wl_max, wl_min
+
+warnings.filterwarnings('ignore')
+
+
+# -----------------------------
+# 1. 通用工具函数（保持不变）
+# -----------------------------
+def make_strictly_increasing(wl, n, k):
+    unique_wl, indices = np.unique(wl, return_index=True)
+    if len(unique_wl) != len(wl):
+        print(f"Removed {len(wl) - len(unique_wl)} duplicate wavelength points")
+        wl, n, k = wl[indices], n[indices], k[indices]
+    is_increasing = np.diff(wl) > 0
+    if not all(is_increasing):
+        valid_indices = np.concatenate([[True], is_increasing])
+        wl, n, k = wl[valid_indices], n[valid_indices], k[valid_indices]
+    return wl, n, k
+
+
+def read_split_data(file_path):
+    with open(file_path, 'r', encoding='utf-8') as f:
+        lines = [line.strip() for line in f if line.strip() and not line.startswith('#')]
+    split_idx = None
+    for i, line in enumerate(lines):
+        if line == 'wl	k':
+            split_idx = i
+            break
+    n_lines = lines[1:split_idx]
+    wl_n, n_list = [], []
+    for line in n_lines:
+        wl, n_val = line.split()
+        wl_n.append(float(wl)), n_list.append(float(n_val))
+    k_lines = lines[split_idx + 1:]
+    wl_k, k_list = [], []
+    for line in k_lines:
+        wl, k_val = line.split()
+        wl_k.append(float(wl)), k_list.append(float(k_val))
+    wl_n, n_list = np.array(wl_n), np.array(n_list)
+    wl_k, k_list = np.array(wl_k), np.array(k_list)
+    assert np.allclose(wl_n, wl_k), "n和k的波长列表不一致！"
+    sorted_idx = np.argsort(wl_n)
+    return wl_n[sorted_idx], n_list[sorted_idx], k_list[sorted_idx]
+
+
+def fresnel_reflectance(n1, k1, n2, k2):
+    m1 = n1 + 1j * k1
+    m2 = n2 + 1j * k2
+    return np.abs((m1 - m2) / (m1 + m2)) ** 2
+
+
+def thin_film_emissivity(n_film, k_film, d, wl, n_air=1.0, k_air=0.0, n_sub=1.5, k_sub=0.0, r=0.0):
+    m_air = n_air + 1j * k_air
+    m_film = n_film + 1j * k_film
+    m_sub = n_sub + 1j * k_sub
+    R12 = np.abs((m_air - m_film) / (m_air + m_film)) ** 2
+    R23 = np.abs((m_film - m_sub) / (m_film + m_sub)) ** 2
+    beta = 2 * np.pi * m_film / wl
+    delta_complex = beta * d
+    alpha = 2 * np.imag(delta_complex)
+    sqrt_term = np.sqrt(R12 * R23 * np.exp(-alpha))
+    cos_term = np.cos(2 * np.real(delta_complex))
+    R_specular = (R12 + R23 * np.exp(-alpha) + 2 * sqrt_term * cos_term) / (
+                1 + R12 * R23 * np.exp(-alpha) + 2 * sqrt_term * cos_term)
+    R_diffuse = 0.05 * r
+    R_total = (1 - r) * R_specular + r * R_diffuse
+    T_total = (1 - R12) * (1 - R23) * np.exp(-alpha) / (1 + R12 * R23 * np.exp(-alpha) + 2 * sqrt_term * cos_term)
+    return np.clip(1 - R_total - T_total, 0, 1)
+
+
+# -----------------------------
+# 2. 新增：局部相似度计算函数（滑动窗口法）
+# -----------------------------
+def calculate_local_similarity(eps1, eps2, wl_common, window_size=0.5):
+    """
+    计算逐波长的局部相似度（滑动窗口法）
+    输入：
+        eps1/eps2: 两个发射率数组
+        wl_common: 统一波长数组
+        window_size: 滑动窗口宽度（μm），默认0.5μm（覆盖5个采样点，保证平滑）
+    输出：
+        local_corr: 逐波长皮尔逊相关系数（局部相似度）
+        local_cos_sim: 逐波长余弦相似度
+        local_mae: 逐波长局部平均绝对误差（相似度的互补指标）
+    """
+    n_points = len(wl_common)
+    local_corr = np.zeros(n_points)
+    local_cos_sim = np.zeros(n_points)
+    local_mae = np.zeros(n_points)
+
+    # 滑动窗口计算局部相似度（窗口中心为每个波长点）
+    for i in range(n_points):
+        # 确定当前窗口的波长范围
+        wl_center = wl_common[i]
+        window_mask = (wl_common >= wl_center - window_size / 2) & (wl_common <= wl_center + window_size / 2)
+        if np.sum(window_mask) < 3:  # 窗口内至少3个点才计算（避免统计无意义）
+            local_corr[i] = np.nan
+            local_cos_sim[i] = np.nan
+            local_mae[i] = np.nan
+            continue
+
+        # 提取窗口内的发射率数据
+        eps1_window = eps1[window_mask]
+        eps2_window = eps2[window_mask]
+
+        # 计算窗口内的相似度指标
+        corr, _ = pearsonr(eps1_window, eps2_window)
+        cos_sim = cosine_similarity(eps1_window.reshape(1, -1), eps2_window.reshape(1, -1))[0][0]
+        mae = np.mean(np.abs(eps1_window - eps2_window))
+
+        # 存储结果（相关系数和余弦相似度归一化到0-1范围）
+        local_corr[i] = (corr + 1) / 2  # 原始相关系数-1~1 → 映射为0~1
+        local_cos_sim[i] = cos_sim
+        local_mae[i] = mae
+
+    return local_corr, local_cos_sim, local_mae
+
+
+# -----------------------------
+# 3. 核心相似性分析函数（新增局部相似度计算）
+# -----------------------------
+def analyze_spectral_similarity(file1, file2, thicknesses=[1.0], n_sub=1.5, k_sub=0.0, r=0.0, window_size=0.5):
+    # Step 1: 数据读取与预处理
+    wl1, n1, k1 = read_split_data(file1)
+    wl2, n2, k2 = read_split_data(file2)
+    wl1, n1, k1 = make_strictly_increasing(wl1, n1, k1)
+    wl2, n2, k2 = make_strictly_increasing(wl2, n2, k2)
+
+    # Step 2: 统一波长范围
+    wl_min = max(wl1.min(), wl2.min())
+    wl_max = min(wl1.max(), wl2.max())
+    if wl_min >= wl_max:
+        raise ValueError("两个文件的波长范围无交集，无法进行相似性分析！")
+    wl_common = np.linspace(wl_min, wl_max, 1000)
+
+    # Step 3: 插值得到统一波长下的n和k
+    cs_n1 = CubicSpline(wl1, n1)
+    cs_k1 = CubicSpline(wl1, k1)
+    cs_n2 = CubicSpline(wl2, n2)
+    cs_k2 = CubicSpline(wl2, k2)
+    n1_common = cs_n1(wl_common)
+    k1_common = cs_k1(wl_common)
+    n2_common = cs_n2(wl_common)
+    k2_common = cs_k2(wl_common)
+
+    # Step 4: 计算发射率和局部相似度
+    emissivity_dict = {}
+    local_similarity_dict = {}
+    for d in thicknesses:
+        eps1 = thin_film_emissivity(n1_common, k1_common, d, wl_common, n_sub=n_sub, k_sub=k_sub, r=r)
+        eps2 = thin_film_emissivity(n2_common, k2_common, d, wl_common, n_sub=n_sub, k_sub=k_sub, r=r)
+        emissivity_dict[d] = (eps1, eps2)
+
+        # 计算局部相似度
+        local_corr, local_cos_sim, local_mae = calculate_local_similarity(eps1, eps2, wl_common, window_size)
+        local_similarity_dict[d] = (local_corr, local_cos_sim, local_mae)
+
+    # Step 5: 全局相似性指标计算
+    similarity_results = {}
+    for d in thicknesses:
+        eps1, eps2 = emissivity_dict[d]
+        pearson_corr, _ = pearsonr(eps1, eps2)
+        cos_sim = cosine_similarity(eps1.reshape(1, -1), eps2.reshape(1, -1))[0][0]
+        mse = mean_squared_error(eps1, eps2)
+        norm_mse = mse / (np.max([np.var(eps1), np.var(eps2)]) + 1e-8)
+        mae = np.mean(np.abs(eps1 - eps2))
+
+        # 大气窗口指标
+        window_corr, window_mae = None, None
+        if wl_min <= 13 and wl_max >= 8:
+            window_mask = (wl_common >= 8) & (wl_common <= 13)
+            eps1_window = eps1[window_mask]
+            eps2_window = eps2[window_mask]
+            window_corr, _ = pearsonr(eps1_window, eps2_window)
+            window_mae = np.mean(np.abs(eps1_window - eps2_window))
+
+        similarity_results[d] = {
+            "pearson_correlation": pearson_corr,
+            "cosine_similarity": cos_sim,
+            "normalized_mse": norm_mse,
+            "mae": mae,
+            "window_pearson_correlation": window_corr,
+            "window_mae": window_mae
+        }
+
+    # Step 6: 可视化（新增相似度曲线）
+    plot_spectral_comparison(wl_common, emissivity_dict, local_similarity_dict, thicknesses, wl_min, wl_max)
+
+    return similarity_results, wl_common, emissivity_dict, local_similarity_dict
+
+
+# -----------------------------
+# 4. 可视化函数（新增相似度曲线子图）
+# -----------------------------
+def plot_spectral_comparison(wl_common, emissivity_dict, local_similarity_dict, thicknesses, wl_min, wl_max):
+    n_plots = len(thicknesses)
+    fig, axes = plt.subplots(n_plots, 3, figsize=(18, 5 * n_plots))  # 新增1列用于相似度曲线
+    plt.rcParams['font.sans-serif'] = ['Arial']
+
+    for idx, d in enumerate(thicknesses):
+        eps1, eps2 = emissivity_dict[d]
+        local_corr, local_cos_sim, local_mae = local_similarity_dict[d]
+        diff = np.abs(eps1 - eps2)
+
+        # 子图1：发射率光谱对比
+        ax1 = axes[idx, 0] if n_plots > 1 else axes[0]
+        ax1.plot(wl_common, eps1, linewidth=2, label='data.txt', color='darkblue')
+        ax1.plot(wl_common, eps2, linewidth=2, label='data2.txt', color='darkred', linestyle='--')
+        if wl_min <= 13 and wl_max >= 8:
+            ax1.axvspan(8, 13, alpha=0.15, color='orange', label='Atmospheric Window (8-13 μm)')
+        ax1.set_xlabel('Wavelength (μm)', fontsize=12)
+        ax1.set_ylabel('Emissivity ε(λ)', fontsize=12)
+        ax1.set_title(f'Emissivity Spectrum Comparison (Thickness = {d} μm)', fontsize=14, fontweight='bold')
+        ax1.grid(True, alpha=0.3)
+        ax1.legend(fontsize=10)
+        ax1.set_ylim(0, 1.05)
+
+        # 子图2：发射率差异
+        ax2 = axes[idx, 1] if n_plots > 1 else axes[1]
+        ax2.plot(wl_common, diff, linewidth=2, color='darkgreen', label='Absolute Difference |ε1 - ε2|')
+        ax2.fill_between(wl_common, 0, diff, alpha=0.3, color='darkgreen')
+        if wl_min <= 13 and wl_max >= 8:
+            ax2.axvspan(8, 13, alpha=0.15, color='orange', label='Atmospheric Window (8-13 μm)')
+        ax2.set_xlabel('Wavelength (μm)', fontsize=12)
+        ax2.set_ylabel('Absolute Difference', fontsize=12)
+        ax2.set_title(f'Emissivity Difference (Thickness = {d} μm)', fontsize=14, fontweight='bold')
+        ax2.grid(True, alpha=0.3)
+        ax2.legend(fontsize=10)
+        ax2.set_ylim(0, np.nanmax(diff) * 1.2)
+
+        # 子图3：相似度曲线（核心新增）
+        ax3 = axes[idx, 2] if n_plots > 1 else axes[2]
+        # 绘制局部皮尔逊相关系数（归一化到0-1）
+        ax3.plot(wl_common, local_corr, linewidth=2.5, color='#4B0082', label='Local Pearson Correlation (0-1)')        # 绘制局部余弦相似度（0-1）
+        ax3.plot(wl_common, local_cos_sim, linewidth=2.5, color='orange', label='Local Cosine Similarity (0-1)',
+                 linestyle='--')
+        # 标注相似度阈值线（0.9为高相似，0.8为中等相似）
+        ax3.axhline(y=0.9, color='red', linestyle=':', linewidth=1.5, label='High Similarity Threshold (0.9)')
+        ax3.axhline(y=0.8, color='orange', linestyle=':', linewidth=1.5, label='Medium Similarity Threshold (0.8)')
+        # 标注大气窗口
+        if wl_min <= 13 and wl_max >= 8:
+            ax3.axvspan(8, 13, alpha=0.15, color='orange', label='Atmospheric Window (8-13 μm)')
+        ax3.set_xlabel('Wavelength (μm)', fontsize=12)
+        ax3.set_ylabel('Local Similarity (0-1)', fontsize=12)
+        ax3.set_title(f'Wavelength-Dependent Similarity Curve (Thickness = {d} μm)', fontsize=14, fontweight='bold')
+        ax3.grid(True, alpha=0.3)
+        ax3.legend(fontsize=10)
+        ax3.set_ylim(0, 1.05)  # 相似度范围0-1
+        ax3.set_xlim(wl_min, wl_max)
+
+    plt.tight_layout()
+    plt.savefig('spectral_similarity_complete.png', dpi=300, bbox_inches='tight')
+    plt.show()
+
+
+# -----------------------------
+# 5. 主程序执行
+# -----------------------------
+if __name__ == "__main__":
+    file1 = '/Users/spasolreisa/IdeaProjects/asiaMath/data.txt'
+    file2 = '/Users/spasolreisa/IdeaProjects/asiaMath/data2.txt'
+    target_thicknesses = [0.5, 1.0, 1.5, 2.0]
+    base_params = {
+        'n_sub': 1.5,
+        'k_sub': 0.0,
+        'r': 0.05
+    }
+    window_size = 0.5  # 滑动窗口宽度（可调整，建议0.3-0.8μm）
+
+    print("=== 开始光谱相似性分析（含相似度曲线） ===")
+    print(f"文件1: {file1}")
+    print(f"文件2: {file2}")
+    print(f"分析厚度: {target_thicknesses} μm")
+    print(f"滑动窗口宽度: {window_size} μm")
+    print("-" * 50)
+
+    results, wl_common, emissivity_dict, local_similarity_dict = analyze_spectral_similarity(
+        file1, file2,
+        thicknesses=target_thicknesses,
+        n_sub=base_params['n_sub'],
+        k_sub=base_params['k_sub'],
+        r=base_params['r'],
+        window_size=window_size
+    )
+
+    # 输出全局指标
+    print("\n=== 全局相似性指标 ===")
+    for d in target_thicknesses:
+        res = results[d]
+        print(f"\n【厚度 {d} μm】")
+        print(f"全局皮尔逊相关系数: {res['pearson_correlation']:.4f}")
+        print(f"全局余弦相似度: {res['cosine_similarity']:.4f}")
+        print(f"归一化均方误差: {res['normalized_mse']:.4f}")
+        print(f"平均绝对误差: {res['mae']:.4f}")
+        if res['window_pearson_correlation'] is not None:
+            print(f"大气窗口全局相关系数: {res['window_pearson_correlation']:.4f}")
+
+    # 输出局部相似度统计（大气窗口内）
+    print("\n=== 大气窗口（8-13μm）局部相似度统计 ===")
+    for d in target_thicknesses:
+        local_corr, local_cos_sim, _ = local_similarity_dict[d]
+        if wl_min <= 13 and wl_max >= 8:
+            window_mask = (wl_common >= 8) & (wl_common <= 13)
+            window_local_corr = local_corr[window_mask]
+            window_local_cos = local_cos_sim[window_mask]
+            # 过滤NaN值
+            window_local_corr = window_local_corr[~np.isnan(window_local_corr)]
+            window_local_cos = window_local_cos[~np.isnan(window_local_cos)]
+            if len(window_local_corr) > 0:
+                print(f"\n【厚度 {d} μm】")
+                print(f"大气窗口局部相关系数均值: {np.mean(window_local_corr):.4f}")
+                print(f"大气窗口局部相关系数最小值: {np.min(window_local_corr):.4f}")
+                print(f"大气窗口局部余弦相似度均值: {np.mean(window_local_cos):.4f}")
+
+    # 结果解读
+    print("\n=== 结果解读 ===")
+    print("1. 相似度曲线（子图3）：值越接近1，对应波长下的发射率越相似；")
+    print("2. 高相似区域（≥0.9）：两文件在该波段的辐射冷却性能几乎一致；")
+    print("3. 低相似区域（<0.8）：需关注该波段的材料差异对冷却效果的影响；")
+    print("4. 大气窗口内相似度：优先关注该区域，直接决定辐射冷却核心性能是否一致。")
\ No newline at end of file
diff --git a/org/chatgpt2/q2_2.py b/org/chatgpt2/q2_2.py
new file mode 100644
index 0000000..475af45
--- /dev/null
+++ b/org/chatgpt2/q2_2.py
@@ -0,0 +1,321 @@
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy.interpolate import CubicSpline
+from scipy.integrate import simpson
+import os
+
+# -----------------------------
+# Configuration (Update File Path!)
+# -----------------------------
+DATA_FILE_PATH = "/Users/spasolreisa/IdeaProjects/asiaMath/data.txt"  # Replace with your data.txt absolute path
+THICKNESSES = [0.5, 1.0, 1.5, 2.0, 2.5, 3.0]  # Expand thickness range for evaluation
+T_AMBIENT = 300  # Ambient temperature (K)
+SOLAR_IRRADIANCE = 1000  # AM1.5 solar irradiance (W/m²)
+CONVECTION_COEFF = 10  # Convection coefficient (W/(m²K))
+SIGMA = 5.67e-8  # Stefan-Boltzmann constant (W/(m²K⁴))
+
+
+# -----------------------------
+# 1. Fixed Data Parsing Function (Critical Fix for "wl" String Error)
+# -----------------------------
+def read_split_data(file_path):
+    """Read and parse split-format data (wl+n followed by wl+k)"""
+    if not os.path.exists(file_path):
+        raise FileNotFoundError(f"File not found: {file_path}")
+
+    # Read all lines, skip empty lines and comments
+    with open(file_path, 'r', encoding='utf-8') as f:
+        lines = []
+        for line in f:
+            stripped = line.strip()
+            if stripped and not stripped.startswith('#'):
+                lines.append(stripped)
+
+    # Step 1: Identify all headers (lines containing "wl" and either "n" or "k")
+    header_indices = []
+    for i, line in enumerate(lines):
+        parts = line.split()
+        # Header must be exactly 2 parts: ["wl", "n"] or ["wl", "k"] (case-insensitive)
+        if len(parts) == 2 and parts[0].lower() == "wl" and parts[1].lower() in ["n", "k"]:
+            header_indices.append(i)
+
+    # Validate: Must have exactly 2 headers (one for n, one for k)
+    if len(header_indices) != 2:
+        raise ValueError(
+            f"Invalid number of headers! Expected 2 (wl+n and wl+k), found {len(header_indices)}.\nCheck data.txt format.")
+
+    # Step 2: Split data into n-block and k-block
+    n_header_idx = header_indices[0]
+    k_header_idx = header_indices[1]
+
+    # Ensure n-header comes before k-header
+    if n_header_idx > k_header_idx:
+        n_header_idx, k_header_idx = k_header_idx, n_header_idx
+
+    # Extract n data (between n-header and k-header)
+    n_lines = lines[n_header_idx + 1: k_header_idx]
+    # Extract k data (after k-header)
+    k_lines = lines[k_header_idx + 1:]
+
+    # Step 3: Parse n data (skip any invalid lines)
+    wl_n, n_list = [], []
+    for line in n_lines:
+        parts = line.split()
+        # Data line must have exactly 2 numeric parts
+        if len(parts) != 2:
+            continue  # Skip lines with wrong column count
+        try:
+            wl_val = float(parts[0])
+            n_val = float(parts[1])
+            wl_n.append(wl_val)
+            n_list.append(n_val)
+        except ValueError:
+            continue  # Skip non-numeric lines
+
+    # Step 4: Parse k data (skip any invalid lines)
+    wl_k, k_list = [], []
+    for line in k_lines:
+        parts = line.split()
+        if len(parts) != 2:
+            continue
+        try:
+            wl_val = float(parts[0])
+            k_val = float(parts[1])
+            wl_k.append(wl_val)
+            k_list.append(k_val)
+        except ValueError:
+            continue
+
+    # Validate: Must have at least 1 data point for n and k
+    if len(wl_n) == 0:
+        raise ValueError("No valid n data found! Check the format between wl+n and wl+k headers.")
+    if len(wl_k) == 0:
+        raise ValueError("No valid k data found! Check the format after wl+k header.")
+
+    # Convert to numpy arrays
+    wl_n, n_list = np.array(wl_n), np.array(n_list)
+    wl_k, k_list = np.array(wl_k), np.array(k_list)
+
+    # Align wavelengths (if n and k have different wavelength points)
+    if not np.allclose(wl_n, wl_k, rtol=1e-6):
+        print("Warning: Wavelengths for n and k do not match. Automatically aligning...")
+        # Use n's wavelengths as reference, interpolate k to match
+        k_list = np.interp(wl_n, np.sort(wl_k), k_list[np.argsort(wl_k)])
+        wl_k = wl_n  # Sync k's wavelengths to n's
+
+    # Sort by wavelength (ascending) to avoid interpolation errors
+    sorted_idx = np.argsort(wl_n)
+    sorted_wl = wl_n[sorted_idx]
+    sorted_n = n_list[sorted_idx]
+    sorted_k = k_list[sorted_idx]
+
+    print(f"Data loaded successfully: {len(sorted_wl)} valid wavelength points")
+    print(f"Wavelength range: {sorted_wl.min():.2f}–{sorted_wl.max():.2f} μm")
+    return sorted_wl, sorted_n, sorted_k
+
+
+# -----------------------------
+# 2. Core Functions (Unchanged)
+# -----------------------------
+def planck_function(wl, T):
+    """Planck's law: Blackbody radiation (W/(m³sr))"""
+    wl_m = wl * 1e-6  # Convert μm to m
+    c1 = 3.7418e8  # First radiation constant (Wμm⁴/m²)
+    c2 = 14388  # Second radiation constant (μmK)
+    return c1 / (wl_m ** 5 * (np.exp(c2 / (wl * T)) - 1))
+
+
+def solar_spectrum_am15(wl):
+    """AM1.5 global solar irradiance (W/(m²μm))"""
+    spectrum = np.zeros_like(wl)
+    mask = (wl >= 0.3) & (wl <= 2.5)
+    wl_masked = wl[mask]
+    # Empirical fit to AM1.5 data (valid for 0.3–2.5 μm)
+    spectrum[mask] = np.where(
+        wl_masked < 0.5, 800 + 400 * wl_masked,
+        np.where(wl_masked < 1.0, 1000 - 200 * (wl_masked - 0.5),
+                 np.where(wl_masked < 1.5, 900 - 100 * (wl_masked - 1.0),
+                          750 - 200 * (wl_masked - 1.5)))
+    )
+    return spectrum
+
+
+def fresnel_reflectance(n1, k1, n2, k2):
+    """Fresnel reflectance (normal incidence, complex refractive index)"""
+    m1, m2 = n1 + 1j * k1, n2 + 1j * k2
+    return np.abs((m1 - m2) / (m1 + m2)) ** 2
+
+
+def thin_film_optical_properties(n_film, k_film, d, wl):
+    """Calculate emissivity (ε), absorptivity (α), transmissivity (T) of thin film"""
+    R12 = fresnel_reflectance(1.0, 0.0, n_film, k_film)  # Air→Film
+    R23 = fresnel_reflectance(n_film, k_film, 1.0, 0.0)  # Film→Air
+    delta = 2 * np.pi * n_film * d / wl  # Phase difference
+    alpha_abs = 4 * np.pi * k_film * d / wl  # Absorption attenuation
+
+    # Total reflectance and transmissivity (multiple-beam interference)
+    R_total = (R12 + R23 * np.exp(-alpha_abs) + 2 * np.sqrt(R12 * R23 * np.exp(-alpha_abs)) * np.cos(2 * delta)) / \
+              (1 + R12 * R23 * np.exp(-alpha_abs) + 2 * np.sqrt(R12 * R23 * np.exp(-alpha_abs)) * np.cos(2 * delta))
+    T_total = (1 - R12) * (1 - R23) * np.exp(-alpha_abs) / \
+              (1 + R12 * R23 * np.exp(-alpha_abs) + 2 * np.sqrt(R12 * R23 * np.exp(-alpha_abs)) * np.cos(2 * delta))
+    alpha_total = 1 - R_total - T_total  # Kirchhoff's law (α=ε for thermal equilibrium)
+    return alpha_total, R_total, T_total  # α=ε for emissivity
+
+
+# -----------------------------
+# 3. Evaluation Model (Unchanged)
+# -----------------------------
+def evaluate_radiative_cooling(wl_all, n_all, k_all, thickness):
+    """Calculate KPIs and comprehensive score for a given PDMS thickness"""
+    cs_n = CubicSpline(wl_all, n_all)
+    cs_k = CubicSpline(wl_all, k_all)
+
+    # KPI 1: Average Emissivity in 8–13 μm (weighted by Planck function)
+    wl_window = np.linspace(8, 13, 500)
+    # Check if data covers the window (otherwise use nearest values)
+    if wl_all.min() > 8 or wl_all.max() < 13:
+        print(f"Warning: Data does not fully cover 8–13 μm window. Extrapolating...")
+        n_window = cs_n(wl_window, extrapolate=True)
+        k_window = cs_k(wl_window, extrapolate=True)
+    else:
+        n_window = cs_n(wl_window)
+        k_window = cs_k(wl_window)
+    eps_window, _, _ = thin_film_optical_properties(n_window, k_window, thickness, wl_window)
+    planck = planck_function(wl_window, T_AMBIENT)
+    eps_avg = simpson(eps_window * planck, wl_window) / simpson(planck, wl_window)
+
+    # KPI 2: Average Solar Absorptivity in 0.3–2.5 μm (weighted by AM1.5)
+    wl_solar = np.linspace(0.3, 2.5, 500)
+    if wl_all.min() > 2.5 or wl_all.max() < 0.3:
+        print(f"Warning: Data does not cover solar spectrum (0.3–2.5 μm). Using default PDMS properties...")
+        n_solar = np.ones_like(wl_solar) * 1.4  # Typical PDMS n in solar range
+        k_solar = np.ones_like(wl_solar) * 1e-6  # Typical PDMS k in solar range
+    else:
+        n_solar = cs_n(wl_solar, extrapolate=True)
+        k_solar = cs_k(wl_solar, extrapolate=True)
+    alpha_solar, _, _ = thin_film_optical_properties(n_solar, k_solar, thickness, wl_solar)
+    solar_irr = solar_spectrum_am15(wl_solar)
+    alpha_avg = simpson(alpha_solar * solar_irr, wl_solar) / simpson(solar_irr, wl_solar)
+
+    # KPI 3: Maximum Cooling Temperature (ΔT_max)
+    def heat_flux(T_film):
+        planck_film = planck_function(wl_window, T_film)
+        eps_eff = simpson(eps_window * planck_film, wl_window) / simpson(planck_film, wl_window)
+        return SIGMA * eps_eff * T_film ** 4 - alpha_avg * SOLAR_IRRADIANCE - CONVECTION_COEFF * (T_film - T_AMBIENT)
+
+    # Newton-Raphson iteration (stable convergence)
+    T_film = T_AMBIENT - 10  # Initial guess
+    for _ in range(50):
+        q = heat_flux(T_film)
+        if abs(q) < 1e-3:
+            break
+        # Numerical derivative (more stable than analytical)
+        dq_dT = (heat_flux(T_film + 1e-4) - heat_flux(T_film - 1e-4)) / (2e-4)
+        T_film -= q / dq_dT
+        # Prevent unrealistic temperatures
+        if T_film < 200 or T_film > T_AMBIENT:
+            T_film = max(200, min(T_AMBIENT - 5, T_film))
+    delta_T = T_AMBIENT - T_film
+
+    # KPI 4: Cooling Efficiency Ratio (η_CR)
+    eta_cr = eps_avg / (alpha_avg + 0.01)  # +0.01 to avoid division by zero
+
+    # Comprehensive Score (0–100)
+    score = 0.0
+    score += 40 * min(eps_avg, 1.0)  # Cap at 1.0 (ideal emissivity)
+    score += 35 * (1 - min(alpha_avg, 1.0))  # Lower absorption = higher score
+    score += 15 * min(delta_T / 40, 1.0)  # ΔT theoretical upper limit = 40K
+    score += 10 * min(eta_cr / 100, 1.0)  # Cap at 100 (ideal ratio)
+
+    return {
+        "thickness": thickness,
+        "eps_8-13": eps_avg,
+        "alpha_0.3-2.5": alpha_avg,
+        "delta_T_max": delta_T,
+        "eta_cr": eta_cr,
+        "comprehensive_score": score
+    }
+
+
+# -----------------------------
+# 4. Main Execution (Unchanged)
+# -----------------------------
+if __name__ == "__main__":
+    try:
+        # Read data (fixed parsing logic)
+        wl_all, n_all, k_all = read_split_data(DATA_FILE_PATH)
+        print("\n" + "-" * 50 + "\n")
+
+        # Evaluate each thickness
+        results = []
+        for d in THICKNESSES:
+            res = evaluate_radiative_cooling(wl_all, n_all, k_all, d)
+            results.append(res)
+            print(f"Thickness: {d} μm")
+            print(f"  - Avg Emissivity (8–13 μm): {res['eps_8-13']:.4f}")
+            print(f"  - Avg Solar Absorptivity (0.3–2.5 μm): {res['alpha_0.3-2.5']:.4f}")
+            print(f"  - Max Cooling Temperature: {res['delta_T_max']:.2f} K")
+            print(f"  - Cooling Efficiency Ratio: {res['eta_cr']:.2f}")
+            print(f"  - Comprehensive Score: {res['comprehensive_score']:.1f}/100\n")
+
+        # Convert results to numpy array for plotting
+        results_arr = np.array([[
+            res["thickness"], res["eps_8-13"], res["alpha_0.3-2.5"],
+            res["delta_T_max"], res["comprehensive_score"]
+        ] for res in results])
+
+        # Plot KPIs vs Thickness
+        fig, axes = plt.subplots(2, 2, figsize=(14, 10))
+        fig.suptitle("PDMS Thin Film Radiative Cooling Performance vs Thickness", fontsize=16, fontweight='bold')
+
+        # Emissivity (8–13 μm)
+        axes[0, 0].plot(results_arr[:, 0], results_arr[:, 1], 'o-', color='darkred', linewidth=2, markersize=6)
+        axes[0, 0].set_xlabel("Thickness (μm)", fontsize=12), axes[0, 0].set_ylabel("Avg Emissivity (8–13 μm)",
+                                                                                    fontsize=12)
+        axes[0, 0].grid(True, alpha=0.3), axes[0, 0].set_ylim(0, 1.05)
+
+        # Solar Absorptivity (0.3–2.5 μm)
+        axes[0, 1].plot(results_arr[:, 0], results_arr[:, 2], 's-', color='darkblue', linewidth=2, markersize=6)
+        axes[0, 1].set_xlabel("Thickness (μm)", fontsize=12), axes[0, 1].set_ylabel(
+            "Avg Solar Absorptivity (0.3–2.5 μm)", fontsize=12)
+        axes[0, 1].grid(True, alpha=0.3), axes[0, 1].set_ylim(0, 0.5)
+
+        # Max Cooling Temperature
+        axes[1, 0].plot(results_arr[:, 0], results_arr[:, 3], '^-', color='darkgreen', linewidth=2, markersize=6)
+        axes[1, 0].set_xlabel("Thickness (μm)", fontsize=12), axes[1, 0].set_ylabel("Max Cooling Temperature (K)",
+                                                                                    fontsize=12)
+        axes[1, 0].grid(True, alpha=0.3)
+
+        # Comprehensive Score
+        axes[1, 1].plot(results_arr[:, 0], results_arr[:, 4], 'd-', color='darkorange', linewidth=2, markersize=6)
+        axes[1, 1].set_xlabel("Thickness (μm)", fontsize=12), axes[1, 1].set_ylabel("Comprehensive Score (0–100)",
+                                                                                    fontsize=12)
+        axes[1, 1].grid(True, alpha=0.3), axes[1, 1].set_ylim(0, 100)
+
+        plt.tight_layout()
+        plt.savefig("PDMS_radiative_cooling_evaluation.png", dpi=300, bbox_inches='tight')
+        plt.show()
+
+        # Highlight optimal thickness
+        optimal = max(results, key=lambda x: x["comprehensive_score"])
+        print("=" * 50)
+        print(f"Optimal PDMS Thickness: {optimal['thickness']} μm")
+        print(f"Best Comprehensive Score: {optimal['comprehensive_score']:.1f}/100")
+        print(
+            f"Key Performance: ε(8-13μm)={optimal['eps_8-13']:.4f}, α(0.3-2.5μm)={optimal['alpha_0.3-2.5']:.4f}, ΔT={optimal['delta_T_max']:.2f}K")
+        print("=" * 50)
+
+    except Exception as e:
+        print(f"\nError: {e}")
+        print("\nTroubleshooting Steps:")
+        print("1. Check data.txt format: Ensure it has exactly two headers (e.g., 'wl n' and 'wl k')")
+        print("2. Example valid format:")
+        print("   wl n")
+        print("   0.40 1.41491")
+        print("   0.41 1.41403")
+        print("   ...")
+        print("   wl k")
+        print("   0.40 1.40E-06")
+        print("   0.41 1.38E-06")
+        print("3. Ensure no extra 'wl' strings in data lines (only numbers)")
+        print("4. Use space or tab as separator (avoid commas)")
\ No newline at end of file