-
Notifications
You must be signed in to change notification settings - Fork 4
/
aptC_plot.py
142 lines (138 loc) · 11.5 KB
/
aptC_plot.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
#!/usr/bin/env python2
import matplotlib
import matplotlib.pyplot as plt
xl = [0.000000,
0.010000,0.020000,0.030000,0.040000,0.050000,0.060000,0.070000,0.080000,
0.090000,0.100000,0.110000,0.120000,0.130000,0.140000,0.150000,0.160000,
0.170000,0.180000,0.190000,0.200000,0.210000,0.220000,0.230000,0.240000,
0.250000,0.260000,0.270000,0.280000,0.290000,0.300000,0.310000,0.320000,
0.330000,0.340000,0.350000,0.360000,0.370000,0.380000,0.390000,0.400000,
0.410000,0.420000,0.430000,0.440000,0.450000,0.460000,0.470000,0.480000,
0.490000,0.500000,0.510000,0.520000,0.530000,0.540000,0.550000,0.560000,
0.570000,0.580000,0.590000,0.600000,0.610000,0.620000,0.630000,0.640000,
0.650000,0.660000,0.670000,0.680000,0.690000,0.700000,0.710000,0.720000,
0.730000,0.740000,0.750000,0.760000,0.770000,0.780000,0.790000,0.800000,
0.810000,0.820000,0.830000,0.840000,0.850000,0.860000,0.870000,0.880000,
0.890000,0.900000,0.910000,0.920000,0.930000,0.940000,0.950000,0.960000,
0.970000,0.980000,0.990000,1.000000]
yl = [0.000000,
0.242842,0.338914,0.409651,0.466865,0.515209,0.557110,0.594036,0.626957,
0.656561,0.683357,0.707735,0.730004,0.750413,0.769169,0.786443,0.802383,
0.817113,0.830743,0.843368,0.855072,0.865931,0.876010,0.885371,0.894066,
0.902145,0.909653,0.916631,0.923115,0.929140,0.934738,0.939937,0.944765,
0.949247,0.953405,0.957262,0.960837,0.964149,0.967216,0.970054,0.972677,
0.975101,0.977338,0.979402,0.981304,0.983054,0.984664,0.986143,0.987499,
0.988743,0.989880,0.990920,0.991869,0.992734,0.993521,0.994235,0.994883,
0.995470,0.995999,0.996477,0.996906,0.997291,0.997636,0.997944,0.998218,
0.998462,0.998677,0.998867,0.999034,0.999181,0.999309,0.999420,0.999516,
0.999599,0.999669,0.999730,0.999781,0.999824,0.999860,0.999889,0.999914,
0.999934,0.999950,0.999962,0.999972,0.999980,0.999986,0.999990,0.999994,
0.999996,0.999997,0.999998,0.999999,1.000000,1.000000,1.000000,1.000000,
1.000000,1.000000,1.000000,1.000000]
yl = [0.0 ,
0.24284189089843749999999999999999999999999999999999999999999999999999999999999999999999999994 ,
0.33891381153616241091128013830897494597501346175061815085181106016791140839106123644489567939 ,
0.4096508118986554964468613586176876286675910249176568910063898986569389242136034043611362893 ,
0.46686464000000000000000000000000000000000000000000000000000000000000000000000000000000000004 ,
0.51520878690160297251424266470557800821360333656077840659182490399871624315849062399794482769 ,
0.55711036564679084226463064132957720619760515698571085504807992813841969684759803375721161952 ,
0.59403586694497658800362008610502523324778590543874537390636537671367905563905898381357396092 ,
0.62695727405708851860034476690992031585847792570867328348262569558774758915624766495755680337 ,
0.65656142894531250000000000000000000000000000000000000000000000000000000000000000000000000042 ,
0.68335708497998774418637183634179332886070075779714266606849987777298155922511603856062010212 ,
0.70773493804111727989792316835045512711884425667708681706322072645727149531959242943004327479 ,
0.73000371263631036412805879027038237381107153386411289288452436708653059246123542862051058887 ,
0.75041306705768700279223654096504325232259994688105736914495222220870131424495073764190801567 ,
0.76916878371148342382905920508248142729758643104785177865664504205039872657201495788502307921 ,
0.7864432137848497988936034369912098416012138396347553611156269933212088243478330208682259075 ,
0.80238267999999999999999999999999999999999999999999999999999999999999999999999999999999999997 ,
0.81711286071375848422156171309553776882662355675402295720883913794315191359386876273687027721 ,
0.83074279409297163612989839239901355790328642749759040271797040235142889794661677568874910432 ,
0.84336791448053669039491494153704053659820418852533323532676626011267281281150806299497728187 ,
0.85507239459591957990686801092284003243249246071544743696119110664498320196389659353607207076 ,
0.86593097989850596594837564628024840742735289184755416868082291691632424902641029664473753592 ,
0.8760104448405395845869865032923832133406950636382132179935302466877259466115059757797386568 ,
0.88537076313537333592619870160765776265166719730652528570716786016906623561279315729986731746 ,
0.89406605863621822202942099271468058160547927843031641414785351377761375367467479456317364216 ,
0.90214538574218749999999999999999999999999999999999999999999999999999999999999999999999999902 ,
0.90965337578828192224035974750769403151427194445248985922715562926098364967600669364954243689 ,
0.91663077695127165195449156354560512107971566758366948486365356259181728453941365487697389883 ,
0.92311490871651899672635366394183484869868472458007637049934880188714360538438574624163388376 ,
0.92914004717195174562814400085624770217681110300114980976688722575217657102992330766207636557 ,
0.93473775383109182738232483397849235115708113892024964839619232601671793110279283985177984147 ,
0.93993715799810172595930388260563714156050126731320344011491712304666760352103323523009109971 ,
0.94476520063771866258056229600564189622469986798494511505191211776314027334871319185115673214 ,
0.94924684613452738234910366497259242496729386538037598394930269126512715830059565398768019205 ,
0.95340526709961673348471574282263489598370997623223600736415494605626469268010581071048617756 ,
0.95726200642161857863846491592899936722162878090139612998779028899532846996400976088964743294 ,
0.96083712000000000000000000000000000000000000000000000000000000000000000000000000000000000119 ,
0.96414930299435167900788691619332117832203512400168261060341487065187992434392786970547752749 ,
0.96721600193925302030934544260076416867109738689998269907065034689700396609076571546028167964 ,
0.97005351468370397389109612097610181624137119694401812560990641803403820584993954947060297708 ,
0.97267707979704137778120966809819881199046366641719513073103277390834214189494668640772391681 ,
0.97510095682434692257637936255653861178871489036595619299674215292371523713278853846663881901 ,
0.97733849856175537024158077522325560964487104131309408816396649311589133447065168996771979626 ,
0.97940221634657789004627294084038683837159432182292799090904360632634034008475795920738963008 ,
0.9813038392115632602529953961135595528850475388389553476210038153837838451036224652716555394 ,
0.9830543676312583121652840547340179758282166241842901521070540913510103651881511682571774964 ,
0.98466412248680660481695164755205342445621773349724104142564909738435646555144635042994768258 ,
0.9861427897900626043509871861466734729434100517557261196279812424541334226121732414696383538 ,
0.98749946163573132350086704721747485065199174319509404344579164508042728288390624424518808893 ,
0.98874267378906249999999999999999999999999999999999999999999999999999999999999999999999999289 ,
0.98988044026456628537950233698955185333477765300432376118104023652818115618917501962501490781 ,
0.99092028520675547756927657869767514625497471254021485454478628559282113802813701799869061029 ,
0.99186927234581715315020389934364678115383828103382096069543059192093474072734364014398248663 ,
0.99273403226835652859588682011216423141238547274935197462333435393415051148821239674768469124 ,
0.99352078771509967766904028098489835265409751044197015622200010393766631479073091943416031178 ,
0.99423537709299594857433268082780666907822943008687584134452568353339088659461681274309139503 ,
0.99488327636794980955364503365590615773829667420436680019148205831377192464193660610071859717 ,
0.99546961948595557538419843474362356824096930343605276306400563004120603535053652346933378663 ,
0.99599921745430571405251231145256860061404212262248291641902847616628477980515490158256556053 ,
0.99647657620045740279705919932863909220955855139220838149475449307751848872585051973063705386 ,
0.99690591331378910624514291390398965982839404694205547876364148739761054786457840322786454437 ,
0.9972911737646196406748698135100625080127368540864953160301387919931645120677780588908225332 ,
0.99763604468529350381466510728306222868014164715618175753152924261043462553842136687017176744 ,
0.99794396928968571140813956219204831945417449876828093682576058041835955268192625333260277607 ,
0.99821815999999999999999999999999999999999999999999999999999999999999999999999999999999997878 ,
0.99846161084310041500069686286281088176356206094864289478995096615288641305409189906448282553 ,
0.99867710917272036884717344321611779205107785241151850615413447918452206921703107191605120959 ,
0.99886724676864275866760426778805652431211058502205780011048754780735697411135842013423038186 ,
0.99903443035926004065291935093271146586408258427491276926821444879325269485186443830143429494 ,
0.99918089160973551859020962271424090694662219601422323338046800613082795098543718783562691518 ,
0.99930869661423702374469424143851337136395694150135586921745548653543127671815050447241678735 ,
0.9994197549273500435075051741724820162796637716652601281502408851559135726469356252886546668 ,
0.99951582816675434577119722618660722087540857259110098630883361254677200551543357329101176599 ,
0.99959853821652717554576793253894404826832644852415634614597491199695067902809375111917474933 ,
0.99966937505798306635275686224574589724951142520477861690150271287938807454938431579263999228 ,
0.99972970425274538241978016296604360999062338769356182887980668047681076401686369327802246513 ,
0.99978077410074174720519835700003364388884892559266683522562601500038765678854176370608271056 ,
0.99982372249400156589974170775368597962625547827351114200523679295335949665671749759653276069 ,
0.9998595834854887227365058419886525983880825432513592177322972465581824201970105375355119382 ,
0.99988929359070843416463989823946904242277023453716937447736095377775350813604694867820523722 ,
0.99991369783846845749177223531490845058315851495928401575083847571644969492208540589564099485 ,
0.99993355558593750000000000000000000000000000000000000000000000000000000000000000000000122951 ,
0.99994954611201544088640575321738489907966380056919005255891523802409670658280248445431341549 ,
0.99996227400199995046511337782573540008555582104927115215602031585422082441605359732093502956 ,
0.99997227433559256658318990015022519191691800101664040193622710118322574777668209306437116156 ,
0.99998001768942563303514828689088446610505485130106710346151407012337075461934256860728117613 ,
0.99998591496450202305723970803777334137528667202800382205467112289147305605797793252820125127 ,
0.99999032204821539313861238459477763315042258387828404717750393297595119785958622465047279352 ,
0.99999354431995369756988632689180527167366509671666910940221426340158158046673309148178192454 ,
0.99999584100867734532735333587968385321923654026367972456361071170090755329403470576479028698 ,
0.99999742941030077062723292142719038847885249493924818563712942813850378689470607380670514651 ,
0.99999848897218789404153066161828214238663254275651009559999867170499953781103513600329377805 ,
0.9999991652515939975493381838959553590807341713452876835440809011642513768697455834881878792 ,
0.99999957375444534736372202112054442747049320616469956985352583501654882933367159670192781897 ,
0.99999980366044024999801358769995743352765561964288093009712980354806082359478156170081000405 ,
0.99999992144007821252063221435914256823196175173011656922068735757539335159080466598914197555 ,
0.9999999743688748715227123241700716666369038775484188504600258363489889024450996330027423859 ,
0.99999999394369698272337153966612807655364499020690284046106281730497065553706150916743739046 ,
0.99999999920585187507528937783928310574064898456039996036175892635729534893096765916598457299 ,
0.99999999997528742191360455635153690385626066481481896129682592406691532660547012408294797547 ,
1.0]
#plt.plot(xl,yl,'-',xl,Yyl,'--')
plt.plot(xl,yl)
plt.ylabel("Ix(0.5,5)")
plt.xlabel("x")
plt.grid(True)
plt.savefig('../img/aptC.pdf')