├── README.md ├── curves ├── BCPF_fiber_rmselb_m30_r10.csv ├── BCPF_fiber_rmselb_m30_r5.csv ├── mu10.mat ├── mu_curve10.m ├── mu_curve10.pdf ├── mu_curve10.png ├── overfitting.m ├── overfitting_ms30_r10.pdf ├── overfitting_ms30_r10.png ├── overfitting_ms30_r5.pdf └── overfitting_ms30_r5.png ├── heat-maps ├── bias10.mat ├── bias20.mat ├── bias60.mat ├── factor2.pdf ├── factor2.png ├── factor3.pdf ├── factor3.png ├── heat_map10.m ├── heat_map10.pdf ├── heat_map10.png ├── heat_map20.m ├── heat_map20.pdf ├── heat_map60.m ├── heat_map60.pdf ├── latent_factors.m └── latent_factors.mat ├── rmse-curves ├── rmse.mat ├── rmse10.mat ├── rmse_curve.pdf ├── rmse_curve.png ├── rmse_curve10.m ├── rmse_curve10.pdf ├── rmse_curve10.png └── rmse_curves.m └── time-series ├── nyc_data_competness.xls ├── nyc_data_completeness.m ├── nyc_data_completeness.mat ├── nyc_data_completeness.pdf ├── nyc_data_completeness.png ├── road1_fiber_ms50_r10.csv ├── speed_curve.m ├── speed_curve.pdf ├── speed_curve.png ├── time_series_speed1.pdf ├── time_series_speed1.png ├── time_series_speed1.py ├── time_series_speed2.pdf ├── time_series_speed2.png └── time_series_speed2.py /README.md: -------------------------------------------------------------------------------- 1 | Academic drawing 2 | ----------------- 3 | 4 | This is a project providing source codes (including Matlab and Python) for presenting experiment results. 5 | 6 | Contents 7 | -------- 8 | 9 | - [Usage](#usage) 10 | - [Our examples](#our-examples) 11 | - [Our publications](#our-publications) 12 | 13 | Usage 14 | -------------- 15 | 16 | > It is not necessary to open each file in this repository because you can follow this readme document to find your needs. 17 | 18 | Our examples 19 | -------------- 20 | 21 | - **Download** 22 | - [mu10.mat](https://github.com/xinychen/academic-drawing/blob/master/curves/mu10.mat) 23 | - [mu_curve10.m](https://github.com/xinychen/academic-drawing/blob/master/curves/mu_curve10.m) 24 | 25 | and evaluate these in Matlab, then, you will see the following picture: 26 | 27 | ![mu_curve10](https://github.com/xinychen/academic-drawing/blob/master/curves/mu_curve10.png) 28 | 29 | 30 | - **Download** 31 | - [BCPF_fiber_rmselb_m30_r5.csv](https://github.com/xinychen/academic-drawing/blob/master/curves/BCPF_fiber_rmselb_m30_r5.csv) 32 | - [BCPF_fiber_rmselb_m30_r10.csv](https://github.com/xinychen/academic-drawing/blob/master/curves/BCPF_fiber_rmselb_m30_r10.csv) 33 | - [overfitting.m](https://github.com/xinychen/academic-drawing/blob/master/curves/overfitting.m) 34 | 35 | and evaluate these in Matlab, then, you will see the following pictures: 36 | 37 | ![overfitting_ms30_r5](https://github.com/xinychen/academic-drawing/blob/master/curves/overfitting_ms30_r5.png) 38 | ![overfitting_ms30_r10](https://github.com/xinychen/academic-drawing/blob/master/curves/overfitting_ms30_r10.png) 39 | 40 | 41 | - **Download** 42 | - [bias10.mat](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/bias10.mat) 43 | - [heat_map10.m](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/heat_map10.m) 44 | 45 | and evaluate these in Matlab, then, you will see the following picture: 46 | 47 | ![heat_map10](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/heat_map10.png) 48 | 49 | 50 | - **Download** 51 | - [latent_factors.mat](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/latent_factors.mat) 52 | - [latent_factors.m](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/latent_factors.m) 53 | 54 | and evaluate these in Matlab, then, you will see the following pictures: 55 | 56 | ![factor2](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/factor2.png) 57 | ![factor3](https://github.com/xinychen/academic-drawing/blob/master/heat-maps/factor3.png) 58 | 59 | 60 | - **Download** 61 | - [rmse.mat](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse.mat) 62 | - [rmse_curves.m](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse_curves.m) 63 | 64 | and evaluate these in Matlab, then, you will see the following picture: 65 | 66 | ![rmse_curve](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse_curve.png) 67 | 68 | 69 | - **Download** 70 | - [rmse10.mat](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse10.mat) 71 | - [rmse_curve10.m](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse_curve10.m) 72 | 73 | and evaluate these in Matlab, then, you will see the following picture: 74 | 75 | ![rmse_curve10](https://github.com/xinychen/academic-drawing/blob/master/rmse-curves/rmse_curve10.png) 76 | 77 | 78 | - **Download** 79 | - [road1_fiber_ms50_r10.csv](https://github.com/xinychen/academic-drawing/blob/master/time-series/road1_fiber_ms50_r10.csv) 80 | - [time_series_speed1.py](https://github.com/xinychen/academic-drawing/blob/master/time-series/time_series_speed1.py) 81 | - [time_series_speed2.py](https://github.com/xinychen/academic-drawing/blob/master/time-series/time_series_speed2.py) 82 | 83 | and evaluate these in Python, then, you will see the following pictures: 84 | 85 | ![time_series_speed1](https://github.com/xinychen/academic-drawing/blob/master/time-series/time_series_speed1.png) 86 | 87 | ![time_series_speed2](https://github.com/xinychen/academic-drawing/blob/master/time-series/time_series_speed2.png) 88 | 89 | 90 | - **Download** 91 | - [road1_fiber_ms50_r10.csv](https://github.com/xinychen/academic-drawing/blob/master/time-series/road1_fiber_ms50_r10.csv) 92 | - [speed_curve.m](https://github.com/xinychen/academic-drawing/blob/master/time-series/speed_curve.m) 93 | 94 | and evaluate these in Matlab, then, you will see the following picture: 95 | 96 | ![speed_curve](https://github.com/xinychen/academic-drawing/blob/master/time-series/speed_curve.png) 97 | 98 | 99 | - **Download** 100 | - [nyc_data_completeness.mat](https://github.com/xinychen/academic-drawing/blob/master/time-series/nyc_data_completeness.mat) 101 | - [nyc_data_completeness.m](https://github.com/xinychen/academic-drawing/blob/master/time-series/nyc_data_completeness.m) 102 | 103 | and evaluate these in Matlab, then, you will see the following picture: 104 | 105 | ![nyc_data_completeness](https://github.com/xinychen/academic-drawing/blob/master/time-series/nyc_data_completeness.png) 106 | 107 | Our Publications 108 | -------------- 109 | 110 | Most of these examples are from our publications: 111 | 112 | - **Xinyu Chen**, Zhaocheng He, Yixian Chen, Yuhuan Lu, Jiawei Wang (2019). **Missing traffic data imputation and pattern discovery with a Bayesian augmented tensor factorization model**. Transportation Research Part C: Emerging Technologies, 104: 66-77. [[preprint](https://xinychen.github.io/paper/BATF.pdf)] [[doi](https://doi.org/10.1016/j.trc.2019.03.003)] [[slide](https://doi.org/10.5281/zenodo.2632552)] [[data](http://doi.org/10.5281/zenodo.1205229)] [[Matlab code](https://github.com/sysuits/BATF)] 113 | 114 | - **Xinyu Chen**, Zhaocheng He, Lijun Sun (2019). **A Bayesian tensor decomposition approach for spatiotemporal traffic data imputation**. Transportation Research Part C: Emerging Technologies, 98: 73-84. [[preprint](https://www.researchgate.net/publication/329177786_A_Bayesian_tensor_decomposition_approach_for_spatiotemporal_traffic_data_imputation)] [[doi](https://doi.org/10.1016/j.trc.2018.11.003)] [[data](http://doi.org/10.5281/zenodo.1205229)] [[Matlab code](https://github.com/lijunsun/bgcp_imputation)] [[Python code](https://github.com/xinychen/transdim/blob/master/experiments/Imputation-BGCP.ipynb)] 115 | 116 | > Please consider citing our papers if you find these codes help your research. 117 | -------------------------------------------------------------------------------- /curves/BCPF_fiber_rmselb_m30_r10.csv: -------------------------------------------------------------------------------- 1 | 12.107,5.0436,4.7282,4.6261,4.5801,4.5563,4.5417,4.5315,4.5236,4.5171,4.5115,4.5068,4.5027,4.4991,4.496,4.4933,4.4909,4.4888,4.4872,4.486,4.4854,4.4851,4.4851,4.4852,4.4852,4.4851,4.4849,4.4846,4.4842,4.4838,4.4834,4.483,4.4826,4.4823,4.482,4.4817,4.4815,4.4813,4.4811,4.481,4.4809,4.4809,4.4808,4.4808,4.4808,4.4809,4.4809,4.481,4.4811,4.4812,4.4813,4.4814,4.4815,4.4816,4.4818,4.4819,4.482,4.4821,4.4822,4.4824,4.4825,4.4826,4.4826,4.4827,4.4828,4.4829,4.4829,4.4829,4.483,4.483,4.483,4.4829,4.4829,4.4828,4.4828,4.4827,4.4826,4.4825,4.4823,4.4822,4.482,4.4819,4.4817,4.4815,4.4813,4.4811,4.4809,4.4807,4.4805,4.4804,4.4802,4.4801,4.48,4.48,4.48,4.48,4.4801,4.4802,4.4803,4.4805,4.4808,4.481,4.4813,4.4817,4.482,4.4824,4.4828,4.4833,4.4837,4.4842,4.4846,4.4851,4.4857,4.4862,4.4868,4.4874,4.4881,4.4888,4.4896,4.4904,4.4913,4.4922,4.4932,4.4942,4.4953,4.4964,4.4976,4.4988,4.5,4.5013,4.5026,4.5039,4.5052,4.5066,4.508,4.5094,4.5108,4.5122,4.5137,4.5151,4.5165,4.518,4.5194,4.5209,4.5223,4.5237,4.5252,4.5266,4.528,4.5294,4.5307,4.5321,4.5334,4.5348,4.5361,4.5374,4.5387,4.5399,4.5412,4.5424,4.5436,4.5448,4.546,4.5472,4.5483,4.5494,4.5505,4.5516,4.5527,4.5537,4.5547,4.5557,4.5567,4.5577,4.5586,4.5596,4.5605,4.5614,4.5622,4.5631,4.5639,4.5648,4.5656,4.5664,4.5672,4.5679,4.5687,4.5694,4.5701,4.5709,4.5716,4.5722,4.5729,4.5736,4.5742,4.5749,4.5755,4.5761,4.5768,4.5774 2 | 15.769,4.942,4.7001,4.6248,4.572,4.5498,4.5417,4.5383,4.5346,4.5302,4.5261,4.5231,4.5211,4.5202,4.52,4.5205,4.5214,4.5227,4.5242,4.526,4.5278,4.5298,4.5318,4.534,4.5362,4.5385,4.5409,4.5434,4.5461,4.5489,4.5518,4.555,4.5584,4.5621,4.5662,4.5706,4.5753,4.5806,4.5863,4.5926,4.5995,4.607,4.6152,4.6241,4.6339,4.6445,4.6559,4.6683,4.6817,4.6961,4.7115,4.728,4.7456,4.7643,4.784,4.8049,4.8269,4.85,4.8741,4.8992,4.9252,4.9522,4.98,5.0086,5.0379,5.0679,5.0984,5.1293,5.1607,5.1924,5.2243,5.2564,5.2886,5.3208,5.353,5.3852,5.4173,5.4492,5.481,5.5125,5.5438,5.5749,5.6058,5.6363,5.6666,5.6966,5.7263,5.7557,5.7848,5.8137,5.8422,5.8704,5.8982,5.9258,5.953,5.9799,6.0065,6.0327,6.0585,6.084,6.1091,6.1338,6.1581,6.1819,6.2054,6.2284,6.251,6.2731,6.2947,6.3158,6.3364,6.3565,6.376,6.3951,6.4135,6.4314,6.4488,6.4656,6.4818,6.4974,6.5124,6.5268,6.5407,6.554,6.5667,6.5788,6.5903,6.6012,6.6116,6.6214,6.6307,6.6394,6.6476,6.6553,6.6624,6.6691,6.6753,6.681,6.6862,6.691,6.6954,6.6994,6.703,6.7062,6.7091,6.7116,6.7139,6.7159,6.7176,6.7191,6.7204,6.7215,6.7225,6.7234,6.7243,6.7251,6.726,6.727,6.7281,6.7294,6.7309,6.7328,6.735,6.7377,6.741,6.7448,6.7493,6.7545,6.7606,6.7675,6.7755,6.7845,6.7946,6.8059,6.8184,6.8322,6.8472,6.8635,6.8811,6.8998,6.9196,6.9403,6.9617,6.9835,7.0054,7.027,7.0478,7.0673,7.0851,7.1007,7.1137,7.1235,7.1301,7.1331,7.1324,7.1282,7.1204,7.1093,7.095,7.0779 3 | 10.914,5.1576,4.7539,4.6333,4.5782,4.5539,4.5397,4.5301,4.5233,4.5182,4.5143,4.5112,4.5084,4.5057,4.5029,4.5005,4.4987,4.4971,4.4952,4.4932,4.4915,4.4901,4.4891,4.4884,4.4879,4.4875,4.4873,4.4871,4.487,4.487,4.4869,4.4869,4.4869,4.487,4.487,4.4871,4.4872,4.4873,4.4874,4.4876,4.4878,4.4881,4.4883,4.4887,4.489,4.4895,4.4899,4.4904,4.4909,4.4915,4.4922,4.4928,4.4935,4.4943,4.4951,4.4959,4.4968,4.4977,4.4986,4.4996,4.5005,4.5015,4.5025,4.5035,4.5046,4.5056,4.5066,4.5076,4.5086,4.5095,4.5104,4.5113,4.5122,4.513,4.5138,4.5146,4.5154,4.5161,4.5168,4.5175,4.5182,4.5189,4.5196,4.5202,4.5209,4.5216,4.5223,4.523,4.5238,4.5245,4.5252,4.526,4.5267,4.5274,4.5282,4.5289,4.5296,4.5303,4.531,4.5317,4.5323,4.533,4.5336,4.5341,4.5347,4.5352,4.5357,4.5362,4.5366,4.537,4.5374,4.5377,4.538,4.5382,4.5385,4.5387,4.5388,4.5389,4.539,4.5391,4.5392,4.5392,4.5391,4.5391,4.539,4.5389,4.5388,4.5387,4.5385,4.5383,4.5381,4.5379,4.5376,4.5374,4.5371,4.5368,4.5365,4.5361,4.5358,4.5354,4.535,4.5346,4.5342,4.5338,4.5334,4.5329,4.5325,4.532,4.5315,4.531,4.5305,4.53,4.5295,4.5289,4.5284,4.5279,4.5273,4.5267,4.5262,4.5256,4.525,4.5244,4.5238,4.5232,4.5226,4.522,4.5214,4.5208,4.5202,4.5196,4.519,4.5184,4.5178,4.5172,4.5166,4.516,4.5154,4.5148,4.5142,4.5136,4.513,4.5124,4.5118,4.5113,4.5107,4.5101,4.5096,4.509,4.5085,4.5079,4.5074,4.5069,4.5063,4.5058,4.5053,4.5048,4.5043,4.5038,4.5033,4.5028 4 | 15.341,5.2625,4.9899,4.7588,4.6361,4.5891,4.5661,4.5512,4.5401,4.5314,4.5246,4.5195,4.5157,4.5129,4.5111,4.5099,4.5093,4.5093,4.5097,4.5105,4.5118,4.5133,4.5151,4.5171,4.5192,4.5212,4.5231,4.5249,4.5263,4.5274,4.528,4.5283,4.5281,4.5276,4.5269,4.526,4.525,4.524,4.5231,4.5222,4.5213,4.5205,4.5198,4.5191,4.5184,4.5178,4.5171,4.5164,4.5158,4.5151,4.5144,4.5137,4.5131,4.5124,4.5117,4.511,4.5103,4.5096,4.509,4.5083,4.5077,4.507,4.5064,4.5059,4.5053,4.5048,4.5043,4.5038,4.5033,4.5029,4.5024,4.502,4.5017,4.5013,4.501,4.5007,4.5004,4.5001,4.4998,4.4996,4.4994,4.4992,4.499,4.4988,4.4986,4.4985,4.4984,4.4982,4.4981,4.498,4.4979,4.4978,4.4977,4.4976,4.4976,4.4975,4.4974,4.4974,4.4973,4.4973,4.4972,4.4972,4.4972,4.4971,4.4971,4.4971,4.497,4.497,4.497,4.497,4.497,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.4969,4.497,4.497,4.497,4.497,4.497,4.497,4.497,4.497,4.497,4.497,4.497,4.4971,4.4971,4.4971,4.4971,4.497,4.497,4.497,4.497,4.497,4.497,4.4969,4.4969,4.4968,4.4968,4.4967,4.4966,4.4965,4.4964,4.4963,4.4962,4.496,4.4959,4.4957,4.4955,4.4953,4.4951,4.4948,4.4946,4.4943,4.494,4.4937,4.4934,4.4931,4.4928,4.4925,4.4922,4.4919,4.4915,4.4912,4.491,4.4907,4.4904,4.4901,4.4899,4.4896,4.4894,4.4891,4.4889,4.4887,4.4885,4.4884,4.4882,4.488,4.4879 5 | 13.858,4.8746,4.6488,4.5712,4.5364,4.5169,4.5041,4.4946,4.4872,4.4811,4.4759,4.4716,4.4678,4.4647,4.462,4.4597,4.4578,4.4561,4.4549,4.4542,4.4544,4.4559,4.4593,4.4646,4.4715,4.4794,4.487,4.4932,4.4971,4.4989,4.499,4.498,4.4966,4.4952,4.4939,4.4928,4.4919,4.4913,4.491,4.491,4.4911,4.4914,4.4919,4.4924,4.493,4.4937,4.4945,4.4953,4.4961,4.497,4.498,4.499,4.5,4.5011,4.5022,4.5034,4.5046,4.5058,4.507,4.5082,4.5095,4.5107,4.5119,4.5132,4.5144,4.5156,4.5168,4.518,4.5192,4.5204,4.5215,4.5227,4.5239,4.525,4.5262,4.5273,4.5285,4.5296,4.5308,4.532,4.5331,4.5343,4.5354,4.5366,4.5378,4.5389,4.5401,4.5413,4.5424,4.5436,4.5447,4.5458,4.547,4.5481,4.5491,4.5502,4.5513,4.5523,4.5532,4.5542,4.5551,4.556,4.5568,4.5576,4.5583,4.559,4.5596,4.5603,4.5608,4.5614,4.5619,4.5624,4.5628,4.5633,4.5638,4.5643,4.5649,4.5655,4.5661,4.5668,4.5676,4.5684,4.5694,4.5703,4.5713,4.5724,4.5735,4.5747,4.5759,4.5771,4.5784,4.5796,4.5809,4.5822,4.5836,4.5849,4.5862,4.5876,4.589,4.5904,4.5918,4.5933,4.5947,4.5963,4.5978,4.5994,4.6011,4.6028,4.6045,4.6064,4.6083,4.6102,4.6123,4.6144,4.6166,4.6189,4.6213,4.6238,4.6264,4.6291,4.632,4.6349,4.638,4.6412,4.6445,4.648,4.6516,4.6553,4.6592,4.6632,4.6674,4.6716,4.6761,4.6806,4.6853,4.6901,4.6951,4.7002,4.7054,4.7107,4.7161,4.7217,4.7274,4.7332,4.7391,4.745,4.7511,4.7573,4.7636,4.77,4.7764,4.783,4.7896,4.7963,4.8031,4.81,4.8169,4.8239,4.8309,4.838 6 | 14.209,5.2518,5.0111,4.8447,4.7513,4.7128,4.6887,4.6716,4.6586,4.648,4.6392,4.6313,4.6239,4.616,4.6069,4.5954,4.5809,4.5641,4.5475,4.5333,4.5219,4.5128,4.505,4.4982,4.4921,4.4865,4.4816,4.4772,4.4734,4.4701,4.4674,4.465,4.4629,4.4612,4.4597,4.4584,4.4572,4.4561,4.4551,4.4542,4.4534,4.4526,4.4518,4.4511,4.4504,4.4498,4.4492,4.4487,4.4482,4.4478,4.4474,4.447,4.4467,4.4465,4.4463,4.4461,4.446,4.4459,4.4458,4.4458,4.4458,4.4459,4.446,4.446,4.4461,4.4462,4.4463,4.4464,4.4465,4.4466,4.4467,4.4467,4.4468,4.4468,4.4468,4.4468,4.4468,4.4467,4.4467,4.4466,4.4466,4.4465,4.4464,4.4462,4.4461,4.446,4.4458,4.4457,4.4455,4.4454,4.4452,4.445,4.4448,4.4447,4.4445,4.4443,4.4441,4.4439,4.4438,4.4436,4.4434,4.4433,4.4431,4.4429,4.4428,4.4427,4.4425,4.4424,4.4423,4.4421,4.442,4.4419,4.4418,4.4417,4.4416,4.4416,4.4415,4.4414,4.4413,4.4413,4.4412,4.4412,4.4412,4.4411,4.4411,4.4411,4.441,4.441,4.441,4.441,4.441,4.441,4.441,4.4411,4.4411,4.4411,4.4411,4.4412,4.4412,4.4413,4.4413,4.4414,4.4414,4.4415,4.4416,4.4416,4.4417,4.4418,4.4419,4.442,4.4421,4.4422,4.4423,4.4424,4.4425,4.4426,4.4427,4.4428,4.443,4.4431,4.4432,4.4434,4.4435,4.4437,4.4438,4.444,4.4441,4.4443,4.4445,4.4446,4.4448,4.445,4.4452,4.4453,4.4455,4.4457,4.4459,4.4461,4.4463,4.4465,4.4467,4.4469,4.4471,4.4474,4.4476,4.4478,4.448,4.4483,4.4485,4.4487,4.4489,4.4492,4.4494,4.4497,4.4499,4.4502,4.4504,4.4507,4.4509,4.4512 7 | 15.927,5.1085,4.6811,4.5991,4.5575,4.5319,4.5167,4.5068,4.4996,4.4938,4.4886,4.484,4.4798,4.4763,4.4733,4.4707,4.468,4.4652,4.4625,4.4599,4.4576,4.4555,4.4536,4.452,4.4506,4.4494,4.4485,4.4478,4.4473,4.447,4.4468,4.4468,4.4469,4.447,4.4473,4.4476,4.448,4.4484,4.4488,4.4493,4.4497,4.4502,4.4507,4.4512,4.4517,4.4522,4.4526,4.4531,4.4536,4.4541,4.4545,4.455,4.4555,4.4559,4.4564,4.4569,4.4574,4.4579,4.4585,4.459,4.4595,4.4601,4.4607,4.4613,4.4619,4.4625,4.4631,4.4637,4.4644,4.4651,4.4657,4.4664,4.4671,4.4678,4.4685,4.4692,4.4699,4.4706,4.4713,4.472,4.4727,4.4734,4.4741,4.4748,4.4755,4.4762,4.4769,4.4776,4.4783,4.4791,4.4798,4.4805,4.4812,4.4819,4.4827,4.4834,4.4841,4.4849,4.4856,4.4863,4.4871,4.4878,4.4885,4.4893,4.49,4.4907,4.4915,4.4922,4.493,4.4937,4.4944,4.4951,4.4959,4.4966,4.4974,4.4981,4.4989,4.4996,4.5004,4.5012,4.5021,4.5029,4.5038,4.5048,4.5058,4.5068,4.5079,4.5091,4.5103,4.5116,4.5129,4.5144,4.5159,4.5175,4.5192,4.5211,4.523,4.525,4.5271,4.5293,4.5316,4.5341,4.5366,4.5392,4.542,4.5448,4.5478,4.5508,4.5539,4.5571,4.5604,4.5637,4.5672,4.5707,4.5742,4.5778,4.5815,4.5852,4.5889,4.5927,4.5965,4.6003,4.6042,4.608,4.6119,4.6157,4.6196,4.6234,4.6273,4.6311,4.6349,4.6387,4.6425,4.6463,4.6501,4.6538,4.6575,4.6612,4.6649,4.6686,4.6722,4.6758,4.6794,4.6829,4.6864,4.6899,4.6934,4.6968,4.7002,4.7036,4.707,4.7104,4.7137,4.717,4.7202,4.7235,4.7267,4.7299,4.7331,4.7363 8 | 10.21,4.9035,4.7275,4.6891,4.6629,4.6433,4.6276,4.6153,4.6061,4.5992,4.5939,4.5895,4.5857,4.5822,4.5788,4.5753,4.5716,4.5678,4.5638,4.5598,4.5558,4.5518,4.5481,4.5445,4.5412,4.5383,4.5357,4.5334,4.5315,4.53,4.5287,4.5278,4.5271,4.5265,4.5262,4.5259,4.5257,4.5255,4.5252,4.5249,4.5245,4.524,4.5233,4.5225,4.5216,4.5206,4.5194,4.5182,4.5169,4.5156,4.5143,4.5129,4.5115,4.5102,4.5089,4.5076,4.5064,4.5052,4.504,4.5029,4.5018,4.5008,4.4998,4.4989,4.498,4.4971,4.4963,4.4955,4.4948,4.494,4.4934,4.4927,4.4921,4.4915,4.4909,4.4903,4.4898,4.4893,4.4888,4.4883,4.4878,4.4873,4.4868,4.4864,4.4859,4.4854,4.485,4.4845,4.484,4.4835,4.483,4.4825,4.4819,4.4814,4.4808,4.4803,4.4797,4.4791,4.4785,4.4778,4.4772,4.4766,4.4759,4.4753,4.4747,4.4741,4.4735,4.473,4.4725,4.472,4.4715,4.4711,4.4707,4.4704,4.47,4.4697,4.4695,4.4692,4.469,4.4688,4.4686,4.4684,4.4682,4.4681,4.4679,4.4678,4.4676,4.4675,4.4673,4.4672,4.4671,4.467,4.4668,4.4667,4.4666,4.4665,4.4664,4.4663,4.4661,4.466,4.4659,4.4658,4.4657,4.4656,4.4655,4.4653,4.4652,4.4651,4.465,4.4648,4.4647,4.4646,4.4644,4.4643,4.4642,4.464,4.4638,4.4637,4.4635,4.4633,4.4632,4.463,4.4628,4.4626,4.4624,4.4622,4.462,4.4617,4.4615,4.4612,4.461,4.4607,4.4605,4.4602,4.4599,4.4596,4.4593,4.459,4.4587,4.4583,4.458,4.4576,4.4572,4.4569,4.4565,4.4561,4.4557,4.4553,4.4549,4.4545,4.454,4.4536,4.4532,4.4528,4.4524,4.4519,4.4515,4.4511,4.4507,4.4503 9 | 11.63,4.9045,4.6813,4.5955,4.5558,4.5362,4.5237,4.5138,4.5047,4.4958,4.4872,4.4794,4.4728,4.4675,4.4634,4.46,4.4568,4.4539,4.4513,4.449,4.447,4.4454,4.4441,4.443,4.442,4.4412,4.4405,4.4399,4.4394,4.439,4.4386,4.4383,4.4381,4.4379,4.4377,4.4376,4.4376,4.4375,4.4376,4.4376,4.4377,4.4379,4.4381,4.4383,4.4387,4.439,4.4395,4.44,4.4407,4.4414,4.4421,4.443,4.4439,4.4449,4.446,4.4471,4.4482,4.4494,4.4507,4.4524,4.455,4.4589,4.4651,4.4734,4.4828,4.4922,4.5016,4.511,4.5209,4.5314,4.5425,4.5543,4.5669,4.5801,4.594,4.6087,4.6241,4.6403,4.6574,4.6754,4.6944,4.7144,4.7355,4.7578,4.7813,4.8061,4.8322,4.8597,4.8886,4.9189,4.9507,4.9839,5.0186,5.0547,5.0923,5.1312,5.1715,5.2131,5.2559,5.2998,5.3448,5.3909,5.4378,5.4856,5.5341,5.5832,5.6329,5.683,5.7335,5.7842,5.8351,5.8861,5.937,5.9879,6.0386,6.089,6.1391,6.1887,6.2378,6.2864,6.3343,6.3815,6.4279,6.4735,6.5181,6.5618,6.6043,6.6458,6.6861,6.7252,6.763,6.7995,6.8346,6.8683,6.9005,6.9312,6.9605,6.9882,7.0145,7.0392,7.0625,7.0844,7.1048,7.124,7.1418,7.1586,7.1742,7.1889,7.2028,7.216,7.2286,7.2408,7.2527,7.2644,7.2762,7.2881,7.3004,7.3131,7.3264,7.3405,7.3554,7.3712,7.3881,7.4061,7.4253,7.4457,7.4674,7.4905,7.5148,7.5405,7.5675,7.5957,7.6252,7.6558,7.6875,7.7202,7.754,7.7886,7.824,7.8602,7.8971,7.9346,7.9726,8.0111,8.0501,8.0894,8.129,8.169,8.2092,8.2496,8.2902,8.3309,8.3718,8.4128,8.4539,8.4952,8.5365,8.5779,8.6193,8.6609 10 | 15.557,4.9395,4.6386,4.565,4.528,4.506,4.4925,4.4843,4.4807,4.4823,4.4886,4.4965,4.5039,4.5097,4.5139,4.5168,4.5186,4.5196,4.52,4.5198,4.5193,4.5186,4.5178,4.5169,4.5161,4.5154,4.515,4.5148,4.515,4.5154,4.5163,4.5175,4.5192,4.5213,4.5239,4.527,4.5306,4.5348,4.5396,4.5449,4.5509,4.5575,4.5648,4.5728,4.5814,4.5908,4.6008,4.6116,4.623,4.6351,4.6479,4.6613,4.6752,4.6897,4.7047,4.7202,4.736,4.7522,4.7687,4.7854,4.8022,4.8192,4.8362,4.8532,4.8701,4.887,4.9037,4.9203,4.9367,4.9529,4.9688,4.9846,5.0001,5.0153,5.0303,5.0451,5.0597,5.0741,5.0883,5.1023,5.1161,5.1299,5.1435,5.1569,5.1704,5.1837,5.197,5.2102,5.2234,5.2367,5.2499,5.2631,5.2763,5.2896,5.3029,5.3162,5.3296,5.343,5.3564,5.3698,5.3833,5.3968,5.4103,5.4239,5.4374,5.4509,5.4644,5.4779,5.4913,5.5047,5.518,5.5313,5.5444,5.5575,5.5705,5.5833,5.5961,5.6087,5.6211,5.6335,5.6456,5.6576,5.6695,5.6811,5.6926,5.7039,5.7151,5.726,5.7368,5.7474,5.7578,5.7681,5.7782,5.7881,5.7979,5.8076,5.8171,5.8265,5.8359,5.8451,5.8542,5.8633,5.8724,5.8814,5.8905,5.8995,5.9086,5.9177,5.9269,5.9362,5.9457,5.9553,5.9651,5.975,5.9852,5.9957,6.0064,6.0175,6.0289,6.0407,6.0529,6.0655,6.0786,6.0922,6.1064,6.1211,6.1365,6.1525,6.1692,6.1867,6.2049,6.2239,6.2438,6.2646,6.2863,6.309,6.3328,6.3576,6.3836,6.4108,6.4391,6.4688,6.4998,6.5322,6.5659,6.6012,6.638,6.6763,6.7163,6.7579,6.8011,6.8462,6.8929,6.9415,6.9919,7.0442,7.0982,7.1542,7.212,7.2717 11 | -1.3276e+07,-1.2108e+07,-1.2024e+07,-1.1999e+07,-1.199e+07,-1.1985e+07,-1.1981e+07,-1.1978e+07,-1.1976e+07,-1.1974e+07,-1.1973e+07,-1.1972e+07,-1.1971e+07,-1.197e+07,-1.1969e+07,-1.1968e+07,-1.1967e+07,-1.1966e+07,-1.1966e+07,-1.1965e+07,-1.1965e+07,-1.1964e+07,-1.1963e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.1961e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07 12 | -1.3617e+07,-1.2111e+07,-1.2021e+07,-1.1997e+07,-1.1985e+07,-1.1979e+07,-1.1975e+07,-1.1972e+07,-1.197e+07,-1.1968e+07,-1.1967e+07,-1.1966e+07,-1.1965e+07,-1.1964e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.1961e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07 13 | -1.3138e+07,-1.2121e+07,-1.2035e+07,-1.2007e+07,-1.1994e+07,-1.1987e+07,-1.1983e+07,-1.198e+07,-1.1977e+07,-1.1976e+07,-1.1974e+07,-1.1973e+07,-1.1972e+07,-1.197e+07,-1.1969e+07,-1.1968e+07,-1.1967e+07,-1.1965e+07,-1.1964e+07,-1.1964e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07 14 | -1.3547e+07,-1.2159e+07,-1.2079e+07,-1.203e+07,-1.2002e+07,-1.1989e+07,-1.1982e+07,-1.1977e+07,-1.1974e+07,-1.1972e+07,-1.197e+07,-1.1968e+07,-1.1967e+07,-1.1966e+07,-1.1965e+07,-1.1965e+07,-1.1964e+07,-1.1964e+07,-1.1963e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.1961e+07,-1.1961e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07 15 | -1.3427e+07,-1.208e+07,-1.2006e+07,-1.1988e+07,-1.1981e+07,-1.1977e+07,-1.1974e+07,-1.1972e+07,-1.197e+07,-1.1968e+07,-1.1967e+07,-1.1966e+07,-1.1965e+07,-1.1964e+07,-1.1964e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07 16 | -1.3446e+07,-1.2175e+07,-1.2101e+07,-1.206e+07,-1.204e+07,-1.2031e+07,-1.2025e+07,-1.2021e+07,-1.2017e+07,-1.2014e+07,-1.2011e+07,-1.2008e+07,-1.2006e+07,-1.2003e+07,-1.1999e+07,-1.1995e+07,-1.1991e+07,-1.1987e+07,-1.1983e+07,-1.198e+07,-1.1977e+07,-1.1974e+07,-1.1972e+07,-1.197e+07,-1.1969e+07,-1.1967e+07,-1.1966e+07,-1.1964e+07,-1.1963e+07,-1.1962e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07 17 | -1.3624e+07,-1.2111e+07,-1.2014e+07,-1.1995e+07,-1.1985e+07,-1.1978e+07,-1.1973e+07,-1.197e+07,-1.1967e+07,-1.1965e+07,-1.1964e+07,-1.1962e+07,-1.1961e+07,-1.196e+07,-1.1959e+07,-1.1958e+07,-1.1957e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07,-1.1935e+07 18 | -1.3039e+07,-1.2084e+07,-1.2031e+07,-1.2018e+07,-1.2009e+07,-1.2002e+07,-1.1997e+07,-1.1992e+07,-1.1989e+07,-1.1986e+07,-1.1983e+07,-1.1981e+07,-1.1978e+07,-1.1977e+07,-1.1975e+07,-1.1973e+07,-1.1972e+07,-1.1971e+07,-1.197e+07,-1.1969e+07,-1.1968e+07,-1.1967e+07,-1.1967e+07,-1.1966e+07,-1.1965e+07,-1.1965e+07,-1.1964e+07,-1.1964e+07,-1.1963e+07,-1.1963e+07,-1.1963e+07,-1.1962e+07,-1.1962e+07,-1.1961e+07,-1.1961e+07,-1.1961e+07,-1.196e+07,-1.196e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07 19 | -1.3204e+07,-1.2079e+07,-1.2006e+07,-1.1991e+07,-1.1985e+07,-1.1981e+07,-1.1978e+07,-1.1976e+07,-1.1974e+07,-1.1972e+07,-1.197e+07,-1.1968e+07,-1.1966e+07,-1.1965e+07,-1.1963e+07,-1.1962e+07,-1.1961e+07,-1.196e+07,-1.1959e+07,-1.1959e+07,-1.1958e+07,-1.1958e+07,-1.1957e+07,-1.1957e+07,-1.1956e+07,-1.1956e+07,-1.1956e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07 20 | -1.3563e+07,-1.2099e+07,-1.2002e+07,-1.1984e+07,-1.1974e+07,-1.1969e+07,-1.1966e+07,-1.1963e+07,-1.1961e+07,-1.196e+07,-1.1958e+07,-1.1957e+07,-1.1956e+07,-1.1955e+07,-1.1954e+07,-1.1954e+07,-1.1953e+07,-1.1952e+07,-1.1952e+07,-1.1951e+07,-1.195e+07,-1.195e+07,-1.195e+07,-1.1949e+07,-1.1949e+07,-1.1948e+07,-1.1948e+07,-1.1948e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1947e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1946e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1945e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1944e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1943e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1942e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.1941e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.194e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1939e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1938e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1937e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07,-1.1936e+07 21 | -------------------------------------------------------------------------------- /curves/BCPF_fiber_rmselb_m30_r5.csv: -------------------------------------------------------------------------------- 1 | 19.58,5.5326,4.9695,4.9248,4.9128,4.9052,4.8979,4.89,4.8817,4.8734,4.866,4.8599,4.8551,4.8514,4.8485,4.8462,4.8444,4.8429,4.8416,4.8404,4.8394,4.8384,4.8375,4.8367,4.836,4.8352,4.8345,4.8339,4.8332,4.8326,4.832,4.8315,4.8309,4.8304,4.8299,4.8294,4.8289,4.8284,4.8279,4.8275,4.827,4.8266,4.8261,4.8257,4.8253,4.8249,4.8245,4.8241,4.8237,4.8233,4.8229,4.8226,4.8222,4.8218,4.8215,4.8211,4.8208,4.8205,4.8201,4.8198,4.8195,4.8192,4.8189,4.8186,4.8183,4.818,4.8177,4.8175,4.8172,4.8169,4.8167,4.8165,4.8162,4.816,4.8158,4.8155,4.8153,4.8151,4.8149,4.8147,4.8145,4.8143,4.8142,4.814,4.8138,4.8136,4.8135,4.8133,4.8132,4.813,4.8129,4.8127,4.8126,4.8124,4.8123,4.8121,4.812,4.8119,4.8118,4.8116,4.8115,4.8114,4.8113,4.8111,4.811,4.8109,4.8108,4.8107,4.8106,4.8105,4.8104,4.8102,4.8101,4.81,4.8099,4.8098,4.8097,4.8096,4.8095,4.8094,4.8093,4.8092,4.8091,4.809,4.8089,4.8088,4.8087,4.8086,4.8085,4.8084,4.8083,4.8082,4.8081,4.808,4.8079,4.8078,4.8077,4.8076,4.8075,4.8074,4.8073,4.8072,4.8071,4.807,4.8069,4.8068,4.8067,4.8065,4.8064,4.8063,4.8062,4.8061,4.806,4.8059,4.8057,4.8056,4.8055,4.8053,4.8052,4.8051,4.8049,4.8048,4.8047,4.8045,4.8043,4.8042,4.804,4.8038,4.8037,4.8035,4.8033,4.8031,4.8029,4.8027,4.8024,4.8022,4.8019,4.8016,4.8013,4.801,4.8006,4.8002,4.7998,4.7993,4.7988,4.7982,4.7975,4.7967,4.7957,4.7945,4.793,4.7911,4.7887,4.7854,4.7809,4.7749,4.7674,4.7594,4.7528,4.7486 2 | 14.368,5.459,5.2584,5.2604,5.2499,5.1767,5.0744,5.0173,4.9907,4.9737,4.9611,4.9515,4.9443,4.9389,4.9347,4.9314,4.9288,4.9265,4.9246,4.9228,4.9211,4.9194,4.9177,4.9159,4.9139,4.9115,4.9088,4.9055,4.9016,4.8972,4.8923,4.8874,4.8828,4.8788,4.8754,4.8727,4.8706,4.8688,4.8673,4.8661,4.865,4.8641,4.8632,4.8625,4.8618,4.8611,4.8605,4.8599,4.8593,4.8588,4.8583,4.8578,4.8573,4.8569,4.8565,4.856,4.8556,4.8552,4.8548,4.8544,4.8541,4.8537,4.8533,4.853,4.8526,4.8523,4.852,4.8516,4.8513,4.851,4.8507,4.8504,4.8501,4.8498,4.8495,4.8492,4.8489,4.8487,4.8484,4.8481,4.8478,4.8476,4.8473,4.8471,4.8468,4.8466,4.8463,4.8461,4.8459,4.8456,4.8454,4.8452,4.8449,4.8447,4.8445,4.8443,4.844,4.8438,4.8436,4.8434,4.8432,4.843,4.8428,4.8426,4.8424,4.8422,4.842,4.8418,4.8416,4.8414,4.8413,4.8411,4.8409,4.8407,4.8405,4.8403,4.8402,4.84,4.8398,4.8396,4.8395,4.8393,4.8391,4.839,4.8388,4.8386,4.8385,4.8383,4.8382,4.838,4.8378,4.8377,4.8375,4.8374,4.8372,4.8371,4.8369,4.8367,4.8366,4.8364,4.8363,4.8361,4.836,4.8359,4.8357,4.8356,4.8354,4.8353,4.8351,4.835,4.8348,4.8347,4.8346,4.8344,4.8343,4.8341,4.834,4.8339,4.8337,4.8336,4.8335,4.8333,4.8332,4.8331,4.8329,4.8328,4.8327,4.8326,4.8324,4.8323,4.8322,4.8321,4.8319,4.8318,4.8317,4.8316,4.8315,4.8313,4.8312,4.8311,4.831,4.8309,4.8308,4.8306,4.8305,4.8304,4.8303,4.8302,4.8301,4.83,4.8299,4.8298,4.8297,4.8296,4.8295,4.8294,4.8293,4.8292,4.8291,4.829 3 | 14.027,5.5108,5.2105,5.1653,5.1511,5.1408,5.1256,5.1047,5.0792,5.0512,5.0226,4.9943,4.9659,4.9373,4.91,4.8863,4.8675,4.8535,4.8435,4.8364,4.8314,4.8278,4.825,4.8227,4.8208,4.8191,4.8176,4.8161,4.8146,4.8132,4.8118,4.8104,4.809,4.8076,4.8063,4.805,4.8038,4.8027,4.8016,4.8007,4.7998,4.799,4.7983,4.7977,4.7971,4.7966,4.7961,4.7957,4.7954,4.7951,4.7948,4.7945,4.7943,4.7941,4.7939,4.7937,4.7935,4.7933,4.7931,4.7928,4.7926,4.7923,4.792,4.7917,4.7913,4.791,4.7905,4.7901,4.7896,4.789,4.7884,4.7878,4.7872,4.7865,4.7858,4.7851,4.7844,4.7836,4.7829,4.7821,4.7814,4.7806,4.7798,4.7791,4.7783,4.7776,4.7769,4.7761,4.7754,4.7747,4.7741,4.7734,4.7727,4.7721,4.7715,4.7709,4.7703,4.7697,4.7691,4.7685,4.768,4.7675,4.7669,4.7664,4.7659,4.7654,4.765,4.7645,4.764,4.7636,4.7631,4.7627,4.7623,4.7618,4.7614,4.761,4.7606,4.7602,4.7598,4.7594,4.759,4.7586,4.7582,4.7578,4.7574,4.757,4.7566,4.7562,4.7558,4.7554,4.755,4.7546,4.7542,4.7538,4.7534,4.753,4.7525,4.7521,4.7516,4.7512,4.7507,4.7502,4.7497,4.7491,4.7486,4.748,4.7474,4.7467,4.746,4.7453,4.7445,4.7437,4.7428,4.7419,4.7409,4.7399,4.7388,4.7377,4.7365,4.7354,4.7343,4.7332,4.7321,4.7311,4.7302,4.7293,4.7284,4.7277,4.727,4.7263,4.7257,4.7252,4.7246,4.7241,4.7237,4.7232,4.7228,4.7224,4.722,4.7216,4.7212,4.7209,4.7205,4.7202,4.7198,4.7195,4.7191,4.7188,4.7185,4.7181,4.7178,4.7175,4.7172,4.7168,4.7165,4.7162,4.7159,4.7155,4.7152,4.7149 4 | 19.071,5.5489,5.1207,5.0667,5.0226,4.9761,4.9332,4.8983,4.8722,4.8533,4.8397,4.8296,4.8218,4.8154,4.8099,4.8047,4.7996,4.7942,4.7884,4.7822,4.7755,4.7685,4.7617,4.7551,4.749,4.7435,4.7387,4.7345,4.7308,4.7277,4.725,4.7227,4.7206,4.7189,4.7173,4.7159,4.7147,4.7136,4.7126,4.7116,4.7108,4.71,4.7093,4.7086,4.7079,4.7073,4.7068,4.7062,4.7057,4.7052,4.7047,4.7043,4.7038,4.7034,4.703,4.7026,4.7022,4.7019,4.7015,4.7012,4.7009,4.7005,4.7002,4.6999,4.6996,4.6993,4.6991,4.6988,4.6985,4.6983,4.698,4.6978,4.6975,4.6973,4.6971,4.6969,4.6967,4.6965,4.6963,4.6961,4.6959,4.6957,4.6955,4.6953,4.6952,4.695,4.6948,4.6947,4.6945,4.6944,4.6942,4.6941,4.694,4.6938,4.6937,4.6936,4.6934,4.6933,4.6932,4.6931,4.693,4.6929,4.6928,4.6926,4.6925,4.6924,4.6924,4.6923,4.6922,4.6921,4.692,4.6919,4.6918,4.6917,4.6917,4.6916,4.6915,4.6914,4.6914,4.6913,4.6912,4.6912,4.6911,4.691,4.691,4.6909,4.6909,4.6908,4.6907,4.6907,4.6906,4.6906,4.6905,4.6905,4.6904,4.6904,4.6903,4.6903,4.6903,4.6902,4.6902,4.6901,4.6901,4.6901,4.69,4.69,4.6899,4.6899,4.6899,4.6898,4.6898,4.6898,4.6897,4.6897,4.6897,4.6896,4.6896,4.6896,4.6896,4.6895,4.6895,4.6895,4.6894,4.6894,4.6894,4.6894,4.6893,4.6893,4.6893,4.6893,4.6892,4.6892,4.6892,4.6892,4.6892,4.6891,4.6891,4.6891,4.6891,4.6891,4.689,4.689,4.689,4.689,4.689,4.6889,4.6889,4.6889,4.6889,4.6889,4.6889,4.6888,4.6888,4.6888,4.6888,4.6888,4.6888,4.6887,4.6887,4.6887 5 | 13.585,5.4415,5.128,5.0059,4.9294,4.8682,4.8228,4.7924,4.7727,4.7593,4.7498,4.7426,4.737,4.7326,4.729,4.726,4.7235,4.7214,4.7195,4.7179,4.7165,4.7152,4.7141,4.7131,4.7122,4.7113,4.7106,4.7098,4.7092,4.7085,4.7079,4.7074,4.7068,4.7063,4.7059,4.7054,4.705,4.7046,4.7042,4.7038,4.7034,4.7031,4.7027,4.7024,4.7021,4.7018,4.7015,4.7012,4.701,4.7007,4.7005,4.7002,4.7,4.6998,4.6996,4.6993,4.6992,4.699,4.6988,4.6986,4.6984,4.6983,4.6981,4.698,4.6978,4.6977,4.6975,4.6974,4.6973,4.6972,4.697,4.6969,4.6968,4.6967,4.6966,4.6965,4.6964,4.6963,4.6963,4.6962,4.6961,4.696,4.6959,4.6959,4.6958,4.6957,4.6957,4.6956,4.6955,4.6955,4.6954,4.6954,4.6953,4.6953,4.6952,4.6952,4.6951,4.6951,4.695,4.695,4.6949,4.6949,4.6948,4.6948,4.6948,4.6947,4.6947,4.6947,4.6946,4.6946,4.6946,4.6945,4.6945,4.6945,4.6944,4.6944,4.6944,4.6944,4.6943,4.6943,4.6943,4.6943,4.6942,4.6942,4.6942,4.6942,4.6942,4.6941,4.6941,4.6941,4.6941,4.6941,4.694,4.694,4.694,4.694,4.694,4.6939,4.6939,4.6939,4.6939,4.6939,4.6939,4.6938,4.6938,4.6938,4.6938,4.6938,4.6938,4.6938,4.6937,4.6937,4.6937,4.6937,4.6937,4.6937,4.6937,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932 6 | 13.934,5.1914,4.9623,4.8941,4.8508,4.8158,4.7878,4.7666,4.7512,4.7403,4.7325,4.7269,4.7227,4.7196,4.7171,4.7152,4.7137,4.7125,4.7115,4.7106,4.7099,4.7092,4.7087,4.7082,4.7078,4.7074,4.707,4.7067,4.7064,4.7061,4.7058,4.7055,4.7053,4.705,4.7048,4.7046,4.7044,4.7042,4.704,4.7038,4.7036,4.7034,4.7032,4.7031,4.7029,4.7027,4.7026,4.7024,4.7022,4.7021,4.7019,4.7018,4.7016,4.7015,4.7014,4.7012,4.7011,4.7009,4.7008,4.7007,4.7005,4.7004,4.7003,4.7001,4.7,4.6999,4.6998,4.6996,4.6995,4.6994,4.6993,4.6991,4.699,4.6989,4.6988,4.6986,4.6985,4.6984,4.6983,4.6982,4.698,4.6979,4.6978,4.6977,4.6976,4.6974,4.6973,4.6972,4.697,4.6969,4.6968,4.6967,4.6965,4.6964,4.6963,4.6961,4.696,4.6958,4.6957,4.6955,4.6954,4.6952,4.6951,4.6949,4.6947,4.6945,4.6944,4.6942,4.694,4.6938,4.6936,4.6934,4.6931,4.6929,4.6927,4.6924,4.6921,4.6918,4.6915,4.6912,4.6909,4.6905,4.6902,4.6897,4.6893,4.6889,4.6884,4.6878,4.6873,4.6867,4.686,4.6853,4.6846,4.6838,4.6829,4.682,4.6811,4.6801,4.6791,4.678,4.677,4.6759,4.6749,4.6739,4.673,4.6721,4.6712,4.6704,4.6697,4.6691,4.6685,4.668,4.6675,4.6671,4.6667,4.6663,4.666,4.6657,4.6655,4.6652,4.665,4.6648,4.6646,4.6645,4.6643,4.6642,4.664,4.6639,4.6638,4.6637,4.6636,4.6635,4.6634,4.6633,4.6632,4.6631,4.6631,4.663,4.6629,4.6628,4.6628,4.6627,4.6626,4.6626,4.6625,4.6625,4.6624,4.6624,4.6623,4.6622,4.6622,4.6621,4.6621,4.662,4.662,4.6619,4.6619,4.6618,4.6618,4.6617 7 | 16.956,5.2727,5.0142,4.9808,4.968,4.9584,4.949,4.9379,4.9234,4.9045,4.883,4.8608,4.8381,4.8146,4.7916,4.7712,4.7549,4.7427,4.7336,4.727,4.7222,4.7188,4.7164,4.7147,4.7134,4.7124,4.7117,4.711,4.7105,4.71,4.7096,4.7092,4.7088,4.7085,4.7082,4.7079,4.7076,4.7073,4.7071,4.7068,4.7065,4.7063,4.7061,4.7058,4.7056,4.7054,4.7052,4.705,4.7048,4.7046,4.7044,4.7042,4.704,4.7038,4.7036,4.7034,4.7032,4.703,4.7029,4.7027,4.7025,4.7023,4.7022,4.702,4.7018,4.7017,4.7015,4.7014,4.7012,4.7011,4.7009,4.7008,4.7006,4.7005,4.7003,4.7002,4.7001,4.6999,4.6998,4.6997,4.6996,4.6995,4.6993,4.6992,4.6991,4.699,4.6989,4.6988,4.6987,4.6986,4.6985,4.6984,4.6984,4.6983,4.6982,4.6981,4.698,4.698,4.6979,4.6978,4.6978,4.6977,4.6976,4.6976,4.6975,4.6975,4.6974,4.6973,4.6973,4.6972,4.6972,4.6971,4.6971,4.697,4.697,4.697,4.6969,4.6969,4.6968,4.6968,4.6967,4.6967,4.6967,4.6966,4.6966,4.6965,4.6965,4.6965,4.6964,4.6964,4.6964,4.6963,4.6963,4.6963,4.6963,4.6962,4.6962,4.6962,4.6961,4.6961,4.6961,4.696,4.696,4.696,4.696,4.6959,4.6959,4.6959,4.6959,4.6958,4.6958,4.6958,4.6958,4.6958,4.6957,4.6957,4.6957,4.6957,4.6956,4.6956,4.6956,4.6956,4.6955,4.6955,4.6955,4.6955,4.6955,4.6954,4.6954,4.6954,4.6954,4.6954,4.6953,4.6953,4.6953,4.6953,4.6952,4.6952,4.6952,4.6952,4.6951,4.6951,4.6951,4.6951,4.695,4.695,4.695,4.6949,4.6949,4.6949,4.6948,4.6948,4.6948,4.6947,4.6947,4.6946,4.6946,4.6945,4.6945,4.6944 8 | 17.339,5.6266,5.5091,5.4607,5.4452,5.4631,5.3519,5.0423,4.9346,4.9025,4.8778,4.8528,4.8272,4.8052,4.7875,4.7729,4.7607,4.7505,4.742,4.735,4.7293,4.7247,4.7209,4.7177,4.7149,4.7126,4.7105,4.7087,4.707,4.7056,4.7042,4.703,4.702,4.701,4.7001,4.6994,4.6987,4.6981,4.6975,4.697,4.6966,4.6962,4.6958,4.6955,4.6953,4.695,4.6948,4.6947,4.6945,4.6944,4.6942,4.6941,4.6941,4.694,4.6939,4.6939,4.6938,4.6938,4.6937,4.6937,4.6937,4.6937,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6936,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6935,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6934,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6933,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932,4.6932,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.6931,4.693,4.693,4.693,4.693,4.693,4.693,4.693,4.693,4.693,4.693,4.693,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6929,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6928,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6927,4.6926,4.6926,4.6926 9 | 13.123,5.5549,5.3761,5.1017,4.9838,4.932,4.9028,4.8825,4.8671,4.8548,4.8451,4.8375,4.8318,4.8276,4.8247,4.8225,4.8209,4.8197,4.8187,4.8179,4.8172,4.8165,4.8159,4.8154,4.8149,4.8144,4.8139,4.8135,4.813,4.8126,4.8122,4.8118,4.8114,4.811,4.8106,4.8102,4.8098,4.8093,4.8089,4.8085,4.808,4.8075,4.807,4.8065,4.806,4.8054,4.8048,4.8041,4.8034,4.8026,4.8018,4.8009,4.7999,4.7988,4.7977,4.7964,4.7949,4.7934,4.7916,4.7897,4.7875,4.7852,4.7826,4.7797,4.7765,4.7729,4.769,4.7647,4.7599,4.7548,4.7494,4.7439,4.7384,4.7332,4.7285,4.7241,4.7202,4.7168,4.7138,4.7112,4.709,4.707,4.7054,4.7039,4.7027,4.7016,4.7006,4.6997,4.699,4.6983,4.6977,4.6972,4.6967,4.6962,4.6958,4.6954,4.6951,4.6947,4.6944,4.6942,4.6939,4.6936,4.6934,4.6931,4.6929,4.6927,4.6925,4.6923,4.6921,4.6919,4.6917,4.6915,4.6914,4.6912,4.691,4.6908,4.6907,4.6905,4.6903,4.6902,4.69,4.6898,4.6896,4.6895,4.6893,4.6891,4.689,4.6888,4.6886,4.6884,4.6882,4.688,4.6878,4.6876,4.6874,4.6872,4.687,4.6868,4.6865,4.6863,4.6861,4.6858,4.6855,4.6852,4.6849,4.6846,4.6843,4.684,4.6836,4.6833,4.6829,4.6825,4.6821,4.6816,4.6811,4.6807,4.6802,4.6796,4.6791,4.6785,4.6779,4.6772,4.6766,4.6759,4.6752,4.6744,4.6737,4.6729,4.6722,4.6714,4.6706,4.6698,4.669,4.6682,4.6674,4.6666,4.6659,4.6651,4.6644,4.6637,4.663,4.6623,4.6616,4.661,4.6603,4.6597,4.6591,4.6585,4.6579,4.6574,4.6568,4.6563,4.6557,4.6552,4.6547,4.6542,4.6537,4.6533,4.6529,4.6524 10 | 17.652,5.6782,5.1402,5.0619,5.0237,4.9961,4.9754,4.9599,4.948,4.9386,4.9308,4.924,4.9179,4.9121,4.9064,4.9006,4.8946,4.8883,4.8818,4.8755,4.8694,4.8639,4.8589,4.8545,4.8506,4.8472,4.8442,4.8416,4.8393,4.8372,4.8354,4.8337,4.8322,4.8308,4.8295,4.8284,4.8273,4.8263,4.8253,4.8245,4.8236,4.8229,4.8221,4.8214,4.8207,4.8201,4.8195,4.8189,4.8183,4.8177,4.8172,4.8166,4.8161,4.8156,4.8151,4.8146,4.8142,4.8137,4.8132,4.8128,4.8123,4.8119,4.8114,4.811,4.8106,4.8102,4.8097,4.8093,4.8089,4.8085,4.808,4.8076,4.8072,4.8068,4.8064,4.8059,4.8055,4.8051,4.8047,4.8042,4.8038,4.8034,4.8029,4.8025,4.8021,4.8016,4.8012,4.8007,4.8003,4.7998,4.7994,4.7989,4.7984,4.798,4.7975,4.7971,4.7966,4.7962,4.7957,4.7953,4.7948,4.7944,4.7939,4.7935,4.7931,4.7927,4.7923,4.7919,4.7915,4.7912,4.7908,4.7905,4.7901,4.7898,4.7895,4.7892,4.789,4.7887,4.7885,4.7882,4.788,4.7878,4.7876,4.7874,4.7872,4.787,4.7869,4.7867,4.7866,4.7865,4.7863,4.7862,4.7861,4.786,4.7859,4.7858,4.7857,4.7856,4.7856,4.7855,4.7854,4.7853,4.7853,4.7852,4.7851,4.7851,4.785,4.785,4.7849,4.7849,4.7848,4.7848,4.7847,4.7847,4.7847,4.7846,4.7846,4.7845,4.7845,4.7844,4.7844,4.7844,4.7843,4.7843,4.7843,4.7842,4.7842,4.7842,4.7841,4.7841,4.7841,4.784,4.784,4.784,4.7839,4.7839,4.7839,4.7838,4.7838,4.7838,4.7838,4.7837,4.7837,4.7837,4.7836,4.7836,4.7836,4.7836,4.7835,4.7835,4.7835,4.7834,4.7834,4.7834,4.7834,4.7833,4.7833,4.7833,4.7833,4.7832 11 | -1.3898e+07,-1.2267e+07,-1.212e+07,-1.211e+07,-1.2107e+07,-1.2105e+07,-1.2103e+07,-1.2101e+07,-1.2099e+07,-1.2097e+07,-1.2095e+07,-1.2094e+07,-1.2093e+07,-1.2092e+07,-1.2091e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2086e+07,-1.2086e+07,-1.2086e+07,-1.2086e+07,-1.2086e+07,-1.2085e+07,-1.2085e+07,-1.2085e+07,-1.2085e+07,-1.2085e+07,-1.2085e+07,-1.2084e+07,-1.2084e+07,-1.2084e+07,-1.2084e+07,-1.2084e+07,-1.2084e+07,-1.2083e+07,-1.2083e+07,-1.2083e+07,-1.2083e+07,-1.2083e+07,-1.2083e+07,-1.2083e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2082e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2073e+07,-1.2073e+07,-1.2072e+07,-1.2071e+07,-1.207e+07,-1.2068e+07,-1.2066e+07,-1.2063e+07,-1.2061e+07,-1.2059e+07,-1.2057e+07 12 | -1.3464e+07,-1.2235e+07,-1.2183e+07,-1.2175e+07,-1.2167e+07,-1.2156e+07,-1.2141e+07,-1.2131e+07,-1.2125e+07,-1.2121e+07,-1.2119e+07,-1.2117e+07,-1.2115e+07,-1.2114e+07,-1.2113e+07,-1.2112e+07,-1.2112e+07,-1.2111e+07,-1.211e+07,-1.211e+07,-1.211e+07,-1.2109e+07,-1.2109e+07,-1.2108e+07,-1.2108e+07,-1.2107e+07,-1.2107e+07,-1.2106e+07,-1.2105e+07,-1.2104e+07,-1.2103e+07,-1.2102e+07,-1.2101e+07,-1.21e+07,-1.2099e+07,-1.2098e+07,-1.2098e+07,-1.2097e+07,-1.2097e+07,-1.2096e+07,-1.2096e+07,-1.2096e+07,-1.2096e+07,-1.2096e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2095e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2094e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2093e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2092e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.2091e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.209e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2089e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2088e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07,-1.2087e+07 13 | -1.3454e+07,-1.2249e+07,-1.2172e+07,-1.2161e+07,-1.2157e+07,-1.2153e+07,-1.2149e+07,-1.2144e+07,-1.2138e+07,-1.2131e+07,-1.2124e+07,-1.2117e+07,-1.2109e+07,-1.2103e+07,-1.2096e+07,-1.2091e+07,-1.2087e+07,-1.2084e+07,-1.2081e+07,-1.208e+07,-1.2079e+07,-1.2078e+07,-1.2077e+07,-1.2077e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2075e+07,-1.2075e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2068e+07,-1.2068e+07,-1.2068e+07,-1.2068e+07,-1.2068e+07,-1.2068e+07,-1.2067e+07,-1.2067e+07,-1.2067e+07,-1.2067e+07,-1.2067e+07,-1.2067e+07,-1.2066e+07,-1.2066e+07,-1.2066e+07,-1.2066e+07,-1.2065e+07,-1.2065e+07,-1.2065e+07,-1.2065e+07,-1.2065e+07,-1.2064e+07,-1.2064e+07,-1.2064e+07,-1.2064e+07,-1.2064e+07,-1.2064e+07,-1.2063e+07,-1.2063e+07,-1.2063e+07,-1.2063e+07,-1.2063e+07,-1.2063e+07,-1.2062e+07,-1.2062e+07,-1.2062e+07,-1.2062e+07,-1.2062e+07,-1.2062e+07,-1.2062e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.2061e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.206e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2059e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2058e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2057e+07,-1.2056e+07,-1.2056e+07,-1.2056e+07,-1.2056e+07,-1.2056e+07,-1.2056e+07,-1.2055e+07,-1.2055e+07,-1.2055e+07,-1.2055e+07,-1.2055e+07,-1.2054e+07,-1.2054e+07,-1.2054e+07,-1.2054e+07,-1.2053e+07,-1.2053e+07,-1.2053e+07,-1.2053e+07,-1.2052e+07,-1.2052e+07,-1.2052e+07,-1.2052e+07,-1.2052e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.2051e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07 14 | -1.3809e+07,-1.2268e+07,-1.215e+07,-1.2136e+07,-1.2125e+07,-1.2113e+07,-1.2103e+07,-1.2094e+07,-1.2088e+07,-1.2083e+07,-1.2079e+07,-1.2076e+07,-1.2074e+07,-1.2072e+07,-1.207e+07,-1.2069e+07,-1.2067e+07,-1.2065e+07,-1.2063e+07,-1.2061e+07,-1.2059e+07,-1.2057e+07,-1.2055e+07,-1.2053e+07,-1.2052e+07,-1.205e+07,-1.2049e+07,-1.2048e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07 15 | -1.3378e+07,-1.2231e+07,-1.214e+07,-1.2114e+07,-1.2097e+07,-1.2083e+07,-1.2073e+07,-1.2066e+07,-1.2062e+07,-1.2058e+07,-1.2056e+07,-1.2054e+07,-1.2053e+07,-1.2052e+07,-1.2051e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07 16 | -1.3431e+07,-1.2172e+07,-1.2108e+07,-1.2094e+07,-1.2085e+07,-1.2076e+07,-1.2069e+07,-1.2064e+07,-1.2059e+07,-1.2057e+07,-1.2054e+07,-1.2053e+07,-1.2052e+07,-1.2051e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2037e+07,-1.2037e+07,-1.2037e+07,-1.2037e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07 17 | -1.3699e+07,-1.2182e+07,-1.2127e+07,-1.2121e+07,-1.2118e+07,-1.2116e+07,-1.2113e+07,-1.2109e+07,-1.2105e+07,-1.2099e+07,-1.2093e+07,-1.2087e+07,-1.2082e+07,-1.2076e+07,-1.207e+07,-1.2065e+07,-1.2061e+07,-1.2058e+07,-1.2055e+07,-1.2054e+07,-1.2052e+07,-1.2051e+07,-1.2051e+07,-1.205e+07,-1.205e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07 18 | -1.3664e+07,-1.227e+07,-1.2205e+07,-1.2183e+07,-1.2171e+07,-1.216e+07,-1.2149e+07,-1.2125e+07,-1.2105e+07,-1.2098e+07,-1.2092e+07,-1.2085e+07,-1.2078e+07,-1.2072e+07,-1.2068e+07,-1.2064e+07,-1.2062e+07,-1.2059e+07,-1.2057e+07,-1.2056e+07,-1.2054e+07,-1.2053e+07,-1.2053e+07,-1.2052e+07,-1.2051e+07,-1.2051e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2048e+07,-1.2048e+07,-1.2047e+07,-1.2047e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07 19 | -1.3342e+07,-1.2271e+07,-1.2209e+07,-1.2143e+07,-1.2114e+07,-1.2103e+07,-1.2096e+07,-1.2092e+07,-1.2088e+07,-1.2085e+07,-1.2083e+07,-1.2081e+07,-1.208e+07,-1.2079e+07,-1.2078e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.2069e+07,-1.2069e+07,-1.2068e+07,-1.2067e+07,-1.2066e+07,-1.2066e+07,-1.2065e+07,-1.2064e+07,-1.2062e+07,-1.2061e+07,-1.206e+07,-1.2058e+07,-1.2057e+07,-1.2055e+07,-1.2054e+07,-1.2053e+07,-1.2051e+07,-1.205e+07,-1.2049e+07,-1.2049e+07,-1.2048e+07,-1.2047e+07,-1.2047e+07,-1.2046e+07,-1.2046e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2045e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2044e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2043e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2042e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.2041e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.204e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2039e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2038e+07,-1.2037e+07,-1.2037e+07,-1.2037e+07,-1.2037e+07,-1.2037e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2036e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2035e+07,-1.2034e+07,-1.2034e+07,-1.2034e+07,-1.2033e+07,-1.2033e+07,-1.2033e+07,-1.2033e+07,-1.2032e+07,-1.2032e+07,-1.2032e+07,-1.2032e+07,-1.2032e+07,-1.2031e+07,-1.2031e+07,-1.2031e+07,-1.2031e+07,-1.2031e+07,-1.203e+07,-1.203e+07,-1.203e+07,-1.203e+07,-1.203e+07,-1.2029e+07,-1.2029e+07,-1.2029e+07,-1.2029e+07,-1.2029e+07,-1.2029e+07,-1.2029e+07 20 | -1.371e+07,-1.2279e+07,-1.2157e+07,-1.2141e+07,-1.2132e+07,-1.2125e+07,-1.2121e+07,-1.2118e+07,-1.2115e+07,-1.2113e+07,-1.2111e+07,-1.211e+07,-1.2108e+07,-1.2106e+07,-1.2105e+07,-1.2103e+07,-1.2101e+07,-1.2099e+07,-1.2097e+07,-1.2096e+07,-1.2094e+07,-1.2092e+07,-1.2091e+07,-1.2089e+07,-1.2088e+07,-1.2087e+07,-1.2086e+07,-1.2085e+07,-1.2084e+07,-1.2084e+07,-1.2083e+07,-1.2083e+07,-1.2082e+07,-1.2082e+07,-1.2081e+07,-1.2081e+07,-1.2081e+07,-1.208e+07,-1.208e+07,-1.208e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2079e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2078e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2077e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2076e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2075e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2074e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2073e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2072e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.2071e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.207e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07,-1.2069e+07 21 | -------------------------------------------------------------------------------- /curves/mu10.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/mu10.mat -------------------------------------------------------------------------------- /curves/mu_curve10.m: -------------------------------------------------------------------------------- 1 | % Load data "mu10.mat". 2 | load mu10; 3 | 4 | % Create a figure. 5 | figure('DefaultAxesFontSize',16); 6 | ax = gca; 7 | 8 | % Plot the mu curves. 9 | xbound = [1,1,1:500,500,500,fliplr(1:500)]; 10 | ybound = [mean(mu(:,1))+std(mu(:,1)),mean(mu(:,1))-std(mu(:,1)),mean(mu)-std(mu),... 11 | mean(mu(:,end))-std(mu(:,end)),mean(mu(:,end))+std(mu(:,end)),fliplr(mean(mu)+std(mu))]; 12 | hold on; 13 | plot(mu','-','Color',[0.76,0.87,0.78]); 14 | hold on; 15 | s = fill(xbound,ybound,'y','EdgeColor','white'); 16 | alpha(s,0.20); 17 | hold on; 18 | h = plot(mean(mu),'Color','r','LineWidth',2); 19 | hold on; 20 | p = plot(39.01*ones(500,1),'-.','Color','magenta','LineWidth',2); 21 | 22 | % Change settings and items. 23 | box on; 24 | grid on; 25 | ax.XTick = [0,50,100,200,500]; 26 | ax.XTickLabel = [0,50,100,200,500]; 27 | ylim([38,40.5]); 28 | ax.YTick = [38,38.5,38.75,39,39.25,39.50,39.75,40.5]; 29 | ax.YTickLabel = ['38.00';'38.50';'38.75';'39.00';'39.25';'39.50';'39.75';'40.50']; 30 | ax.GridLineStyle = '-.'; 31 | xlabel('Iteration'); 32 | ylabel('Speed (km/h)'); 33 | legend([p,h],{'Mean of Observations', 'Model Parameter \mu'}); 34 | 35 | % Save the figure as "mu_curve10.pdf". 36 | set(gcf, 'PaperSize', [6 4.5]); 37 | set(gcf, 'PaperPositionMode', 'manual'); 38 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 39 | saveas(gcf,'mu_curve10','pdf'); -------------------------------------------------------------------------------- /curves/mu_curve10.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/mu_curve10.pdf -------------------------------------------------------------------------------- /curves/mu_curve10.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/mu_curve10.png -------------------------------------------------------------------------------- /curves/overfitting.m: -------------------------------------------------------------------------------- 1 | %% overfitting #1 2 | data = importdata('BCPF_fiber_rmselb_m30_r5.csv'); 3 | rmse = data(1:10,:); 4 | lb = data(11:20,:); 5 | 6 | color = [55,126,184]/255; 7 | gray = [211,211,211]/255; 8 | orange = [ 0.91,0.41,0.17]; 9 | blue = [0,0.4470,0.7410]; 10 | light_blue = [0.5843,0.8157,0.9882]; 11 | light_red = [255,204,204]/255; 12 | 13 | figure('DefaultAxesFontSize',16); 14 | ax = gca; 15 | 16 | yyaxis left; 17 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 18 | ybound = [max(rmse(:,1)),min(rmse(:,1)),min(rmse),min(rmse(:,end)),max(rmse(:,end)),fliplr(max(rmse))]; 19 | hold on; 20 | s = fill(xbound,ybound,color,'EdgeColor',gray,'linewidth',0.5); 21 | alpha(s,0.10); 22 | hold on; 23 | plot(rmse','-','color',light_blue,'linewidth',0.5); 24 | hold on; 25 | h1 = plot(mean(rmse),'-','color',blue,'linewidth',2,'DisplayName','rmse'); 26 | 27 | box on; 28 | grid on; 29 | ax.XTick = [0,25,50,100,200]; 30 | ax.XTickLabel = [0,25,50,100,200]; 31 | ylim([4.60,5.60]); 32 | ax.YTick = [4.60,4.70,4.80,4.90,5.10,5.60]; 33 | ax.YTickLabel = ['4.60';'4.70';'4.80';'4.90';'5.10';'5.60']; 34 | ax.GridLineStyle = '-.'; 35 | xlabel('Epoch'); 36 | ylabel('RMSE (km/h)'); 37 | 38 | yyaxis right; 39 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 40 | ybound = [max(lb(:,1)),min(lb(:,1)),min(lb),min(lb(:,end)),max(lb(:,end)),fliplr(max(lb))]; 41 | hold on; 42 | s = fill(xbound,ybound,gray,'EdgeColor',gray,'linewidth',0.5); 43 | alpha(s,0.10); 44 | hold on; 45 | plot(lb','-','color',light_red,'linewidth',0.5); 46 | hold on; 47 | h2 = plot(mean(lb),'-','color',orange,'linewidth',2,'DisplayName','lb'); 48 | ylim([-1.24e7,-1.19e7]); 49 | ylabel('Lower bound'); 50 | 51 | legend([h1,h2],{'Average of RMSE';'Average of lower bound'},'location','north'); 52 | 53 | set(gcf, 'PaperSize', [6 4.5]); 54 | set(gcf, 'PaperPositionMode', 'manual'); 55 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 56 | saveas(gcf,'overfitting_ms30_r5','pdf'); 57 | 58 | %% overfitting #2 59 | data = importdata('BCPF_fiber_rmselb_m30_r10.csv'); 60 | rmse = data(1:10,:); 61 | lb = data(11:20,:); 62 | 63 | figure('DefaultAxesFontSize',16); 64 | ax = gca; 65 | 66 | yyaxis left; 67 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 68 | ybound = [max(rmse(:,1)),min(rmse(:,1)),min(rmse),min(rmse(:,end)),max(rmse(:,end)),fliplr(max(rmse))]; 69 | hold on; 70 | s = fill(xbound,ybound,color,'EdgeColor',gray,'linewidth',0.5); 71 | alpha(s,0.10); 72 | hold on; 73 | plot(rmse','-','color',light_blue,'linewidth',0.5); 74 | hold on; 75 | h1 = plot(mean(rmse),'-','color',blue,'linewidth',2,'DisplayName','rmse'); 76 | 77 | box on; 78 | grid on; 79 | ax.XTick = [0,25,50,100,200]; 80 | ax.XTickLabel = [0,25,50,100,200]; 81 | ylim([4.35,8.60]); 82 | ax.YTick = [4.35,4.60,4.90,5.20,6.00,8.60]; 83 | ax.YTickLabel = ['4.35';'4.60';'4.90';'5.20';'6.00';'8.60']; 84 | ax.GridLineStyle = '-.'; 85 | xlabel('Epoch'); 86 | ylabel('RMSE (km/h)'); 87 | 88 | yyaxis right; 89 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 90 | ybound = [max(lb(:,1)),min(lb(:,1)),min(lb),min(lb(:,end)),max(lb(:,end)),fliplr(max(lb))]; 91 | hold on; 92 | s = fill(xbound,ybound,gray,'EdgeColor',gray,'linewidth',0.5); 93 | alpha(s,0.10); 94 | hold on; 95 | plot(lb','-','color',light_red,'linewidth',0.5); 96 | hold on; 97 | h2 = plot(mean(lb),'-','color',orange,'linewidth',2,'DisplayName','lb'); 98 | ylim([-1.205e7,-1.19e7]); 99 | ylabel('Lower bound'); 100 | 101 | legend([h1,h2],{'Average of RMSE';'Average of lower bound'},'location','north'); 102 | 103 | set(gcf, 'PaperSize', [6 4.5]); 104 | set(gcf, 'PaperPositionMode', 'manual'); 105 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 106 | saveas(gcf,'overfitting_ms30_r10','pdf'); -------------------------------------------------------------------------------- /curves/overfitting_ms30_r10.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/overfitting_ms30_r10.pdf -------------------------------------------------------------------------------- /curves/overfitting_ms30_r10.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/overfitting_ms30_r10.png -------------------------------------------------------------------------------- /curves/overfitting_ms30_r5.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/overfitting_ms30_r5.pdf -------------------------------------------------------------------------------- /curves/overfitting_ms30_r5.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/curves/overfitting_ms30_r5.png -------------------------------------------------------------------------------- /heat-maps/bias10.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/bias10.mat -------------------------------------------------------------------------------- /heat-maps/bias20.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/bias20.mat -------------------------------------------------------------------------------- /heat-maps/bias60.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/bias60.mat -------------------------------------------------------------------------------- /heat-maps/factor2.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/factor2.pdf -------------------------------------------------------------------------------- /heat-maps/factor2.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/factor2.png -------------------------------------------------------------------------------- /heat-maps/factor3.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/factor3.pdf -------------------------------------------------------------------------------- /heat-maps/factor3.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/factor3.png -------------------------------------------------------------------------------- /heat-maps/heat_map10.m: -------------------------------------------------------------------------------- 1 | % Load data "bias10.mat". 2 | load bias10; 3 | 4 | % Create a figure. 5 | figure('DefaultAxesFontSize',16); 6 | ax = gca; 7 | 8 | % Create a heatmap. 9 | mat = bias{2}.*ones(61,144)+bias{3}'.*ones(61,144); 10 | imagesc(mat); 11 | c = colorbar; 12 | colormap hot; 13 | ylabel(c,'Speed (km/h)'); 14 | ax.XTick = [1,13,25,37,49,61,73,85,97,109,121,133,144]; 15 | ax.XTickLabel = ['00:00';'02:00';'04:00';'06:00';'08:00';'10:00';'12:00';... 16 | '14:00';'16:00';'18:00';'20:00';'22:00';'23:50']; 17 | xtickangle(45); 18 | ax.YTick = [1,6,11,16,21,26,31,36,41,46,51,56,61]; 19 | ax.YTickLabel = ['08-01';'08-06';'08-11';'08-16';'08-21';'08-26';'08-31';... 20 | '09-05';'09-10';'09-15';'09-20';'09-25';'09-30']; 21 | 22 | % Save the figure as "heat_map10.pdf". 23 | set(gcf, 'PaperSize', [6 4.5]); 24 | set(gcf, 'PaperPositionMode', 'manual'); 25 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 26 | saveas(gcf,'heat_map10','pdf'); -------------------------------------------------------------------------------- /heat-maps/heat_map10.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/heat_map10.pdf -------------------------------------------------------------------------------- /heat-maps/heat_map10.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/heat_map10.png -------------------------------------------------------------------------------- /heat-maps/heat_map20.m: -------------------------------------------------------------------------------- 1 | % Load data "bias20.mat". 2 | clc; 3 | clear; 4 | load bias20; 5 | 6 | % Create a figure. 7 | figure('DefaultAxesFontSize',16); 8 | ax = gca; 9 | 10 | % Create a heatmap. 11 | mat = bias{2}.*ones(61,144)+bias{3}'.*ones(61,144); 12 | imagesc(mat); 13 | c = colorbar; 14 | colormap hot; 15 | ylabel(c,'Speed (km/h)'); 16 | ax.XTick = [1,13,25,37,49,61,73,85,97,109,121,133,144]; 17 | ax.XTickLabel = ['00:00';'02:00';'04:00';'06:00';'08:00';'10:00';'12:00';... 18 | '14:00';'16:00';'18:00';'20:00';'22:00';'23:50']; 19 | xtickangle(45); 20 | ax.YTick = [1,6,11,16,21,26,31,36,41,46,51,56,61]; 21 | ax.YTickLabel = ['08-01';'08-06';'08-11';'08-16';'08-21';'08-26';'08-31';... 22 | '09-05';'09-10';'09-15';'09-20';'09-25';'09-30']; 23 | 24 | % Save the figure as "heat_map20.pdf". 25 | set(gcf, 'PaperSize', [6 4.5]); 26 | set(gcf, 'PaperPositionMode', 'manual'); 27 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 28 | saveas(gcf,'heat_map20','pdf'); -------------------------------------------------------------------------------- /heat-maps/heat_map20.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/heat_map20.pdf -------------------------------------------------------------------------------- /heat-maps/heat_map60.m: -------------------------------------------------------------------------------- 1 | % Load data "bias60.mat". 2 | clc; 3 | clear; 4 | load bias60; 5 | 6 | % Create a figure. 7 | figure('DefaultAxesFontSize',16); 8 | ax = gca; 9 | 10 | % Create a heatmap. 11 | mat = bias{2}.*ones(61,144)+bias{3}'.*ones(61,144); 12 | imagesc(mat); 13 | c = colorbar; 14 | colormap hot; 15 | ylabel(c,'Speed (km/h)'); 16 | ax.XTick = [1,13,25,37,49,61,73,85,97,109,121,133,144]; 17 | ax.XTickLabel = ['00:00';'02:00';'04:00';'06:00';'08:00';'10:00';'12:00';... 18 | '14:00';'16:00';'18:00';'20:00';'22:00';'23:50']; 19 | xtickangle(45); 20 | ax.YTick = [1,6,11,16,21,26,31,36,41,46,51,56,61]; 21 | ax.YTickLabel = ['08-01';'08-06';'08-11';'08-16';'08-21';'08-26';'08-31';... 22 | '09-05';'09-10';'09-15';'09-20';'09-25';'09-30']; 23 | 24 | % Save the figure as "heat_map60.pdf". 25 | set(gcf, 'PaperSize', [6 4.5]); 26 | set(gcf, 'PaperPositionMode', 'manual'); 27 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 28 | saveas(gcf,'heat_map60','pdf'); -------------------------------------------------------------------------------- /heat-maps/heat_map60.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/heat_map60.pdf -------------------------------------------------------------------------------- /heat-maps/latent_factors.m: -------------------------------------------------------------------------------- 1 | load latent_factors; 2 | 3 | % figure #1 (factor2) 4 | figure('DefaultAxesFontSize',16); 5 | ax = gca; 6 | 7 | imagesc(factor2); 8 | colormap hot; 9 | ax.YTick = [1,6,11,16,21,26,31,36,41,46,51,56,61]; 10 | ax.YTickLabel = ['08-01';'08-06';'08-11';'08-16';'08-21';'08-26';'08-31';... 11 | '09-05';'09-10';'09-15';'09-20';'09-25';'09-30']; 12 | 13 | set(gcf, 'PaperSize', [3 10]); 14 | set(gcf, 'PaperPositionMode', 'manual'); 15 | set(gcf, 'PaperPosition', [0 0 3 10]); 16 | saveas(gcf,'factor2','pdf'); 17 | 18 | % figure #2 (factor3) 19 | figure('DefaultAxesFontSize',16); 20 | ax = gca; 21 | 22 | imagesc(factor3); 23 | c = colorbar; 24 | colormap hot; 25 | ax.YTick = [1,13,25,37,49,61,73,85,97,109,121,133,144]; 26 | ax.YTickLabel = ['00:00';'02:00';'04:00';'06:00';'08:00';'10:00';'12:00';... 27 | '14:00';'16:00';'18:00';'20:00';'22:00';'23:50']; 28 | 29 | set(gcf, 'PaperSize', [3 10]); 30 | set(gcf, 'PaperPositionMode', 'manual'); 31 | set(gcf, 'PaperPosition', [0 0 3 10]); 32 | saveas(gcf,'factor3','pdf'); -------------------------------------------------------------------------------- /heat-maps/latent_factors.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/heat-maps/latent_factors.mat -------------------------------------------------------------------------------- /rmse-curves/rmse.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse.mat -------------------------------------------------------------------------------- /rmse-curves/rmse10.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse10.mat -------------------------------------------------------------------------------- /rmse-curves/rmse_curve.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse_curve.pdf -------------------------------------------------------------------------------- /rmse-curves/rmse_curve.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse_curve.png -------------------------------------------------------------------------------- /rmse-curves/rmse_curve10.m: -------------------------------------------------------------------------------- 1 | load rmse10; 2 | 3 | % Set colors. 4 | color = [55,126,184]/255; 5 | 6 | % Create a figure. 7 | figure('DefaultAxesFontSize',16); 8 | ax = gca; 9 | 10 | % Plot the RMSE curves. 11 | xbound = [1,1,1:500,500,500,fliplr(1:500)]; 12 | ybound = [max(rmse10(:,1)),min(rmse10(:,1)),min(rmse10),min(rmse10(:,end)),max(rmse10(:,end)),fliplr(max(rmse10))]; 13 | hold on; 14 | plot(rmse10'); 15 | hold on; 16 | s = fill(xbound,ybound,'y','EdgeColor',color,'linewidth',1); 17 | alpha(s,0.1); 18 | hold on; 19 | h = plot(mean(rmse10),'color','r','linewidth',3,'DisplayName','method #2'); 20 | 21 | % Change settings and items. 22 | box on; 23 | grid on; 24 | ax.XTick = [0,50,100,200,500]; 25 | ax.XTickLabel = [0,50,100,200,500]; 26 | ylim([4.10,4.40]); 27 | ax.YTick = [4.10,4.12,4.14,4.16,4.18,4.20,4.25,4.40]; 28 | ax.YTickLabel = ['4.10';'4.12';'4.14';'4.16';'4.18';'4.20';'4.25';'4.40']; 29 | ax.GridLineStyle = '-.'; 30 | xlabel('Iteration'); 31 | ylabel('RMSE (km/h)'); 32 | legend([h,s],{'Averaged Curve';'Curves Boundary'}); 33 | 34 | % Save the figure as "rmse_curve10.pdf". 35 | set(gcf, 'PaperSize', [6 4.5]); 36 | set(gcf, 'PaperPositionMode', 'manual'); 37 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 38 | saveas(gcf,'rmse_curve10','pdf'); -------------------------------------------------------------------------------- /rmse-curves/rmse_curve10.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse_curve10.pdf -------------------------------------------------------------------------------- /rmse-curves/rmse_curve10.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/rmse-curves/rmse_curve10.png -------------------------------------------------------------------------------- /rmse-curves/rmse_curves.m: -------------------------------------------------------------------------------- 1 | % Load data "rmse.mat". 2 | load rmse; 3 | 4 | % Set colors. 5 | color1 = [55,126,184]/255; 6 | color2 = [1,0.6,0.78]; 7 | color3 = [255,0,255]/255; 8 | 9 | % Create a figure. 10 | figure('DefaultAxesFontSize',14); 11 | ax = gca; 12 | 13 | % Plot the first RMSE curve. 14 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 15 | ybound = [max(rmse{1}(:,1)),min(rmse{1}(:,1)),min(rmse{1}),min(rmse{1}(:,end)),max(rmse{1}(:,end)),fliplr(max(rmse{1}))]; 16 | s = fill(xbound,ybound,'y','EdgeColor',color1,'linewidth',0.5); 17 | alpha(s,0.1); 18 | hold on; 19 | h1 = plot(mean(rmse{1}),'color','r','linewidth',2,'DisplayName','method #1'); 20 | 21 | % Plot the second RMSE curve. 22 | xbound = [1,1,1:200,200,200,fliplr(1:200)]; 23 | ybound = [max(rmse{2}(:,1)),min(rmse{2}(:,1)),min(rmse{2}),min(rmse{2}(:,end)),max(rmse{2}(:,end)),fliplr(max(rmse{2}))]; 24 | hold on; 25 | s = fill(xbound,ybound,color2,'EdgeColor',color1,'linewidth',0.5); 26 | alpha(s,0.1); 27 | hold on; 28 | h2 = plot(mean(rmse{2}),'color',color3,'linewidth',2,'DisplayName','method #2'); 29 | 30 | % Change settings and items. 31 | grid on; 32 | ax.XTick = [0,25,50,100,200]; 33 | ax.XTickLabel = [0,25,50,100,200]; 34 | ylim([3.95,4.15]); 35 | ax.YTick = [3.97,3.98,3.99,4.01,4.02,4.03,4.15]; 36 | ax.YTickLabel = [3.97,3.98,3.99,4.01,4.02,4.03,4.15]; 37 | ax.GridLineStyle = '-.'; 38 | xlabel('Iteration'); 39 | ylabel('RMSE (km/h)'); 40 | legend([h1,h2],{'method #1','method #2'}); 41 | 42 | % Save a figure in a format of "rmse_curve.pdf". 43 | set(gcf, 'PaperSize', [6 4.5]); 44 | set(gcf, 'PaperPositionMode', 'manual'); 45 | set(gcf, 'PaperPosition', [0 0 6 4.5]); 46 | saveas(gcf,'rmse_curve','pdf'); -------------------------------------------------------------------------------- /time-series/nyc_data_competness.xls: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/nyc_data_competness.xls -------------------------------------------------------------------------------- /time-series/nyc_data_completeness.m: -------------------------------------------------------------------------------- 1 | clear; 2 | data = load('data_completeness.mat'); 3 | data = data.data; 4 | 5 | orange = [ 0.91,0.41,0.17]; 6 | blue = [0,0.4470,0.7410]; 7 | figure('DefaultAxesFontSize',14); 8 | ax = gca; 9 | xlabel('Month'); 10 | 11 | yyaxis left; 12 | h1 = plot(data(:,1),'-','color',blue,'linewidth',2,'DisplayName','rmse'); 13 | ylabel('Data amount (\times 10^6)'); 14 | ylim([0,1.6e+6]); 15 | ax.YTick = 0:0.4e+6:1.6e+6; 16 | ax.YTickLabel = ['0 ';'0.4';'0.8';'1.2';'1.6']; 17 | yyaxis right; 18 | h2 = plot(data(:,2),'-','color',orange,'linewidth',2,'DisplayName','lb'); 19 | ylabel('Completeness'); 20 | ax.YTick = 0:0.25:1; 21 | ax.YTickLabel = 0:0.25:1; 22 | 23 | box on; 24 | grid on; 25 | ax.XTick = 1:4:48; 26 | ax.XTickLabel = ['Apr.';'Aug.';'Dec.';'Apr.';'Aug.';'Dec.';'Apr.';'Aug.';'Dec.';'Apr.';'Aug.';'Dec.';]; 27 | ax.GridLineStyle = '-.'; 28 | set(gcf, 'PaperSize', [8 3]); 29 | set(gcf, 'PaperPositionMode', 'manual'); 30 | set(gcf, 'PaperPosition', [0 0 8 3]); 31 | saveas(gcf,'nyc_data_completeness','pdf'); -------------------------------------------------------------------------------- /time-series/nyc_data_completeness.mat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/nyc_data_completeness.mat -------------------------------------------------------------------------------- /time-series/nyc_data_completeness.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/nyc_data_completeness.pdf -------------------------------------------------------------------------------- /time-series/nyc_data_completeness.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/nyc_data_completeness.png -------------------------------------------------------------------------------- /time-series/speed_curve.m: -------------------------------------------------------------------------------- 1 | clear; 2 | data = importdata('road1_fiber_ms50_r10.csv'); 3 | data = data.data; 4 | 5 | figure('DefaultAxesFontSize',16); 6 | ax = gca; 7 | 8 | hold on;plot(1:144,data(1:144,4),'k-.','linewidth',0.5); 9 | hold on;plot(144+1:2*144,data(144+1:2*144,4),'color',[0,0.4470,0.7410],'linewidth',2); 10 | hold on;plot(2*144+1:3*144,data(2*144+1:3*144,4),'k-.','linewidth',0.5); 11 | hold on;plot(3*144+1:5*144,data(3*144+1:5*144,4),'color',[0,0.4470,0.7410],'linewidth',2); 12 | hold on;plot(5*144+1:6*144,data(5*144+1:6*144,4),'k-.','linewidth',0.5); 13 | hold on;plot(6*144+1:7*144,data(6*144+1:7*144,4),'color',[0,0.4470,0.7410],'linewidth',2); 14 | 15 | ax.XTick = [144,2*144,3*144,4*144,5*144,6*144,7*144]; 16 | ax.XTickLabel = [' ';' ';' ';' ';' ';' ';' ']; 17 | ax.YTick = [0,60]; 18 | ax.YTickLabel = [' ';' ']; 19 | x = [144,2*144,2*144,144]; 20 | y = [10,10,50,50]; 21 | 22 | v1 = [144,0; 2*144,0; 2*144,60; 144,60]; 23 | v2 = v1; v2(:,1)=v1(:,1)+2*144; 24 | v3 = v2; v3(:,1)=v2(:,1)+144; 25 | v4 = v3; v4(:,1)=v3(:,1)+2*144; 26 | f = [1,2,3,4]; 27 | patch('faces',f,'vertices',v1,'edgecolor',[0.5 0.5 0.5],'facecolor','green','facealpha',0.1); 28 | patch('faces',f,'vertices',v2,'edgecolor',[0.5 0.5 0.5],'facecolor','green','facealpha',0.1); 29 | patch('faces',f,'vertices',v3,'edgecolor',[0.5 0.5 0.5],'facecolor','green','facealpha',0.1); 30 | patch('faces',f,'vertices',v4,'edgecolor',[0.5 0.5 0.5],'facecolor','green','facealpha',0.1); 31 | xlim([10,50]); 32 | xlim([0,7*144]); 33 | 34 | grid on; 35 | box off; 36 | 37 | set(gcf, 'PaperSize', [6 2]); 38 | set(gcf, 'PaperPositionMode', 'manual'); 39 | set(gcf, 'PaperPosition', [0 0 6 2]); 40 | saveas(gcf,'curve1','pdf'); -------------------------------------------------------------------------------- /time-series/speed_curve.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/speed_curve.pdf -------------------------------------------------------------------------------- /time-series/speed_curve.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/speed_curve.png -------------------------------------------------------------------------------- /time-series/time_series_speed1.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/time_series_speed1.pdf -------------------------------------------------------------------------------- /time-series/time_series_speed1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/time_series_speed1.png -------------------------------------------------------------------------------- /time-series/time_series_speed1.py: -------------------------------------------------------------------------------- 1 | import pandas as pd 2 | import matplotlib.pyplot as plt 3 | import matplotlib.patches as patches 4 | import numpy as np 5 | 6 | obs = pd.read_csv('road1_fiber_ms50_r10.csv') 7 | data = obs.values 8 | plt.style.use('classic') 9 | fig = plt.figure(figsize=(13, 2.5)) 10 | ax = fig.add_axes([0.05,0.22,0.92,0.75]) 11 | 12 | plt.plot(data[0:14*144-1,0],color="#006ea3",linewidth=0.4,label="Actual value") 13 | plt.plot(data[0:14*144-1,2],color="#f79b5f",linewidth=1.8,label=r"Global parameter + biases") 14 | plt.plot(data[0:14*144-1,3],color="#e3120b",linewidth=1.0,label="Estimated value") 15 | 16 | ax.set_xlim([0, 14*144-1]) 17 | ax.set_ylim([0, 60]) 18 | ax.grid(color='gray', linestyle='-', linewidth=0.1, alpha=0.2) 19 | 20 | for i in range(14): 21 | if data[144*i,1] > 0: 22 | someX, someY = i*144, 0 23 | currentAxis = plt.gca() 24 | ax.add_patch(patches.Rectangle((someX, someY), 144, 60, 25 | alpha=0.1, facecolor='green')) 26 | 27 | plt.xticks(np.arange(0, 14*144, 72),["00:00", "12:00", "00:00", "12:00", 28 | "00:00", "12:00", "00:00", "12:00", 29 | "00:00", "12:00", "00:00", "12:00", 30 | "00:00", "12:00", "00:00", "12:00", 31 | "00:00", "12:00", "00:00", "12:00", 32 | "00:00", "12:00", "00:00", "12:00", 33 | "00:00", "12:00", "00:00", "12:00"], rotation=30) 34 | plt.yticks(np.arange(10,60,10),[10,20,30,40,50,60]) 35 | ax.set_ylabel("Speed (km/h)") 36 | plt.legend(ncol=3, loc=4) 37 | 38 | fig.savefig("time_series_speed1.pdf") -------------------------------------------------------------------------------- /time-series/time_series_speed2.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/time_series_speed2.pdf -------------------------------------------------------------------------------- /time-series/time_series_speed2.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/xinychen/academic-drawing/927d729e3f9115d7c7285d97c63cbb6c32cea449/time-series/time_series_speed2.png -------------------------------------------------------------------------------- /time-series/time_series_speed2.py: -------------------------------------------------------------------------------- 1 | import pandas as pd 2 | import matplotlib.pyplot as plt 3 | import matplotlib.patches as patches 4 | import numpy as np 5 | 6 | obs = pd.read_csv('road1_fiber_ms50_r10.csv') 7 | data = obs.values 8 | plt.style.use('classic') 9 | fig = plt.figure(figsize=(13, 2.5)) 10 | ax = fig.add_axes([0.05,0.22,0.92,0.75]) 11 | 12 | plt.plot(data[42*144:56*144-1,0],color="#006ea3",linewidth=0.4,label="Actual value") 13 | plt.plot(data[42*144:56*144-1,2],color="#f79b5f",linewidth=1.8,label=r"Global parameter + biases") 14 | plt.plot(data[42*144:56*144-1,3],color="#e3120b",linewidth=1.0,label="Estimated value") 15 | 16 | ax.set_xlim([0, 14*144-1]) 17 | ax.set_ylim([0, 60]) 18 | ax.grid(color='gray', linestyle='-', linewidth=0.1, alpha=0.2) 19 | 20 | for i in range(42,56): 21 | if data[144*i,1] > 0: 22 | someX, someY = (i-42)*144, 0 23 | currentAxis = plt.gca() 24 | ax.add_patch(patches.Rectangle((someX, someY), 144, 60, 25 | alpha=0.1, facecolor='green')) 26 | 27 | plt.xticks(np.arange(0, 14*144, 72),["00:00", "12:00", "00:00", "12:00", 28 | "00:00", "12:00", "00:00", "12:00", 29 | "00:00", "12:00", "00:00", "12:00", 30 | "00:00", "12:00", "00:00", "12:00", 31 | "00:00", "12:00", "00:00", "12:00", 32 | "00:00", "12:00", "00:00", "12:00", 33 | "00:00", "12:00", "00:00", "12:00"], rotation=30) 34 | plt.yticks(np.arange(10,60,10),[10,20,30,40,50,60]) 35 | ax.set_ylabel("Speed (km/h)") 36 | plt.legend(ncol=3, loc=4) 37 | 38 | fig.savefig("time_series_speed2.pdf") --------------------------------------------------------------------------------