yfunction fn=f1(y0)
y=[
130
195
278
283
400
421
451
458
594
603
670
681
683
685
712
713
718
723
739
765
799
799
805
832
834
880
883
896
908
910
923
928
929
947
954
975
1005
1019
1041
1054
1073
1084
1103
1169
1169
1188
1194
1207
1207
1224
1228
1257
1259
1278
1282
1286
1296
1309
1327
1378
1394
1418
1462
1464
1486
1512
1538
1544
1545
1564
1581
1595
1598
1598
1623
1647
1652
1687
1691
1705
1735
1767
1779
1781
1786
1803
1805
1831
1842
1846
1852
1868
1887
1889
1894
1897
1909
1914
1935
1938
1939
1998
2017
2019
2020
2024
2031
2031
2035
2039
2049
2061
2074
2087
2096
2110
2121
2129
2141
2152
2155
2156
2176
2183
2214
2225
2255
2270
2285
2294
2299
2302
2315
2322
2333
2335
2336
2345
2352
2371
2372
2381
2405
2406
2416
2439
2439
2440
2473
2476
2486
2492
2497
2512
2512
2516
2523
2527
2528
2533
2567
2570
2612
2623
2636
2638
2640
2641
2650
2652
2669
2688
2721
2746
2762
2762
2779
2787
2796
2806
2827
2831
2834
2840
2850
2858
2861
2870
2881
2900
2918
2926
2933
2936
2945
2948
2949
2956
2958
2994
2998
2998
3008
3027
3044
3093
3095
3103
3118
3132
3137
3175
3175
3179
3190
3197
3258
3269
3299
3301
3315
3333
3336
3338
3342
3344
3349
3357
3361
3361
3379
3383
3384
3420
3439
3441
3459
3469
3505
3518
3536
3537
3554
3559
3560
3561
3626
3635
3639
3645
3661
3685
3686
3702
3722
3724
3747
3751
3770
3785
3801
3805
3806
3841
3843
3850
3864
3907
3908
3921
3921
3921
3924
3931
3938
3950
3952
3972
4029
4033
4045
4054
4077
4079
4087
4099
4109
4120
4134
4144
4147
4161
4164
4166
4166
4167
4170
4177
4177
4187
4188
4192
4196
4212
4231
4237
4247
4249
4251
4279
4280
4286
4305
4306
4316
4317
4346
4378
4403
4422
4443
4445
4467
4482
4486
4488
4501
4518
4546
4560
4561
4567
4571
4572
4583
4645
4656
4656
4688
4701
4710
4710
4730
4731
4750
4770
4785
4785
4807
4812
4821
4833
4841
4900
4902
4917
4937
4940
4943
4989
5006
5022
5039
5082
5082
5128
5145
5168
5183
5206
5224
5231
5240
5256
5272
5315
5323
5324
5330
5346
5356
5381
5407
5424
5426
5457
5469
5476
5508
5539
5541
5550
5564
5566
5587
5636
5640
5703
5737
5744
5835
5841
5930
5943
5967
5995
6006
6031
6099
6134
6150
6196
6249
6261
6263
6306
6309
6313
6315
6318
6387
6430
6431
6440
6464
6487
6550
6589
6716
6744
6819
6861
7255
7282
7368
7445
7509
7530
7538
7572
7583
7644
7716
7803
7861
7927
7961
7978
7990
8006
8103
8280
8396
8420
8422
8431
8479
8520
8560
8658
8676
8695
8894
8911
8937
9016
9105
9413
9582
9612
9672
9826
9929
10233
10301
10500
10784
11242
11424
11481
11640
12250
12581
12665
13399
13685
13687
14370
14638
15389
15418
15964
16873
17299
20362
20469
20473
24235
25092
30038
31243
31485
33284
34000
36789
40901
41185
45694
50000];
f=[];n=509;hn=707;
for i=1:1:509;
    f(i)=(1/(2*pi)^(1/2))*exp(-((y0-y(i))./hn).^2/2);
end
fn=1/(n*hn)*sum(f(i));怎么求fn的积分