Untersuchung über die Bestrahlung auf irgendeinerschiefen Ebene für die Zeitdauer länger als ein Tag.

Sci. Bull. Fac. Ednc., Nagasaki Univ., No.27 pp.197-215 (1976) 197

Untersuchung über die Bestrahlung auf irgendeiner schiefen Ebene für die Zeitdauer länger als ein Tag.

Takao Sato

Fakültät von Erziehung, Nagasaki Universität

Ubersicht

Für die Berechnung des Sonnenscheindauers beziehungsweise der Bestrahlung auf einer beliebig geneigten Ebene für die Zeitdauer langer als ein Tag muß man vier fundamentale Formeln(5), (6), (7), (8)zu vier leicht integrierbaren Formeln bezogen auf die Sonnenlange verwandeln.

1 Problem-Stellung

Es sei nunw die mittlere Bestrahlungs-Starke für eine Minute auf irgendeinem Breitenkreis. Dann ist der Wert für einen Tag auf irgendeinem Punkt auf jenem Breitenkreis.

W W T (1)

Deshalb wird der Wert fürm Tage

Wm=-. (2)

,hierbei t1 der Zeitpunkt der unteren Kulmination der Sonne in erstem Tag ist.

Bezeichnet man sogar der Abstand zwischen Erde und Sonne in einer astronomischen Einheit mit p, die Sonnenlange mit x,, die Exzentrizitat mit P und ein Jahr mit T, dann ist

offenbar

(3) Folglich hat man

Wm=

,hier ist{ }irgendeins von vier fundamentalen Formeln, sind. Das heißt :

die in Reference 1 angegeben

198 Takao Sato

2 sinasin8 • Ih +2 cosa cos8 sinro (5)

(to +11.0 —A)sina sin 8+cos a cos8 • (sin (to —A) +siniko) (6) 2 (lr 0 ± t O - 7r) sina sinb+ 2 (sin Ito +sin to cosA)cosa cos 8 (7)

2 to sina sina+ 2 cosa cos8 sin to cos A (8)

Zur Umreohnung von diesen Formeln muß man derartig so verwandeln, daß man diese leicht integrieren kann, bezogen auf X.

'Nun setzen wir cos t

o = —tgg) tg8, cos —tga tg8 (9)

und weiter n = tg k= sin E, p = tga (10)

Dann wird

und daraus in 8=k sin, (11)

cos to = —ntga, cos = —p tg8 (12)

sin'tko =1/1— (1 + p 2)sin 2 8/cos 8 (13)

cos Uhilka =1/1— (1 +p2)sin28 (14)

Die letztere Formel laßt sich entwickeln, wenn (1 +p2)sin28 <1

Da aber (1+p2)sin28=sec2asin28 (15)

,so ist die Bedingung der Entwicklung

— (16)

Der Erfolg wird nun mit

111 cos 8 sin*

0 = 1 —2—(1+p2)k2sin2X--(1+p2)2k4sin4X--(1+p2) 2324 3k6sin6x

—1.3.5.2 m • (2m-2)k2rnsin2mX(1+132)m— M! (17)

Nun folgt1(—ptg8) =-2--+ sin- (p.tg8) (18)

Es laßt sich entwickeln solange-1<ptg8<lfolgendermassen.

..—-1-ptg8+-2--•tg8+-j3 113313„g •t8+

(2m) v.+ 1 tg62m+ 1

22m• (m!) 2 • (2m±1)(19)

Da aber

tg8 = sin 8

V—cos28=ksinX (1 — k2sin2x,)-+(20) ,so folgt

osin8 = 2—7rksinX+ pk 2sin 2 x(1— k 2 sin 2 X) -+

Untersuchung Uber die Bestrahiung auf ivgendeiner schiefen Ebene für die 199 zeitdauer länger als ein tag.

1 •1 -1)3 k4 sin4X (1 —k2sin2X) 3 +

(2m-2)! P2m-1 2m-1

+

22m-2 ((rn_i)/) 2 2m —1k2ln sinuax (1 —k2Sin2X) 2

+22m (2 •(1110 2 2m+1m)!p2m+12m+1 k2.+2sinz.+2x(i_kzsinzxy2

(21)

P sin8 + cos 8 sin o =-7`2pksinX+E am k2m Sin2mX (22

,hierin

ao=1

1 a l---2(1+P2)+P2

111

a2-_„, 22 (1 +p2)2+n2+1-14 .2!\'2"23'-

1.3 (,P2

)31.31 3_4j_1•3•1,

a3—22 .3!"-22.2r2.3.21'2.4. 5P

1.3. 51.3. 513.5

a 4=-24

.4!(1+p2) 4+22.3!P2+2.3•22.2!p4

1.3 • 151.3.51 .n6+ . (23)

2.4 5 2'2.4.6 7'

Weiter das allgemeine Glied zu ausdrücken ist das folgende Verfahren empfehlich.

Das erste Glied ist naturlich der Koeffizient vonk2m sin2mx, welches aus (17) sich entsteht. Das zweite und weitere stammen aus (21).

Nun fur Kurz setzen wir

(24) ,so folgt bei Km

1 3

pe.o sin3'2't2 71. (1 _12)1—32+•____1n4 t4 (1___12) +

21-1

!p21—2 • t21 • (1 —t2) +

221-2 • Ge 1).0 2 22-1

21+1

2/ (22 )! P21±2 (2E+1) t2.1+(21 .—t2)- +

2•(V) 2

200 . Takao -Sato

_2m+1

+2'(m!)

(2m)!2P2m+2(2m+1) •2m+2

^••

• t • (1 —42)(25) 2 Da am der Koeffizient von tm ist, so wird das letzte Glied von am

(2m-2)! P2m

22m-2 ((m— 1).02 • 2m-1(26)

Zur Ausdrucken des allgemein enGliedes, ist es. nur nötig den Koeffizient von t2m-1 aus der Serie

(1-12)-2111 •

, die in dem t2' einschließenden Glied sich befindet, zu suchen. Das heißt (1 —t2)-n—l+nt2+ (-11)(-n-1) 2!! (-t2)2+(-n)(-n-31)(-n-2) ( t2) 3

+ (- - 1) (q!n-2)--(-n-q+1)( t2)q

=E ( -n) ( n 1) ( n q+1)( t2)q (27) q=0Cl!

, wobei setzen wir nun

q=m—e, n=22-1 2

und sogar mußen wir jeden Koeffizient in dem t2' einschließenden Glied, aussch- ließlich (1- t2) ,multiplizieren. Folglich

(22 -2)! P2/ t2t .t2m-2t(-1)m-1 (22-1\( 22-1..)(22-1

221-2 (U— 1).02 22-1(m-2)!• \ 2I2 -I- 2

(22

2-1m+2+1) (22-2)!

1-.P2' • t2m ( —1) 2m-21(2i-1(22+4(22+3...(22-3+2(m-2)) 222 • ((i —1) .02•2.2- 1 (m 2) ! \ 2 )\ 2I\ 2 /2

(22-2)!p21t2m 1 (2e-1)(22+1)(22+3)-•(2m-3) —

2"-2 • (i —1) !(Q-1) ! • 2i - 1 • (m- ! •

(22-2) ! 1 1 (22-1).22.(22+1)(22+2)•-(2m-3)(2m-2)

-

2m+1-2 • (22-1) • (m-2)! • ((2-1)!) 2 • 22(22+2) (2m-2)

p2t t2m (2m-2) !

• (22-1) • (m-2)! - 1) ! 2 • 2's • 2(2+1) ( .2- 2) (m -1)

1 1 1 1 (2m-2)! ,21tym

22m-2 • 22-1 • (m-2)! • (.2- 1) ! •Cm-1)! • •(28)

Das ist denn das allgemeine Glied von am. Das letzte Glied von am (26) laßt sich

ausdrucken, indem man in dem Koeffizient vont2m, i-m setzt, und sogar der aus (17)

sich ergebende Koeffizient laßt sich folgendermassen umschreiben.

Untersuchung Uber die Bestrahiung auf ivgendeiner schiefen Ebene für die

zeitdauer länger als ein tag. 201

_1.3.5. (2m-3)2 (2m-3) (2m-2)

2m • m!+P"--2m"

(2m - 2)!(2m -2)!

2m • m! • 2m-1 • (m-1)!(14-P2)m-22'1 • m I (m-1)+P" (29)

Dann ist der Ausdruck von am folgendermassen geschrieben (2m-2)

-

22,1m (m-1) ! (1 +p2)m+

mz1 1111 (2m -2)!

,24_ (2m -2)! p2m

t ---122m-2•22-1•(M-0!(E-1)!•(m--1)!22m-2•[(M-1).1)2•2m-1

(2m-2)! 2\ •(2m-2)!m11 1 p2/

--

22m-1.m./- 1)!"m"P22m-2h2E-1•(m-4!•-!• (rn-! (3°)

ao

--1

2-(1+p2)+p2

2! a 2==23

.2.1.1!(1+p2)2+22p2+222!.3p4

(.14_ 4! 7,2+24 4! r,4, 4! n

a3- 25

.3! 2!s-+P2-324.2! 2!.3.2.1'424.5.2! 2.16'- 6! .4._6!6! a4- -

274!3"(1P2»26.1.3! 0! 3!/32+ 26.3.2! 1! 3r426.5.1! 2! 3!P6+ 26.7.0! 3! 3!138

82,88! ( p2 p4 + p6 7,8+

4!p1 °/

a5=29

.5! 4.1(1+p)+28.4!1.4! 0!3.3! 1!5.2! 2!7.1! 3r9.0! 4!

10! (14 ..,2) 6+ 10! ( 1 n2+ n4+ p6

a6=-211

.6! 5!21°.5!1.510r3.4! 1!"5.3! 2!

17,84 _1 ni 04_ 1

7.2! 3r9.1! 4!'11.0! 5!"/

12!^12! ( 1 T,2 _4_ 1 "44 ps

a7et 2"6!1.6! . OP'3.5!1."5.4! 2!

17„84 _ 1 niO4_ n1.2.4_ 1 n14)

+7

.3! 3IF9.2! 4!"11.1! 5!"13.0! 6!'

14! (14_p2)8+ 14! ( 1 ,24_ rd+ p6+ p8

a8,_. 215

.8!7.1214.7!k1.7! Or3.6!5.5! 2!7.4! 3!

1 ni 0 1 n124 . 1 nI,44 1 n16)

+9

.3! 4!"11.2! 5!"13.1! 6!'15.0! 7!"

202 Täkao Sato 16!

ao--(1+1)2)9+21.66 2179!8!1.8! .( 1 p2+ 0!3.7! 1r,4+5.6! 2!'7.5! n6+ n 3!8'

1"io4. 1 ,124. 1 ,04.4_ 1 „16.4_ 1 nis)

+9 .4! 4!' 11.3! 51' 13.2! 6!'15.1! 7!'17.0! 8!'

18!1! (1 n , ,

ai0=

219.10! 9/(14_p2)10+ 218.89!\1.9! Or2+3! lr4+5.7! 2!'8+7.6! 3r

1 n10, 1 n12 _1_ 1 fo.4 1 n16 +

9.5! 4rs11.4! 5!'13.3! 61'15.2! 7!'

1ol8+ 1

+17.1! 81-19.0! 9!P20) (31)

Nun in (22) laßt sich sin2mx folgendermassen umschreiben :

1 m

sin2mX= 2m)+2 .-1m—E( —1)s (2ms)cos2sx} (32)

, so folgt

•

cosa(plko sin8+cos 8sin 1,>0)x+ 2 E a. k2msin2m m-o

_1(2mam k2m

=COS a(lr pk sinX+

2 m-02'\ m

m

+ E E am k2m(—1),s (2m)cos 2sk.) (33)

„0m_s

Wenn man diesen Ausdruck in der Form

a' sinX+ E (-1)s bs cos2 sx (34 umschreiben, hat man

a' =-!--rsina sin 6 2

bcosa(2m o = E)k2m

am .-o2m

b. = E cosa (2m)k2m

.-022n1-1\111—s

1 - cosat2r+ 2s\ )a r+sk2r+2s (35) 22s-1

r.--0 22r k r

Damit ist die fur Integration bequeme Form gemacht worden.

Es handelt sich noch weiter um die Verwandeung von

cos 8 sin to =-1/1— (l+n2)k2 sin2X. (36)

Untersuchung über die Bestrahlung auf irgendeiner schiefen Ebene für die Zeitdauer länger als ein Tag.

Diese Form laßt sich fur alle X entwickeln, wenn sin2 8 • (1+n2)< 1 d. h.

( 7r 8) <9<--8

(37)

(38)

203

2. Das Fortran-Programm

Der Verfasser hat auf Grund dem bisherigen Erfolg das folgende Fortran -Programm geschafft.

(Unvollendet Fortsetzung fo;gt)

Reference 1 Berechnungsmethoden der Bestrahlung und Sonnenscheindauer auf der schiefen Ebene. Science Bulletin of the Faculty of Education Nagasaki

University No. 25

204 Takao Sato

9 6

87

94

93

910

97 98 9

101 500

102

REAL (1-).P(11)‹t511,51),B8(11951)93(11,51),G(11951),R(26.2 10),F(51),J-C11951)914(11751) aZ(11)981(11225),E(51),D(50,50)

--73.11.1593/280.0 EK=SPV23445*Q) 00 1 M=1,40 E(1)=1.0 FM=FLOAT(M

F(M41)=E(M)*(2.0e=M-1.1')/2.0 CONT1N.E

WRITE(be 96) (E(M),1=1.50) F('RMAT'IP '10E1245)

F(1)=2.':

DO 87 V.,9'1 FM=FLrl-(M'

F(M+1)=J="2.03:F4+?. 1.+7403eFM)/((1.0+FM)X2) CONT14 E

tAjTfi-,n1')

WRITEtte94) (m,M=111 ) rnelA7(1H '10112) WP1TE(''.600)

WRITE(6993) (F(4),"1=1,50) P0R1!A'11H 910E124‘) WRITE(e9600)

00. 9 i=J925 R(11/)71.0 FI=FIL,T(I)

DO :=1.25 FL=FLOAT'L'

R(L+1,1)=R(L,I)2.0e(F,_+-1,(2.0FL+2.v»F1-1.0)/(FU“2.0.xFj+FL)) CONTIF

WRITE(69 WR1TEt6, FoRMAT(IH F3RMATtiH CONTINUE DO 100 J=3

Q7) 1,(R(L,i',L=1,15) (3) 1,(( _1,121-'214,25)

,1391X,1HR91Xe13E9,4) 1391X,14P,1X•12E9.4) DO 100 J=3e11

AL(J)=-90+15*(J-1) P(J)=TAKCAL(J)341) C(171)=P(J)**2

WRITE(69101) AL(J) FORMAT(11-15,2HALIF10.0) WR1TE(6,600)

FORMAT(IHO) WRITE(6,600) 00 2 11=1,50 C(1,M+1)=C(19M)/FLOAT(M) WR1TE(6,102) M

FORMAT(1H,I6)

WRITEt6,600)

00 3 L=1,M

FL=FLCAT(L)

3

103

4

2

5

11 10

100

9]

90 92

89

Untersuchung über die Bestrahlung auf irgendeiner schiefen Ebene für die

Zeitdauer länger als ein Tag. 205

OCL,m)=(2.::**L-1.0)*rLUAlcm-Lietp(J)**2/((2.0*FL+1.0)*L) C(L+1,m)=C(L,M)*0(L,M)

CONTINUE

WR1TE(6,103) (C(L,m),,=1,M) WRITE(6,1n3) (D(‘9M),L=1,M) FORMAT(1H 910E12.4)

5=0.0

DO 4 L=1,M S=S+C•L,M) CONTINUE 813(J,M)=5 B(JIM)=E(m)*8B(J,M) FM=FLOAT(m)

G(.1,1)=(1.0+P(J)>*2)/2.0

G(J,M+1)=G(J,M)g(2.0FM-1.0)e(1.0+P(J)>J*2)/(2.0*(1.c+FM)) A(J,M)=-G(J,M)+8(J,m)

CONTINLIE 5=0.0

DO 5 m=1,5'

H(J,m)=F(m)gA(J,M)EK—(224)/(2.0;ak(20)) S=S+H(J,*°)

cONTIUF

BZ(J)=(1.V+5>COS(4L(J)) DO 10 1=1,25

s=0.0

DO 11 L=1,25

S=S+R(L,I)i:A(J,L+1-1)*EK**(2n)/(2.0e*(2*(L-1))) CONTINUE

BI(J,1)=53;EKeae(2*1-2)>ZOS'AL(J)*2)/(2.0g02*1-1)) CONTTuE

CONTImuE WRITE(6,600)

WRITE(6,91) (BZ(/),1=1,13) FÜRMAT(1H ,13E10.4)

DO 92 1=1,25

WRITE(6,90) (BI(JII),J=1,13) FORMAT:1H '13E10.4)

CONTANUE WRITE(6,600) DO 89 M=1,50

WRITE(6,90) (B(J9M),J=1.13) WRITE(6,90) (G(J9M)9J=1,13) WRITE(6990) (A(J9M),J=1›13)

WRITE(69600)

CONTINUE

STOP

END

206 Takao Sato

1

INTEGER AL(11)

REAL AA1(11).AA2(11),AA3(11).AA4(11).AA5(11).AA6(11).AA7(11) AA8(11),AA9(11).AA10(11)

D10=-0.5 A1=1.0 D11=A1

020=-44125 A2=0.5 D21=1.0*A2 D22=A2/3,0 A3=6.0/2.0**4 030=-1,0/6.0 031=A3 D32=A3*2.0/3.0

D33=A3*0.2 A4=120.0/2,0**6 C40=—A4/48.0

D41=A4/6.0

D42=A4/6.0 D43=A4*0.1 C44=A4/42.0 A5=30,0*56,0/2.0**8 D50=•..A5/240.0

D51=A5/24.0 052=A5/18.0 053=A5/20.0 054=A5/42.0 D55=A5/(9.0*24,0) A6=540.0*56.0/2.0**10

D60=—A6/1440,0 D61=A6/120.0 D62=A6/72.0 D63=A6/60.0 064=A6/84.0 D65=A6/(9,0*24,0)

D66=A6/(11,0*120,0)

A7=12,0*11.0*10,0*9,0*8.0*7.0/2**12 D70=—A7 /(2.0*5040.0)

071=A7/720,0 D72=A7/ 360.0

D73=A7/240.0

074=A7/252.0 D75=A7/(18.0*24.0)

D76=A7/(11.0*120.0) C77=A7/(13.0*720.0)

A8=8,0*9.0*10.0*11.0*12,0*13.0*14.0/2.0**14 D80=•••A8/(2.0*40320,0)

081=A8/504090

D82=A8/(3,0*720.0)

Untersuchung über die Bestrahlung auf irgendeiner schiefen Ebene für die

Zeitdauer länger als ein Tag 207

0049 0050 0051.

0052 0053 0054 0055 0056 0057 0058 0059 0060 0061 0062 0063 0064 0065 0066 0067 0068 0069 0070 0071.

0072 0073 0074 0075 0076 0077 0078 0079 0080 0081 0082 0083 0084 0085 0086 0087 0088 0089 0090 0091 0092 0093 0094 0095 0096 0097 0098

D83=A8/(5.0*240.0) C84=A8/(7.0*24.0*6.0) D85=A8/(9.0*24.0*6.0) D86=A8/(11.0*240.0) C87=A8/(13.0*720.0) D88=A8/(15.0*5040.0)

A9=9.0*10.0*11.0*12.0*13.0*14.0*15.0*16.0/2.0**16 D90=—A9/(2.0*362880.0)

C91=A9/40320.0 D92=A9/(3.0*5040.0) D93=A9/(5.0*1440.0)

D94=A9/(7.0*6.0*120.0) D95=A9/(9.0*24.0**2) C96=A9/(11.0*6.0*120.0) 097=A91(13.0*2.0*720.0) D98=A9/(15.0*5040.0)

C99=A9/(17.0*40320.0)

A10=11.0*12.0*13.0*14.0*15.0*16.0*17.0*18.0/ 2.0**18 D100=—A10/(2.0*3628800.0)

D101=A10/362880.0 E102=A10/(3.0*40320.0) D103=A10/50400.0 D104=A10/(7.0*6.0*720.0) D105=A10/(9.0*120.0*24.0) D106=A10/(11.0*24.0*120.0) 0107=A10/(13.0*6.0*720.0) D108=A10/(30.0*5040.0)

D109=A10/(17.0*40320.0) D1010=410/(19.0*362880.0)

hRITE(6,2)

R1TE(6*/) D10,D119020 tD21 eD229D309D31 t03290339040 WRITE(6,3)

WRITE(6,1) 0419D422D439044,050,D51205290531D549D55 WRITE(6t4)

WRITE(6,1) D609D6190622D63,D649065906690709D719D72 hRITE(695)

WRITE(621) D739D74,D759D762D77,D80,D819D829D83,D84 WRITE(6t6)

WRITE(6,1) 0859D869D879088,0909D91,D92,D93,094,D95 NRITE(697)

WRITE(6,1) D96,D97,D98,D99,D100,D101,D102,D103,D104,D105 kRITE(6v8)

WRITE(691) D106,D107,D108,D109,D1010

WRITE(6, 9) A1tA2,A39A4A5,A6,A7,A87A9,A10 C=3.141593/180.0

DC 100 J=1,11 NRITE(6,600) hRITE(6,600) AL(J)=.-90^15*J

INRITE(6,101) AL(J),AL(J),AL(J),AL(J),AL(J)

208

0099 0100 0101 0102 0103 0104 0105 0106 0107 0108 0109 0110 0111 0112 0113 0114 0115 0116 0117 0118 0119 0120 0121 0122 0123 0124 0125 0126 0127 0128 01.29 0130 0131 0132 0133 0134 0135 0136 0137 0138 0139 0140 0141 0142 0143 0144 0145 0146 0147 0148

Takao Sato

P=7AMAL(J)*0) PP=1•0+P**2

El0=DIO*PP El1=D11*P**2 AA). (J) =E10+Ell E20=D20*PP**2 E21=D21*P**2

E22=D22*P**4

WRITE(6210) AL(J) pE109EllsAAI(J) AA2 =E20+E21+E22

WRITE(6,11) AL(J),E20,E219E22,AA2 (J) E304)30*PP**3

E31=D31*P**2 E32=D32*P**4

E33=D33*P**6 AA3(J)=E30+E31+E32+E33

WRITE(6,12) AL (J) 2E309E312E32 sE339AA3(J) E40=D40*PP**4

E41=D41*P**2

E42=D42*P**4 E43=D43*P**6

E44=D44*P**8

AA4(.1)=E40+E41+E42+E43+E44

hRITE(6,13) AL(J) 2E409E419E422E43,E44,AA4 (J) E 50=D50*PP**5

E51=D51*P**2 E52=D52*P**4 E53=D53iiP**6 E54=D54*P**8 E55=D55*P**10

AA5 (J)=E50+E51+E52+E53+E54+E55

INPITE(6214) AL(J),E509E511E529E539E549E559AA5 (J) E60=D60*PP**6

E61=D61*P**2 E62=D62*P**4 E63=D63*P**6 E64=D64*P**8 E65=D65*P**10 E66=D66*P**12

AA6 (J)=E60+E61+E62+E63+E64+E65+E66

WRITE(6,15) AL(J),E609E612E629E639E649E659E669AA6 (J)

ET0=D70*PP**7 E71=071*P**2 E72=1)72*P**4

E73=D73*P**6 E74=D74*P**8 E75=D75*P**10

E76=076*P**12 E77=D77*P**14

AA7(J)=E70+E71+E72+E73+E74+E75+E76•E77

0149 0150 0151 0152 0153 0154 0155 0156 0157 0158 0159 0160 0161 0162 0163 0164 0165 0166 0167 0168 0169 0170 0171 0172 0173 0174 0175 0176 0177 0178 0179 0180 0181 0182 0183 0184 0185 0186 0187 0188 0189 0190

0191

0192

Untersuchung über die Bestrahlung auf irgendeiner schiefen Ebene für die Zeitdauer länger als ein Tag

20

100 2

3

4

WRITE(6,16) AL(J),E70,E71,E72,E73,E749E75,E76,E77,AA7(J)

E80=D80*PP**8 E81=D81*P**2 E82=D82*P**4 E83=D83*P**6 E84=D84*P**8 E85.D85*P**10

E86=D86*P**12 E87.1387*P**14 E88.D88*P**16

Ate(J)=E804£814.E824.E83+E844.E85+E86+E874-E88

WRITE(6,17) E801E819E82,E83,E84,E85,E86,E87,E889AA8(J) E90=090*PP**9

E91=D91*P**2 E92=D92*P**4 E93=D93*P**6

E94=094*P**8 E95=0)95*P**10 E96=D96*P**12 E97=D97*(2**14 E98=D98*P**16 E99=D99*P**18

AA9(J)=E904.E91+E924,E93+E94+E95+E96+E97+E984.E99

WRITE(6,18) E90,E91,E92,E939E94,E95,E969E97,E98,E999AA9(J)

E100=D100*PP**10 E101=D101*P**2 E102=D102*P**4 E103=D103*P**6 E104=D104*P**8 E105=0105*P**10 E106=D106*P**12

E107=D107*P**14 E108=DI08*P**16 E109=D109*P**18 E1010=D1010*P**20

AA10(J).E1004.E101+E102+E103+E1044.E105.E106+E107.E1084.E109+E1010 WRITE(7,20) AA1(J),AA2(J),AA3(J),AA4(J),AA5(J),AA6(J),AA7(J), IAA8(J).AA9(J),AA10(J)

FCRMAT(5E16.8)

WRITE(6,19) E100,E101,E102,E103,E104,E105,E1069E107,E1089E109,E10 110,AA10(J)

CONTINUE FORMAT(1H0910E12.5) 1

2 3 4

FCRMAT(IHOt12H 021912H

D32,12H FCRMAT(IH0,12H

D44,12H 053,12H FORMAT41H0,12H

D10,12H D22,12H D33,12H

041,12H 056,12H 054,12H

060912H

Dllt12H 030,12H D40)

042,12H 1 D51,12H D55)

061,12H

020,12H D31t12H

D43112H D52,12H

062,12H

209

0193

0194

0195

0196 0197 0198 0199 0200 0201 0202 0203 0204 0205 0206 0207 0208 0209 0210 0211

5

6

7

8

9 10 11 12 13 14 15 16 17 18 19 101 600

5 6 7 8

9 A B c

O63,12H 070,12H FORMAT(1H0d2H

076,12H 082,12H FORMAT(11101,12H

O88,12H 093,12H FCRMAT11H0,12H

099.12H 0103.12H FORMAT(1H0,12H

0109.12H

064,12H 071,12H

073,12H 077,12H 083,12H

085,12H 090.12H 094,12H

096,12H 0100,12H 0104,12H

0106,12H 01010 ) fORMAT(10E13.6)

FCRMATI1H0, 4HALPH,15‚ 3E14.7) FORMAT(1H0,4HALPH115.4E14.7) FORMAT4IH0.4HALPH45.5E14.7) FCRMATt1H0,4HALPH.15.6E14.7) FORMAT(1H0,4HALPH,I5,7E14.7) FORMAT(1H0,4HALPH,15,8E14.7) FCRMAT11H0,4HALPH,15,9E13.6)

FORMAT(1-10. 10E13.6)

FCRMAT(IHO, 11E12.5) FORMAT(1H0.12E11.4) FORMAT(1H0.516)

FORMAT(1H0) STCP [ND

1 065,12H D72)

074,12H 080,12H D84)

086,12H 091,12H D95)

097.12H 0101.12H 0105)

0107,12H

066112H 075,12H 081,12H

087,12H 092,12H

098,12H 0102,12H

0108,12H

ts..., i--, 0

A) 0

D10 -.50000E+00

D41 .31250E+00

Dll D20

.10000E+01 -.12500E+00 D42 .31250E+00

D43 .18750E+00

021 .50000E+00

D44

D22 030

.16667E+00 -.16667E+00 D50 .44643F-01 -.27344E-01

1 051 .27344E+00

D31 .37500E+00

D52 .36458E+00

032 .25000E+00

D53 .32813E+00

D33 D40

.75000E-01 -.39063E-01 D54 .15625E+00

D55 .30382E-01

C/") D ,-t- a 0

D60 061 D62 D63 064 1 D65 D66 070 D71 D72

-.20508E-01 .24609E+00 .41016E+00 .49219E+00 .35156E+00 .13672E+00 .22372E-01 -.16113E-01 .22559E+00 .45117E+00

D73 .67676E+00

D85 .81462E+00

v96 .99976E+00

D106 .21245E+00 .100000E+01

D74 .64453E+00

086 .39990E+00

097

.42297E+00 D107

D75 076

.37598E+00 .12305E+00

087 088

D77 D80

.17353E-01 -.13092E-01 D90 .11279E+00 .13965E-01 -.10910E-01

Des D,9 01,J

.10474E+00 0108

.11552E-01 -.92735E-03

D109 01010

091 .19638E+00

0101

.18547E-01

.11984E+00 .44513E-01 .98190E-02 .97616E-03

.500000E+00 .375000E+00 .187500E+01 .656250E+01

081 .20947E+00

092 .52368E+00

U1U2 .55641E-,01

.295313E+02

D82 .48877E+00

D93 .10997E+01

0103 .13354E+00

.162422E+03

D83

.87979E+00 D94 .15710E+01

0104 .22256E+00

.105574E+04

D84 .10474E+01

095 .15274E+01

D105 .25966E+00

.791807E+04 .673036E+04

-75 ALPH ALPH ALPH ALPH ALPH ALPH ALPH

-75 -75 -75 -75 -75 -75 -75

-75 -75 -75

-.7464110E+01 -.2785647E+02 -.5544633E+03 -.1939957E+04 -.2027207E+05 -.2269694E+06 -.266220E+07

-75 .1392822E+02 .6964110E+01 .5223082E+01 .4352569E+01 .3808497E+01 .3427648E+01

LPH -75 -.266220E+07 .314201E+01 .875252E+02 -.322903E+08 .291758E+01 .948190E+02 .237719 -.40170E+09 .27352E+01 .10159E+03 .29715E+04 -.5097E+09 .2583E+00 .1079E+02 .3608E+03 .8

-60 ALPH ALPH ALPH ALPH ALPH ALPH ALPH

-60 -60 -60 -60 -60 -60 -60

-60 -60 -60

-.2000001E+01 -.2000002E+01 -.1066668E+02 -.1000002E+02 -.2800006E+02 -.8400020E+02 -.264001E+03 -.858003E+03

-.28600E.04 -.9724E+03

-60 .3000002E+01 .1500001E+01 .1125001E+01 .9375005E+00 .8203129E+00 .7382816E+00

.6464110E+01

.3233255E+02 .1144019E+02

.4849882E+02 .2026507E+03 -.2980907E+03

.6062353E+02 .5066267E+03 .1680097E+04 .3117430E+03

.7072745E+02 .8865967E+03 .5880340E+04 .1592552E+05 .2494919E+04

.7956838E+02 .1329895E+04 .1323076E+05 .7166483E+05 .1633358E+06 .2267488E+05 75252E+02 .182861E+04 .242564E+05 .197078E+06 .898347E+06 .176457E+07 .223971E+06

.237719E+04 .394167E+05 .427003E+06 .291963E+07 .114697E+08 .197788E+08 .234679E+07 29715E+04 .59125E+05 .80063E+06 .72991E+07 .43011E+08 .14834E+09 .22788E+09 .25700E+08 8E+03 .8376E+04 .1361E+06 .1551E+07 .1219E+08 .6305E+08 .1937E+09 .2682E+09 .2913E+08

64001E+03 .676758E+00 .406055E+01 .628418E+00 .439893E+01 .237542 .58914E+00 .47131E+01 .29693E+02 .5564E-01 .5008E+00 .3606E+01 .1

.1000001E+01

.1500002E+01 .1000001E+01

.2250002E+01 .2025003E+01 -.5266673E+01

.2812503E+01 .5062508E+01 .3616079E+01 .2428575E+01

.3281253E+01 .8859389E+01 .1265628E+02 .7382832E+01 .5000009E+01

.3691410E+01 .1328908E+02 .2847662E+02 .3322274E+02 .1630936E+02 .1172730E+02 D6055E+01 .182725E+02 .522071E+02 .913625E+02 .897015E+02 .379506E+02 .302308E+02

.237542E+02 .848366E+02 .197952E+03 .291530E+03 .246679E+03 .916237E+02 .834003E+02 29693E+02 .12725E+03 .37116E+03 .72882E+03 .92505E+03 .68718E+03 .22738E+03 .24182E+03 6E+01 .1803E+02 .6310E+02 .1549E+03 .2621E+03 .2921E+03 .1933E+03 .5764E+02 .7282E+02

^-2.

•(1) • '"1

12,), crq crq

cr r a CD CD C1-' •-• •

•(1)

1-3 gi

•ci) a

n Z't

aa S2O

1•-^ • an C

D f2•^ cD

›•-• • CD 0 ^-t (J) ("6.

CD tri cr ro

(D

a.

rn

C.D

ALPH -45 ...6250004E+00 .3125001E+00 .3125002E+00 .1875002E+00 .4464292E-01 .2321430E+00

ALPH -.8750008E+00 .2734376E+00 .3645836E+00 .3281253E+00 .1562502E4.00 .3038200E-01 .2777780E+00

ALPH -45 -.1312501E+01 .2460938E+00 .4101565E+00 .4921880E+00 .3515630E+00 .1367190E+00 .2237221E-01 .3465912E+00

ALPH -45 -.206250E401 .225586E+00 .451172E+00 .676759E400 .644532E400 .375977E+00 .123047E400 .173528E-01 .451924E+00

-.335157E+01 .209473E+00 .488770E4.00 .879786E+00 .104736E+01 .814617E00 .399903E+00 .112793E+00 .139649E-01 .615105E+00

-.55859E+01 •19638E+00 .52368E+00 .10997E+01 .15710E+01 .15274E+01 .99976E+00 .42297E+00 .10474E+00 .11552E-01 .87132E+00

-.9496E+00 .1855E-01 .5564E-01 .1335E+00 .2226E+00 .2597E+00 .2124E+00 .1198E+00 .4451E-01 .9819E-02 .9762E-03 .1279E+00

ts.., ^--^

N

-30 ALPH

ALPH ALPH ALPH ALPH ALPH ALPH

-30 -30 -30 -30 -30 -30 -30

-30 -30 -30 -30

-.6666667E+00 _.2222223E+0o _.3950618E+00 -.1234568E+00 -.1152264E+00 -.1152264E+00 _.120713E+00 -.130773E+00

-.14530E+00 -.1647E-01

.3333334E+00 .1666667E+00 .1250000E+00 .1041667E+00 .9114586E-01 .8203127E-01 .751953E-01

.698242E-01 .543078E-01 .32 .65460E-01 .58187E-01 .40731 .6182E-02 .6182E-02 .4946E-02

-.3333333E+00

.1851853E-01 -.3703701E-01

.2777779E-01 .2777780E-02 -.2395062E+00

.3472224E-01 .6944450E-02 .5511470E-03 .2292771E-01

.4050928E-01 .1215279E-01 .1929014E-02 .1250287E-03 .3063559E-01

.4557294E-01 .1822918E-01 .4340282E-02 .5626293E-03 .3068888E-04 .3554061E-01 .501302E-01 .250651E-01 .795718E-02 .154723E-02 .168789E-03 .793452E-05 .393585E-01 )1 .325847E-01 .129304E-01 .335233E-02 .548564E-03 .515744E-04 .212847E-05 .428289E-01

.40731E-01 .19396E-01 .62856E-02 .13714E-02 .19340E-03 .15964E-04 .58689E-06 .46337E-01 946E-02 .2748E-02 .1069E-02 .2914E-03 .5480E-04 .6784E-05 .4989E-06 .1653E-07 .5013E-02

$2a 0 (1)

0

-15 ALPH ALPH ALPH ALPH ALPH ALPH

-15 -15 -15 -15 -15 -15

-15 -15 -15

-.5358984E+00 -.1435935E+00 -.2052041E+00 -.5154776E-01 -.3867411E-01 -.3108809E-01

-15 .7179679E-01 .3589839E-01 .2692379E-01 .2243650E-01 .1963193E-01 .1766874E-01

-.4641016E+00 .8591298E-03

.1288695E-02 .1610868E-02 .1879346E-02 .2114265E-02

-.1068360E+00 .2775724E-04

.6939310E-04 .1214379E-03 .1821569E-03

-.1769639E+00 .1186238E-05 .4151835E-05 .9341628E-05

-.2742982E-01

.5796163E-07 -.1703718E-01

.2608273E-06 .3064347E-08 -.1111332E-01

ALPH -15 -.z61801t-01 .161963E-01 .232569E-02 .250466E-03 .171263E-04 .717275E-06 .168539E-07 .170649E-09 -.738973E-02 -.227985E-01 .150395E-01 .251950E-02 .325605E-03 .278303E-04 .155410E-05 .547752E-07 .110922E-08 .985998E-11 -.488453E-02

-.20363E-01 .14099E-01 .26995E-02 .40701E-03 .41745E-04 .29139E-05 .13694E-06 .41596E-08 .73950E-10 .58559E-12 -.31121E-02 -.1855E-02 .1332E-02 .2868E-03 .4942E-04 .5914E-05 .4954E-06 .2910E-07 .1179E-08 .3143E-10 .4978E-12 .3553E-14 -.1808E-03

ALPH ALPH ALPH ALPH ALPH ALPH ALPH -.1

0 0 0 0 0 0 0 -.130920E-01

-.10910E-01 -.9274E-03

ALPH ALPH ALPH ALPH ALPH ALPH ALPH 15

0 0 0 0

-.5000000E+00 .0

-.1250000E+00 .0

-.1666667E+00 .0

-.3906250E-01 .0

-.2734375E-01 .0

-.2050781E-01 .0

-.161133E-01 .0

-01 .0 .0

01 .0 .0

.0 .0

-.227985E-01 -.20363E-01 ",1855E-02

30 30

15 15 15 15 -.5358984E+00 15 -.1435935E+00 15 -.2052041E+00 15 -.5154776E-01 15 -.3867411E-01 15 -.3108809E-01 15 -.261801E-01 , 85E-01 .1503951-C 3E-01 .14099E-01 E-02 .1332E-02

30 30 3o

.0

15 .7179679E-01 .3589839E-0I .2692379E-01 .2243650E-01 .1963193E-01 .1766874E-01

-.5000000E+00

.0 .0

.o .o .0 .0

.0 .0

DIE-01 .161963E-01 .232569 15U395E-01 .251950E-OZ .32 4099E-01 .26995E-02 .40701 32E-02 .2868E-03 .4942E-04

30 3c 3c

.0 .0

.0

-.1250000E+00 .0

.0 .0 .o

.0

.0 .0

.0

-.1666667E+00 .0 .o .o

.0

.0 .0

.0

-.3906250E-01 .0

.o

.0

.0 .0

.0

.0

-.2734375E-01

.0 -.2050781E-01

.0 -.161133E-01

.0 .0 -.130920E-01

.0 .0 -.10910E-01

.0 .0 -.9274E-03

-.4641016E+00

.8591298E-03 -.1068360E+00

•1288695E-02 .2775724E-04 -.1769639E+00

.1610868E...02 .6939310E-04 .1186238E-05 -.2742982E-01

.1879346E-OZ .1214379E-03 .4151835E-05 .5796163E-07 -.1703718E-01

.2114265E-02 .1821569E-03 .9341628E-05 .2608273E-06 .3064347E-08 -.1111332E-01 .232569E-02 .250466E-03 .171263E-04 .717275E-06 .168539E-07 .110649E-09 -.738973E-02 92- .325605E...03 .278303E-04 .155410E-05 .547752E-07 .110922E..08 .985998E-11 ..488453E..02

.40701E-03 .41745E...04 .29139E-05 .13694E-06 .41596E..08 .73950E...10 .58559E-12 -.31121E-02 942E-04 .5914E-05 .4954E-06 .2910E-07. .1179E-08 .3143E-10 .4978E-12 .3553E-14 -.1808E-03

N 2 . a .

.'-t

2 ^CD .-S r)

4-5 n cm O

. A) CD cr cn "

cc^ :2-` ?E:

)-3gi

pa ci) O t1 '—'- ,—t

S=D Z O1

ra Z , • CrC1

CD CD ,-.•

CD ,-t w ö fp.

CD trj CZI"

(D CD ...-, 1-1 SZL ei

ti

ts., b—•

CND

ALPH ALPH ALPH ALPH ALPH ALPH ALPH

30 30 30 30 30 30 30

.-.6666667E+00 -.2222223E+00 -.3950618E+00 -.1234568E+00 -.1152264E+00 -.1152264E+00 -.120713E+00 -.130773E+00 .698242E-01

-.14530E+00 .65460E-01 -.1647E-01 .6182E-02 .6

45 ALPH

ALPH ALPH ALPH ALPH ALPH ALPH

45 45 45 45 -45 45 45

45 45 45

-.1000000E+01 -.5000002E+00 -.1333334E+01 -.6250004E+00 -.8750008E+00 -.1312501E+01 -.206250E+01

.3333334E.00 -.3333333E4,00

.1666667E+00 .1851853E-01 -.3703701E-01

.1250000E+00 .2777779E-01 .2777780E-02 -.2395062E+00

.1041667E+00 .3472224E-01 .6944450E-02 .5511470E-03 .2292771E-01

.9114586E-01 .4050928E-01 .1215279E-01 .1929014E-02 .1250287E-03 .3063559E-01

.8203127E-01 .4557294E-01 .1822916E-01 .4340282E-02 .5626293E-03 .3068888E-04 43554061E-01 .751953E-01 .501302E-01 .250651E-01 .795718E-02 .134723E-02 .168789E-03 .793452E-05 .393585E-01 01 .543078E-01 .325847E-01 .129304E-01 .335233E-02 .548564E-03 .515744E-04 .212847E-05 .428289E-01

.58187E-01 .40731E-01 .19396E-01 .62856E-02 .13714E-02 .19340E-03 .15964E-04 .58689E-06 .46337E-01 .6182E-02 .4946E-02 .2748E-02 .1069E-02 .2914E-03 .5480E-04 .6784E-05 .4989E-06 .1653E-07 .5013E-02

45 .1000000E+01 .5000002E+00 .3750001E+00 .3125001E+00 .2734376E+00 .2460938E+00

LPH 45 -.206250E+01 .225586E+00 .451172E+00 -.335157E+01 .209473E+00 .488770E+00 .879786 -.55859E+01 .19636E+00 .57368E+00 .10997E+01 -.9496E+00 .1855E-01 .5564E-01 .1335E+00 .2

.1729932E-06

.1666668E+00 .1666668E+00

.2500002E+00 .7500008E-01 -.6333336E+00

.3125002E+00 .1875002E+00 .4464292E-01 .2321430E+00

.3645836E+00 .3281253E+00 .1562502E+00 .3038200E-01 .2777780E+00

.4101565E+00 .4921880E+00 .3515630E+00 .1367190E+00 .2237221E-01 .3465912E+00 451172E+00 .676759E+00 .644532E+00 .375977E+00 .123047E+00 .173528E-01 .451924E+00 0 .879786E+00 .104736E+01 .814617E+00 .399903E+00 .112793E+00 .139649E-01 .615105E+00 .10997E...01 .15710E+01 .15274E+01 .99976E+00 .42297E+00 .10474E+00 .11532E-01 .87132E+00 35E+00 .2226E+00 .2597E440 .2124E+00 .1198E+00 .4451E...01 .9819E-02 .9'762E-03 .1279E+00

t^D i--•

.t..

su

0

W rü

«-i-

0

ALPH ALPH ALPH ALPH ALPH ALPH ALPH

60

-.858003E+03 -.28600E+04 -.9724E+03

ALPH ALPH ALPH ALPH ALPH ALPH ALPH

75 75 75 75 75 75 75 75

60 60 60

60 ...2000001E+01 60 •.2000002E+01 60 -.1066668E+02 60 .-.1000002E+02 60 -.2800006E+02 60 -.8400020E+02 60 ...284001E+03

60 .3000002E+01 .1500001E+01 .1125001E+01 .9375005E+00 .8203129E+00 .7382816E+00

.676758E+00 .406055E+01

.628418E+00 .439893E+01 .237542E+02 .58914E+00 .47131E+01 .29693E+02 .1

.5564E01 .5008E+00 .3606E+01 803E

75 75 75

.7464110E+01 ..2785647E+02 -.5544633E+03 ....1939957E+04 -.2027207E+05 -.2269694E+06 ,266220E+07

75 .1392822E+02 .6964110E+01 .5223082E+01 .4352569E+01 .3808497E+01 .3427648E+01

.1000001E+01

.1500002E+01 .1000001E+01

.2250002E+01 .2025003E+01 ...5266673E+01

.2812503E+01 .5062508E+01 .3616079E+01 .2428575E+01

.3281253E+01 .8859389E+01 .1265628E+02 .7382832E+01 .5000009E+01

.3691410E+01 .1328908E+02 .2847662E+02 .3322274E+02 .1630936E+02 .1172730E+02 06055E+01 .182725E+02 .522071E+02 •9136251+02 .897015E+02 .379506E+02 .302308E+02

542E+02 .848366E+02 .197952E+03 .291530E+03 .246679E+03 .916237E+02 .834003E+02 +02 .12725E+03 .37116E+03 .72882E+03 .92505E+03 .68718E+03 .22738E+03 .24182E+03

.1803E+02 .6310E+02 .1549E+03 .2621E+03 .2921E+03 .1933E+03 .5764E+02 .7282E+02

LPH 75 ii266220E+07 .314201E+01 .875252E+02 - .322903E+08 .291758E+01 .948190E+02 .237719 - .40170E+09 .27352E+01 .10159E+03 .29715E+04 - .5097E+09 .2583E+00 .1079E+02 .3608E+03 .8

.6464110E+01

.3233255E+02 .1144019E+02

.4849882E+02 .2026507E+03 -.2980907E+03

.6062353E+02 .5066267E+03 .1680097E+04 .3117430E+03

.7072745E+02 .8865967E+03 .5880340E+04 .1592552E+05 .2494919E+04

.7956838E+02 .1329895E+04 .1323076E+05 .7166483E+05 .1633358E+06 .2267488E+05 75252E+02 .182861E+04 .242564E+05 .197078E+06 .898347E+06 .176457E+07 .223971E+06

.237719E+04 .394167E+05 .427003E+06 .291963E+07 .114697E+08 .197788E+08 .234679E+07 29715E+04 .59125E+05 .80063E+06 .72991E+07 .43011E+08 .14834E+09 .22788E+09 .25700E+08 3E+03 .8376E+04 .1361E+06 .1551E+07 .1219E+08 .6305E+08 .1937E+09 .2682E+09 .2913E+08

N 2 ' . 5: CD A) 2

CD 0 n

crq•crq O. . --: •

cr Al co

a^ el . m

1-3 gi

z3cn crq ict

s2u

Cri:1 m

,--, ... .--t pp a (t

, CD 1-t (i) 0 FD.

^-,-.

CD tll CD CD ,i pi.

rD.

t^., .-+

01