Skip to content

Instantly share code, notes, and snippets.

@relyks
Created May 29, 2017 06:44
Show Gist options
  • Save relyks/4e63e4fbe73e2a0ad9d681624eb96bc7 to your computer and use it in GitHub Desktop.
Save relyks/4e63e4fbe73e2a0ad9d681624eb96bc7 to your computer and use it in GitHub Desktop.
Formac Example Code
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