//JOB(ERCC07,STEPHENS,T=200,L=10000) //OPTIONS(QUOTES,OPT) //ALGOL 'COMMENT' *************** * * * ALGG042 * * * *************** ; 'BEGIN' 'COMMENT' !***FY31****BENDDME************************ ; 'BEGIN' 'COMMENT' A6,A30; 'REAL' E1,E2,E3,Q1,Q2,Q3,T1,T2,T3,QQ,TT,D,ANS,R,AT,BT,CT; 'INTEGER' X,U,V,P,S,I,J; 'ARRAY' T[1:10],Q[1:18],E[1:10,1:18],A[1:3],B[1:3],C[1:3],ET [1:3,1:3]; 'FOR' I:=1 'STEP' 1 'UNTIL' 10 'DO' T[I]:=READ ; 'FOR' J:=1 'STEP' 1 'UNTIL' 18 'DO' Q[J]:=LN(READ); 'FOR' I:=1 'STEP' 1 'UNTIL' 10 'DO' 'FOR' J:=1 'STEP' 1 'UNTIL' 18 'DO' E[I,J]:=LN(READ); QQ:=LN(READ); TT:=READ; 'IF' TT>T[9]-(T[9]-T[8])/2 'THEN' 'BEGIN' U:=9; 'GOTO' SS1; 'END'; 'IF' TTT[P] 'AND' TT-T[P]<(T[P+1]-T[P])/2 'THEN' U:=P 'ELSE' 'G OTO' RPT 1; SS1: 'IF' QQ>Q[17]-(Q[17]-Q[16])/2 'THEN' 'BEGIN' V:=17; 'GOTO' SS2; 'END'; 'IF' QQQ[S] 'AND' QQ-Q[S]<(Q[S+1]-Q[S])/2 'THEN' V:=S 'ELSE' 'GOTO' RPT 2; SS2: Q1:=Q[V-1]; Q2:=Q[V]; Q3:=Q[V+1]; T1:=T[U-1]; T2:=T[U]; T3:=T[U+1]; 'FOR' I:=1 'STEP' 1 'UNTIL' 3 'DO' 'FOR' J:=1'STEP' 1 'UNTIL' 3 'DO' ET[I,J]:=E[U-2+I,V-2+J]; D:=Q1**2*(Q2-Q3)-Q2**2*(Q1-Q3)+Q3**2*(Q1-Q2); 'FOR' X:=1 'STEP' 1 'UNTIL' 3 'DO' 'BEGIN' A[X]:=(ET[X,1]*(Q2-Q3)-ET[X,2]*(Q1-Q3)+ET[X,3]*(Q1-Q2)) /D; B[X]:=(-ET[X,1]*(Q2**2-Q3**2)+ET[X,2]*(Q1**2-Q3**2) -ET[X,3]*(Q1**2 - Q2**2))/D; C[X]:=(ET[X,1]*Q2*Q3*(Q2-Q3)-ET[X,2]*Q1*Q3*(Q1-Q3)+ET[X,3]*Q1*Q2 * (Q1-Q2))/D; 'END'; E1:=A[1]*QQ**2+B[1]*QQ+C[1]; E2:=A[2]*QQ**2+B[2]*QQ+C[2]; E3:=A[3]*QQ**2+B[3]*QQ+C[3]; R:=T1**2*(T2-T3)-T2**2*(T1-T3)+T3**2*(T1-T2); AT:=(E1*(T2-T3)-E2*(T1-T3)+E3*(T1-T2))/R; BT:=(-E1*(T2**2-T3**2)+E2*(T1**2-T3**2)-E3*(T1**2-T2**2))/R; CT:=(E1*T2*T3*(T2-T3)-E2*T1*T3*(T1-T3)+E3*T1*T2*(T1-T2))/R; ANS:=EXP(AT*TT**2+BT*TT+CT); WRITE TEXT('(' ANSWER = ')'); PRINT(ANS,0,4); NEWLINES(1); SS3: 'END'; 'END' //RUN 500;1000;1500;2000;2500;3000;3500;4000;4500;5000; 0.1;0.12;0.15;0.2;0.25;0.3;0.4;0.5;0.6;0.8;1.0;1.2;1.5;2;3;5;10;20; 0.132;0.093;0.073;0.054;0.039;0.029;0.021;0.015;0.011;0.008; 0.145;0.104;0.082;0.061;0.045;0.034;0.025;0.018;0.014;0.010; 0.170;0.117;0.094;0.072;0.054;0.041;0.030;0.023;0.017;0.013; 0.190;0.138;0.112;0.086;0.066;0.051;0.039;0.030;0.023;0.018; 0.210;0.155;0.126;0.100;0.077;0.060;0.047;0.036;0.028;0.022; 0.225;0.166;0.141;0.112;0.088;0.068;0.054;0.042;0.033;0.026; 0.255;0.192;0.163;0.130;0.105;0.083;0.066;0.053;0.042;0.034; 0.275;0.212;0.184;0.150;0.120;0.097;0.078;0.063;0.050;0.040; 0.300;0.230;0.200;0.163;0.133;0.107;0.087;0.071;0.058;0.047; 0.325;0.260;0.230;0.189;0.155;0.127;0.105;0.085;0.070;0.058; 0.350;0.280;0.250;0.210;0.174;0.143;0.118;0.096;0.080;0.066; 0.370;0.302;0.270;0.228;0.190;0.157;0.130;0.108;0.090;0.075; 0.400;0.328;0.295;0.250;0.210;0.170;0.145;0.123;0.103;0.087; 0.430;0.360;0.330;0.280;0.238;0.200;0.170;0.144;0.120;0.104; 0.470;0.405;0.375;0.326;0.280;0.236;0.200;0.170;0.144;0.123; 0.540;0.460;0.435;0.380;0.330;0.280;0.240;0.205;0.175;0.150; 0.610;0.540;0.515;0.460;0.395;0.340;0.295;0.255;0.220;0.190; 0.690;0.600;0.590;0.530;0.460;0.400;0.350;0.305;0.268;0.232; 20; 3000;