5. 4 0 . 1149 6 8 2 E ー 0 ろこ 0 第当 1 4 1 0 4 5 E ー 0 0 を 160-6 6 ろ E ニ 0 一 . 5 日 9955nE ー℃当 分割 C による計算結果 出力例 5.3 ( その 3 ) いろいろな分割で計算してみる 167 変位 1 2 ろ 4 5 6 7 8 9 11 1 2 1 ろ 14 15 16 17 18 19 2 0 2 2 2 ろ 2 4 2 5 2 6 2 7 28 2 9 ろ 0 ろ 1 ろ 2 ろ 4 ろ 5 ろ 6 5 7 点番号 工成分 D X 0 . 2 42 0 7 5 5 E ー 0 ろ 0 . 2 6 2 7 9 2 2 E ー 0 ろ 0 . 2 7 810 5 7 E - 0 ろ 0 . 1 ろ 7 2 5 0 4 E ー 0 ろ 0 . 15 6 019 ろ E- 0 . 1819 ろろ 7 E ー 0 ろ 0 . 19 5 5 6 8 4 E ー 0 ろ 0 . 7 7 0 7 1 1 5 E ー 0 4 0 . 8 7 6 4 7 4 5 E ー 0 4 0 . 9 8 ろ 6 418 E ー 0 4 0 . 1 ろ 5 9 0 2 6 E ー 0 ろ 0 . 17 2 5 46 2 E ー 0 ろ 0 . 215 9 2 ろ 6 E ー 0 4 0 . 2 5 4 2 8 5 4 E ー 0 4 0 . ろ 0 6 2 6 0 2 E ー 0 4 成分 し Y - 籖当 -0 4 0 弖生 9- E ニ 0 互 ー 0 . 8 9 4 0 2 2 7 E - 0 5 0- 亠 00 0 0 0 0 0 F ー 0 . 4 7 5 7 4 8 5 E ー 0 4 0 ーー℃ . 2 R ろ 9 2 5 a F ー 4 ー 0 . 1 0 5 9 98 7 E ー 0 4 0 上立 000 ー 00.0. E 一 00 - 0 . 6 6 5 0 8 2 5 F ー 0 4 ー 0 亠 56 2 6 216 E -0 当 ー 0 . 4 4 6 2 8 0 5 E ー 0 4 ー 0 . 18 716 2 4 E ー 0 4 -0 . ろ -8 0-6-606E ュ 0 生 0 ・ 1 6 2 4 0 8 1 E ー 0 ろ 0 . 1 5 7 0 0 7 ろ E ー 0 ろ 14 0 2 4 8 4 E ー 0 ろ 0 を 1 17 5 4 ろ 7 E ー 0 ろ 0 . 8 3 984 7 1 E ー 0 4 0 . 4 5 0 61 ろ 9 E ー 0 4 0 を 14 5 42 ろ 5 E ー 0 4 0 . 16 ろ 8 5 8 5 E ー 0 ろ 0 . 14 7 8 0 6 5 E ー 0 ろ 0 . 12 4 5 912 E ー 0 ろ 0 . 9 ろ 0 0 17 6 E ー 0 4 0 . 6 4 416 8 8 E ー 0 4 0 . 8 ろ 0 9 2 2 8 E ー 0 4 0 0 0 0 0 0 0 0 0 0 ー 0 亠 242 150 丘こ℃ 5 ー 0 . 9 ろ 6 8 5 0 ろ E ー 0 ー 0 亠 3-964028E ー 0 4 ー 0 . 6 0 7 0 ろ 48E ー 0 4 ー 0 . 7 ] ろ 3 5 弖 2 E - 0 ー 0 . 4 6 6 ろ 617 E ー 0 4 ー 0 - -57 7 9 0 92E ニ℃ 4 ー 0 亠 258 0 ろ 66- E ニ 04 ー 0 . 4 1 9 2 5 9 9 E ー 0 4 ー 0 . 4 0 1 1 ろ 7 8 E ー 0 4 ー 0 . 今 9 7 0 17 2 E ー C 4 ー 0 . 5 8 0 8 ろ 6 2 E ー 0 4 ー 0 第ち 65 2 9 E ー 0 4 ー 0 . 7 ろ 5 15 4 2 E ー 0 / を ー 0 亠 717 豆 00 ら E ー 0.4 ー 0 . 7 4 5 ろ 0 8 9 E ー 0 4 ー 0 8 2 8 ろ 9 5 ろ E - 0 ー 0 . 9 0 6 7 9 9 7 E - 0 4 ー 0 ュ 987 今 595E ー 0 4 ー 0 . 1 18 5 5 4 ろ E ー 0 4 0 0 0 0 - 0 . 12 8 4 9 ろろ E ー 0 ー 0 ュ ~ 日生旦 0 弖 E ー 0 4

140 4. 実習用プログラム SUBROUTINE JNEXPR く VA しリ計 CO M M 0 N / I B U S / K , L ( 81 ) , I D E N T ( 2 5 6 ) CAL L SK I PB ( 0 ) A R E G : 0 . 0 S G N = 1 , 0 I F ( ( L ( K ) ・ E Q ・ 0 ) . 0 R . ( し ( K ) . EQ . - 6 ) . OR . ( K . GE . 80 ) ) TO 2 TO 1 1 I F ( K . GE . 8 0 ) G 0 T 0 C I F B L A N K I F ( L ( K ) . E Q . - 6 ) GO C T H E F I R S T T E R M I F ( L ( K ) を NE . - 2 ) GO S 6 N : - 1 . 0 K = K + 1 2 CONTINUE K : K + 1 1 0 P : ー L ( K ) C BINARY OPERATOR 1 G 0 T 0 1 ろ CONTINUE C TEST FOR RESTRICTOR A R E G : S G N V A L C A L L 1 N T E R M ( V A L ) 1 F ( 1 0 P . NE . ろ ) GO TO 5 I F ( L ( K ) . NE . - ろ ) GO K : K + 1 1 0 P : 5 5 C 0 N T I N IJ E C T H E N E X T T E R M CALL INTERM(DREG) TO 5 C A R E 6 0 P D R E G I F ( 1 0 P . E Q . 1 ) I F ( I OP . E Q . 2 ) I F ( I OP . EQ . ろ ) I F ( I OP . EQ . 4 ) 1 F ( I OP . E Q . 5 ) G 0 T 0 ろ 1 C 0 N T I N IJ E E N D R E T IJ R N V A L し E : A R E G A R E G = A R E G + D R E G A R E G : A R E 6 - D R E G A R E G : A R E G 朞 D R E G A R E G : A R E G / D R E G A R E G : A R E G 畏 D R E 6

4. 6 プログラムのリスト I F ( L ( K ) . NE . - 6 ) CALL K = K + 1 CALL INEXPR(FINALV) D X N 0 W = 1 . 0 CAL L S K I PB ( 0 ) I F ( K . GE . 8 0 ) GO TO 9 I F ( し ( K ) . NE . - 6 ) CALL K : K + 1 CALL INEXPR(DXNOW) 137 D U M P ( 2 H D 0 , K ) DUMP(2HDO'K) I F ( DX NOW . EQ . 0 . 0 ) CALL DUMP ( 4HDX : 0 , 0 ) 9 ( 0 N T I N U E D X ( L E V E L ) : D X N 0 ー = ( ( F I N A LV- R ( L OC ) ) / D X N 0 ) 益 1 . 0 0 0 0 1 ー F ( ( I, GT . 10 0 0 0 ) , 0 R . ( I. L T . 0 ) ) CALL DUMP ( 2HDO , い M A X K N T ( L E V E L ) : I K 0 U N T ( L E V E L ) : 0 GO TO 1 DATA CARD 5 C 0 N T 1 N U E L T Y P E : 0 10 CONTINUE DO 11 J = 1 , N 0 L D V A L = D A T A ( J ) ( A L L I N E X P R ( V A L IJ E ) D A T A ( J ) : V A L し E I F ( K . GE . 80 ) GO TO 12 11 C 0 N T I N U E 12 C 0 N T I N U E L 00 P C 0 N T R 0 L I F ( L E V E L . L E . 0 ) R E T U R N 1 ろ C 0 N T I N IJ E K L : K 0 U N T ( L E V E L ) MX=MAXKNT(LEVEL) I F ( KL . GE . MX ) GO TO 14 R ( L 0 C ) : S T A R T V ( L E V E L ) + F L K D X ( L E V E L ) FLK=KOUNT(LEVEL) L 0 C = I X ( L E V E L ) KOUNT(LEVEL)=KL + 1 ( NOT YET LEVEL=MAXLEV R E T U R N C E N D 0 F A L 0 0 P 14 ( 0 N T I N U E L 0 C : I X ( L E V E L ) R ( L 0 C ) = S T A R T V ( L E V E L ) K 0 U N T ( L E V E L ) : 0 L E V E L : L E V E L ー 1 ー F ( LEVEL . GE . 1 ) GO TO L E V E L = 0 M A X L E V = 0 R E T U R N E N D 1 ろ

136 4. 実習用プログラム SUBROUTINE INCAÄb(DATAtN'LTYPE)l, D I ME NS I ON DATA ( N ) COMMON /FBUS/ R(256),OLDVAL ( 0 M M 0 N / ー 8 US / K , L ( 81 ) , I DE N T ( 2 5 6 ) C S P E ( I F I C A T 1 0 N S 1 X ( 7 ) , 0 X ( 7 ) , S TA RTV ( 7 ) , K 0 IJ N T ( 7 ) , M A XK NT ( 7 ) WR I T E ( L 0 IJT , 6 01 ) 7 3 C 0 N T I N U E 6 0 T 0 5 5 0 1 F 0 R M A T ( 8 0 A 1 ) READ ( L I N , 5 0 1 , E ND : 7 ろ ) ( L ( い , I : 1 , 8 0 ) 1 C 0 N T I N U E ( R E A D I N I F ( LE VEL . GT . 0 ) GO TO 2 ( C H E C K I F リ N D 0 ー M 0 D E を D A T A K A R D N 0 / 1 / D A T A L E V E L / 0 / , L I N / 5 / , L 0 IJ T / 6 / C D 0 5 T A R T V ( L E V E L ) : R ( L 0 ( ) I X ( L E V E L ) : L 0 C ( A L L L E T ( L 0 C ) IF(LEVEL .6T物7) CALL DUMP(4HLEVL'LEVEL) M A X L E V : L E V E L 8 L E V E L : L E V E L + 1 6 0 T 0 1 ( A L L L E T ( L 0 C ) 7 ( 0 N T I N U E C L E T 6 0 TO 1 0 L T Y P E : N A M E 4 1 5 ) 60 T 0 8 I F ( NAM E . E Q . I F ( NA M E . E Q . 12 0 5 2 0 ) 60 TO 7 ( A L L I N N A M E ( N A M E ) I F ( L ( K ) . LE . 0 ) 60 TO 1 I F ( L ( K ) . 6 E . ろ 0 ) CA L L DIJM P ( 4H CAR D , K ) ー F ( K . GE . 8 0 ) GO TO 1 CALL SK I PB ( 0 ) C ( 0 N T R 0 L C A R D I F ( し ( 1 ) . NE . 5 ) GO TO 5 2 ( 0 N T I N U E C T H E F 1 R S T ( O L U M N CALL TCODE C ( 0 D E ( 0 N V E R S 1 0 N K A R D N 0 : K A R D N 0 + 1 6 0 2 F 0 R M A T ( 1 H , 8 0 A 1 , 10 X 戸 CARD NO , 1 4 ) ろ WR I TE ( LOIJT , 6 02 ) ( L ( い , I : 1 , 8 0 ) , KARDNO ( E ( H 0 P R I N T I N G R E T U R N LTYPE=-I 601 FORMAT い畏益朞 EOF つ

130 8 7 610 1 0 0 16 18 4. 実習用プログラム C 0 N T I N U E D ( 1 , 1 ) : D 1 D ( 1 , 2 ) : D 2 D ( 2 , 1 ) : D 2 D ( 2 , 2 ) : D 1 D ( ろ , ろ ) : D 3 CAL L MULT ( D , B , DB , ろ , 3 , 6 ) ( A L L WUL T ( B , DB , E K , 6 , ろ , 6 ) C A L L S IJ L T ( E S , E K , 6 ) ME : 2ÅMAX 0 ( I ABS ( I -J ) , I ABS ( J-K ) , 1 A8S ( K - い ) + 2 I F ( ME , GT . MM ) GO TO 9 0 0 M : MAX 0 ( M , M E ) LL ( 1 ) : 2 益 I - 1 L L ( 2 ) : 2 I L L ( 3 ) : 2 益 J -1 L L ( 4 ) : 2 朞 J L L ( 5 ) : 2 K -1 L L ( 6 ) : 2 K I F ( LST R . EQ . 1 ) DO 7 1 E : 1 , 6 D 0 7 J E : 1 , 6 I T : L L ( I E ) J T : L L ( J E ) G 0 T 0 1 0 0 I F ( I T . LT . JT ) GO TO 7 1 8 M : ー T - J T + 1 T K ( I B M , JT ) : TK ( I B M , JT ) + E K ( I E , J E ) ( 0 N T I N U E ー F ( DUMP を GE . 1 . 0 ) WR I TE ( 6 , 610 ) L , E K F 0 R MA T ( 1 H 0 , , E L E ME NT ・ , 1 5 / 6 ( 1 H 6 F 15 . 5 / ) ) R E T IJ R N ( 0 N T I N U E DO 16 1 E : 1 , 6 1 T= LL (I E ) EU ( I E ) : TU ( I T ) CAL L VU L T ( E K , E IJ , E F , 6 , 6 ) CAL L V UL T ( D B , E U , D B U , ろ , 6 ) DO 18 1 E = 1 , 6 I T : L L ( I E ) BC い T ) : BC ( I T ) + EF ( I E ) S X = D B U ( 1 ) SY=DBU(2) T X Y : D B U ( ろ )

4.5 入力データの作成法と結果の読みかた プログラム 4.1 D D L の使用例ーー図 5.4 の ( a ) ~ ( e ) を計算するための 117 益 BUNKATU SUU NO SITEI L E T L Y : 2 0 0 L E T L X : 6 0 0 CANTILEVER BEAM データを記述したもの Y ロ IJ S 0 N 0 T E I G I I N P 1 + J + 1 ー畏 D X J 畏 D Y DO J : 0 , N kD 0 1 : 0 , M 朞 S E T T E N Z A H Y 0 U N P N E N 日 N F 朞 L E T N F : N P 1 L E T N B : N P 1 L E T N E : M N L E T N P : M P 1 朞 N P 1 畏 MA ー S UU N ロ S I TE I L E T D Y : L Y / N ÅLET DX=LX/M ÅLET NPI=N + 1 kL E T M P 1 : M + 1 *LET N=I L E T M : ろ ろ 1 L E T N U : 0 . ろ L E T E : 2 1 0 0 0 畏 L ET H : 2 . 0 朞 L E T P : N P 1 QUAD DO I : 1 , M 朞 DO J : 1 , N I 畏 N - N + J kkKOTEI DO I = 0 , N ZYOUKEN I P ー P + J L E T N F = N P 1 畏 L E T N B : N P 1 kLET NE=M%N *LET NP=MPI%NPI L E T D Y : L Y / N ÅLET DX=LX/M kLET NPI=N + 1 L E T M P 1 = M + 1 し E T N : 2 L E T M : 6 6 2 N P 1 0 D F / 2 1 0 D F / 2 L E T D F : 1 / N k%KAZYlJU I 畏 P + J I 畏 P + J + 1 I 畏 P ー P + J + 1 E N U H

4. 6 プログラムのリスト 141 1 2 ろ 1 2 2 1 SUBROUTINE INTERM(V) ( 0 M M 0 N / I B U S / K , L ( 81 ) , I D E N T ( 2 5 6 ) COMMON /FBUS/ R(256)'OLDVAL 1 F ( ( L ( K ) . GE . 1 ) . AND . ( L ( K ) . L E . 2 6 ) ) I F ( L ( K ) . E Q . -10 ) GO TO 2 CALL INREAL(VAL) V = V A L R E T IJ R N C 0 N T I N U E ( A L L I N N A M E ( N A M E ) C A L L L 0 C A T N ( N A M E , L ロ ( ) V : R ( L 0 C ) R E T U R N ( 0 N T I N U E K : K 申 1 V=OLDVAL R E T U R N E N D SUBROUTINE INNAME(NAM) C 0 M M 0 N / ー B U S / K , L ( 81 ) , I D E N T ( 2 5 6 ) 1 F ( ( L ( K ) . L E . 0 ) . OR . ( L ( K ) . G E . 2 7 ) ) 1 C A L L D U M P ( 4 H N A M E , K ) N A M : N A M 1 0 0 + L ( K ) DO 1 1 = 1 , ろ N A M : 0 G 0 T 0 1 K : K + 1 I F ( L ( K ) . L E . 0 ) C 0 N T I N U E D 0 2 ー : K , 8 0 I F ( L ( K ) . L E . 0 ) C 0 N T I N U E K : 81 R E T U R N RETURN END R E T IJ R N 6 0 T 0 ろ 5U8R0 リ〒 'tNE ー NREAL( VALVE メ ( 0 M M 0 N / I B IJ S / K , L ( 81 ) , 1 D E N T ( 2 5 6 ) K I N T : K S I G N : 1 . 0 E X P S 6 N : 1 . 0 V= 0 . 0 E X P : 0 . 0 I F ( L ( K ) . EQ . -1 ) I F ( L ( K ) . E Q . ー 2 ) 6 0 T 0 ろ S 1 6 N = ー 1 . 0 K : K + 1 ( 0 N T I N U E I F ( L ( K ) . LT . ろ 0 ) 6 0 T 0 G 0 T 0 G 0 T 0 1 2 5

4.6 プログラムのリスト 88 CONTIN し E CALL S MTRM ( L , J , K , Q , EE ( L ) , pp ( L ) , HH ( L ) , RHO ( L ) ) CAL L SMTRM ( L , Q , J , E E ( L ) , pp ( L ) , HH ( L ) , RHO ( L ) ) 89 CONTINUE R E T IJ R N END S リ 8R0 リロ NE ' SMTRM(L'I'J'K,E'P'H'R) 、 CO M M 0 N X X ( 2 5 0 ) , YY ( 2 5 0 ) 129 1 , B ( ろ , 6 ) , D ( ろ , ろ ) , DB ( ろ , 6 ) , L L ( 6 ) , E K ( 6 , 6 ) , EU ( 6 ) , EF ( 6 ) , DB IJ ( ろ ) 2 , TK ( 7 0 , 5 0 0 ) , TU ( 5 0 0 ) , TF ( 5 0 0 ) , BC ( 5 0 0 ) ろ , ST ( 2 , 2 ) , T ( 2 , 2 ) , RT ( 2 , 2 ) , STT ( 2 , 2 ) 4 , I I ( 2 5 0 ) , JJ ( 2 5 0 ) , KK ( 2 5 0 ) , QQ ( 2 5 0 ) , E E ( 2 5 0 ) , p p ( 2 5 0 ) , L Q ( 2 5 0 ) 5 B ( 1 0 0 ) 冂 X ( 1 0 0 ) , I Y ( 1 0 0 ) 6 , M M , N N , M , N , S 7 , H H ( 2 5 0 ) , RHO ( 2 5 0 ) , E M ( 2 5 0 ) 8 , L S T R , D IJ M P D A T A P I , R T 0 D / 3 . 1415 9 , 5 7 . 2 9 5 8 / X ー =XX ( I ) X J : X X ( J ) Y J : Y Y ( J ) X K : X X ( K ) Y K : Y Y ( K ) S : ( X J - X い ( Y K - Y い - ( YJ - Y い益 ( XK - X い E S : A B S ( S ) / 2 , 0 1 F ( 5 . EQ . 0 . 0 ) GO TO 9 0 2 ー F ( H . NE . 0 . 0 ) ES= ESkH E M I : E S 益 R E M ( い : E M ( い + E M I E M ( J ) : E M ( J ) + E M I E M ( K ) : E M ( K ) + E M I B ( 1 , 1 ) : ( YJ-YK ) / 5 B ( 1 , ろ ) : ( YK-Y I ) / 5 B ( 1 , 5 ) : ( Y ー -YJ ) /S B ( 2 , 2 ) : ( XK-XJ ) / 5 B ( 2 , 6 ) : ( XJ-X I ) / 5 DO 6 J B : 1 , 5 , 2 B ( 3 , JB ) : B ( 2 , JB + 1 ) 6 B ( ろ , JB + 1 ) =B ( 1 , JB ) ー F ( H . EQ . 0 . 0 ) GO TO 5 D 1 : E / ( 1 . 0 - P kP ) D 2 : D 1 益 P D 3 : 0 . 5 E / ( 1 . 0 + P ) G 0 T 0 8 5 C 0 N T I N IJ E D 1 : E% ( 1 . 0 -P ) / ( ( 1 . 0 + P ) 朞 ( 1 . 0 - 2 . 0 %P ) ) D 2 = E 益 P / ( ( 1 . 0 + P ) 畏 ( 1 . 0 - 2 . 0 ) ) D 3 : 0 . 5 畏 E / ( 1 . 0 + p )

124 1 F ( LTYPE . EQ . 0 ) 60 TO 62 C C CALL INCARD(CARD'12'LTYPE) 5 9 ( ロ N T ー N IJ E I F ( LTYP E . E Q . 17 2 101 ) L QUA D = 1 1 F ( L T Y P E . E Q . 2 0 18 0 9 ) L Q U A D : 0 I F ( L T Y P E . E Q . 12 0 914 ) L Q U A D : - 1 4. 実習用プログラム G 0 T 0 5 9 6 2 ( 0 N T I N U E N 0 : ( A R D ( 1 ) 1 F ( L . NE . NO ) L Q ( L ) = L Q U A D 6 0 T 0 910 I I ( L ) = CARD ( 2 ) JJ ( L ) : CARD ( ろ ) KK ( L ) : C AR D ( 4 ) 6 0 ろ F 0 R MAT ( / / / 1 H 1 , 4 X , 1 H L , 4X WR I TE ( 6 , 6 0 ろ ) 111 ( 0 N T I N U E 61 C 0 N T I N U E G 0 TO 6 1 HH ( L ) : 0 . 0 Q Q ( L ) : 0 R HO ( L ) = C A R D ( 6 ) P P ( L ) : CAR D ( 5 ) E E ( L ) =CARD ( 4 ) 6 6 ( 0 N T I N U E L I N E A R E L E M E N T G 0 T 0 61 RHO ( L ) =CARD ( 9 ) HH ( L ) : CA RD ( 8 ) P P ( L ) : C A RD ( 7 ) EE ( L ) = CARD ( 6 ) QQ ( L ) : CARD ( 5 ) 6 0 C 0 N T I N U E Q IJ A D 1 E L E M E N T GO T 0 61 QQ ( L ) : 0 R H 0 ( L ) = ( A RD ( 8 ) HH ( L ) : CAR D ( 7 ) PP ( L ) : CAR D ( 6 ) E E ( L ) = CAR D ( 5 ) T R I A N G IJ L A R E L E M E N T I F ( LQUAD . EQ . 1 ) GO TO 6 0 I F ( L Q UAD . E Q . - 1 ) G 0 T 0 6 6 2 , 1 1 X , 1 HH , 1 0 X , 3H R HO ) 1 , 1 H I , 4 X , 1 H J , 4X , 1 H K , 4X , 1 HL , 11 X , 1 HE , 18X , 1 H P

3. 構 6 0 9 F 0 R M AT ( 5 X , 1 H L , 7 X , 2 H 5 x , 1 ろ X , 2 HS Y , 1 2 X , ろ H 1 x y 1 , 1 ろ X , 2 H S 1 , 1 ろ X , 2 H S 2 , 12X , ろ HT H D , 1 ろ X , 2 H S E ) 造 解 析 16 1 7 1 9 2 0 21 2 ろ 1 0 1 D 0 D 0 D 0 D 0 圧 ( I E . LE . ろ ) D8 ( I E , JF ) *C8 い E , JF , L ) EK ( I E , JE ) =CK い E , JE , L ) 1 7 JE : 1 , 6 LL ( I E ) : LC ( I E , L ) 1 7 に : 1 , 6 EU ( I E ) : TU ( I T ) I T : LL ( I E ) 1 6 1 E : 1 , 6 1 0 1 L : 1 , N E CALL VUL T ( D B , E U , D 8 U , ろ , 6 ) CAL L VUL T ( EK , EU , E F , 6 , 5 ) C 0 N T I N U E S X : 0 B U ( 1 ) S Y : D B リ ( 2 ) T X Y : D B IJ ( ろ ) T H : P I / 4 . 0 I F ( ABS ( S X-SY ) . LT . 1 . 0 E ーろ 0 ) R : 2 . 0 T X Y / ( S X ー S Y ) T H : 0 . 5 A T A N ( R ) C 0 N T I N IJ E 5 =S ー N ( TH ) C : C 0 S ( T H ) G 0 TO 21 S 1 ニ C 畏 C 畏 S X + 2 . 0 畏 C 畏 S 畏 T Y + S 畏 S 畏 S Y 5 2 : S S 畏 S X ー 2 . 0 C 畏 S 畏 T Y Y + 0 に 5 Y 5 E : S Q R T ( 5 1 2 + S 2 Ⅲ ( 2 ー S 1 畏 S 2 ) T H D : T H R T 0 D WR I TE ( 6 , 2 ろ ) L , SX , SY , T XY , S 1 , S 2 , THD , SE F OR M A T ( 1 H , 1 5 , 7 G 15 . 7 ) 5 T 0 P C 0 N T I N U E C E R R 0 R M E S 5 A 6 E 910 WR I T E ( 6 , 911 ) 9 1 1 F 0 R M A T ( 1 5 H 5 E Q U E N C E 1 0 ろ S T 0 P E N D E P R ロ R