Created
May 29, 2017 06:44
-
-
Save relyks/4e63e4fbe73e2a0ad9d681624eb96bc7 to your computer and use it in GitHub Desktop.
Formac Example Code
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Partial program for calculating f and g. | |
$IBFMC TS02 | |
SYMARG | |
DIMENSION F(100), G(100), OUTPUT(22) | |
C FORMAC DATA DECLARATIONS. | |
ATOMIC MU, SIGMA, EPSI, T | |
DEPEND (MU, SIGMA, EPSI/T) | |
LET F(1) = 1. | |
LET G(1)=0. | |
DO 13 I=l, 100 | |
C CODE TO GENERATE F(I+l) AND G(I+l) FROM F(I) AND G(I). | |
5 LET F(I+l) = SUBST FMCDIF(F(I),T,l)-MU*G(I), LIST | |
LET F(I+1) = EXPAND F(I+1) 10 LET G(I+l) = SUBST(FMCDIF(G(I),T,l)+F(I)), LIST | |
10 LET G(I+l) = EXPAND G(I+l) | |
LIST PARAM (FMCDIF(MU,T,S), -3. *MU*SIGMA), (FMCDIF(SIGMA,T,l), EPSI-2. *SIGMA**2), (FMCDIF(EPSI,T,l), -SIGhlA*MU-SIGMA*2. EPSI) | |
13 CONTINUE | |
Output of Newcomb operators programs. | |
BCDPRINT input | |
*XS21, LS48* N, K*XS21, LS10, DL15, 27, DL3, LS48* 6, 1*LS10, XS3* = *XS2*0.424 91319K+0.13270400K*LS36*2*LS10*-0.58263887K*LS36*3*LS10*-0.52940538K*LS3 6*4*LSlO*-O.l7751735K*LS36*5*LSlO**VP2, XS22*-0.26041666X10*LS36*--1*LS10* K*LS36*6*LS10*-0.13888888X10*LS36*-2*LS1O*K*LS36*7*LS1O**VP2, XS22*f(-0.2 | |
Phototypeset output Zs,lnsk=0.42491319k+0.1327040Ok2-O.58263887k~-O.52940538k4-O.17751735k6-O.26O41666X1Oo1k6-O.13888888X10-2k7 | |
+ (-0.22146267-0.15382270X10k-0.13371744X10k2-0.21294488k3+0.11219618k4 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment