NATL INST, OF STAND & TECH

A111D7 E373A1

Hi

i

NtST

PUBLICATIONS

WBPl

^^B

Bi

Bl

BJ Bfll ^h

jam

BBBBBBBBBSBI Hi

HB



National Bureau of St»

UBrtty. Ml Mm6^ <- NBS TECHNICAL NOTE 550

A Systems

Programmer's Guide for

Implementing OMNITAB II

National Bureau of Standards

Library, E-Ql Admin. Bide. UNITED STATES DEPARTMENT OF COMMERCE

Maurice H. Starts, Secretary

NATIONAL BUREAU OF STANDARDS • Lewis M. Branscomb, Director

NBS TECHNICAL NOTE 550

ISSUED NOVEMBER 1970

Nat. Bur. Stand. (U.S.), Tech. Note 550, 43 pages (Nov. 1970) CODEN: NBTNA

A Systems Programmer's Guide for

Implementing OMNITAB E

Sally T. Peavy, Ruth N. Varner,

and

Shirley G. Bremer

Statistical Engineering Laboratory

Applied Mathematics Division

Institute for Basic Standards

National Bureau of Standards

Washington, D.C. 20234

J?

-

*"»f*u of J

NBS Technical Notes are designed to supplement the Bureau's regular publications program. They provide a means for making available scientific data that are of transient or limited interest. Technical Notes may be

listed or referred to in the open literature.

For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C, 20402.

(Order by SD Catalog No. C 13.46:550). Price 50 cents.

A Systems Programmer's Guide For Implementing OMNI TAB II

Sally T. Peavy, Ruth N. Vainer, and Shirley G. Bremer

OMNITAB II is a general-purpose program which permits direct use of

a computer without prior knowledge of computer languages . Every effort

has been made to produce a system as machine independent as possible to make implementation on any large computer configuration relatively easy.

However, there are a few modifications which may have to be made.

This Technical Note provides assistance to the systems programmer,

with the task of implementing OMNITAB II, by pointing out where

difficulties may occur and how to cope with them. It furthermore

outlines a method for segmenting the OMNITAB II system which is very

large. It is a partial documentation of the OMNITAB program.

OMNITAB II is a large system requiring a large computer. Overlay

and segmentation are virtually essential. A method for segmenting

OMNITAB II is outlined. The method should be useful for many computers.

Keywords: ANSI FORTRAN, double precision, general-purpose computer

program, implementation of OMNITAB II, labeled common, machine

independent, OMNITAB II, overlay, segmentation, system parameters,

transportable computer programs.

OMNITAB II is a general-purpose interpretive program which permits direct use of a computer without a prior knowledge of computer languages. The OMNITAB II system allows

users to utilize a computer to solve their problems in an effective and accurate manner.

OMNITAB was first conceived by Joseph Hilsenrath in the early 1960 's to permit

scientists, engineers and laboratory technicians the use of a large high-speed computer in

the solution of their problems without the assistance of computer programmers. The system

was developed and implemented further by his co-workers, Hilsenrath et al. (1966), for a specific computer and therefore contained subprograms written in symbolic machine language

and FORTRAN. With the advent of third generation computers, it became necessary to rewrite

the OMNITAB program. Walter J. Gilbert, who undertook this task in 1966, initiated a number

of major changes and new ideas. Since 1968 the Statistical Engineering Laboratory of the

Applied Mathematics Division has been responsible for maintaining and developing OMNITAB.

OMNITAB II (Hogben et al. (1970)) is the result of this effort.

This paper is one of four which constitute the documentation for OMNITAB II. A complete

program listing is given in Peavy et al. (1970). A user's guide is given in Hogben et al.

(1970). Test problems with results are given in Varner et al. (1970). The material in this

note, without this introduction, is also stored on the magnetic tape which contains the

OMNITAB II master program and is available from National Technical Information Service

(formerly Clearinghouse), U.S. Department of Commerce, Springfield, Va. 22151. No attempt

is made here to carefully define terms, such as worksheet, which are unique with OMNITAB.

It is assumed the reader is familar with Hogben et al. (1970).

Extensive effort has been exerted to produce virtually machine independent subprograms

in order to make OMNITAB II transportable. The OMNITAB II program has been rewritten in the
American National Standards Institute (ANSI) FORTRAN language. We have tried to avoid the use of any ANSI FORTRAN statement which can not be successfully compiled on any specific

computer. All output is in 120 character per line format. The maximum number of

alphanumeric characters per variable has been limited to three. In some instances, the

steps taken have resulted in a loss of efficiency. For example, use of the UNIVAC FORTRAN V

function FID would have simplified the programming of the subroutine RFORMT. However,

transportability was considered far more important than machine efficiency and a few

compromises had to be made.

Special attention has been given to make it easy to change pertinent system parameters

for implementation on computers with different memory size, word length and logical inputoutput units. All input-output FORTRAN statements use an integer variable in referencing an

I/O unit rather than a specific logical unit. These variables are defined in one subprogram

(SETUP) and only this subroutine needs modification if the present assigned values are not

compatible with the user's particular computer configuration. There are a few instances

where alphanumeric characters must be packed to the full capacity of the computer word

length (six characters per word are assumed) . Again, the systems programmer needs to modify

only one subprogram (PREPAK) , see XREF3, section 2. These features have enabled systems

programmers to successfully implement the OMNITAB II system on many different large

computers (e.g. IBM 360/50 up, GE 625, CDC 3800 and 6600, Burroughs 5500 and UNIVAC 1108)

with a minimum of effort.

The OMNITAB II system is large and, therefore, segmentation is necessary. XREF3,

section 3, contains an outline of the overlay plan used with the National Bureau of

Standards computer. Careful attention was given to the segmentation in order to minimize

the flip-flopping of segments in and out of the computer. Statistics shaving the most

frequently used commands, which were obtained over a three month period, were employed in

determining what subprograms were to be in residence at all times. Furthermore, subprograms

which execute related instructions, such as those for matrix and array operations, were

grouped together. In some cases, in order to solve the overlay problem, it was necessary to duplicate a subroutine in another segment under a different subprogram name. This was avoided as much as possible and was kept to a minimum. Overlay and/or segmentation

techniques vary greatly with different computer software systems. In some cases, systems

programmers may have to expend considerable effort to segment OMNITAB on a particular

computer, a problem which will be investigated soon. If the computer has a large memory

storage and operation is in the batch mode, overlay may not be required. The present size

of OMNITAB II, assuming a floating-point number occupies one memory location, is

approximately 90K words.

Needless to say, such a large scale general-purpose program necessitates extensive

documentation for the user as well as for the systems programmer implementing OMNITAB II. The purpose of the information presented in this publication is to assist in the

implementation procedure. The documentation is in the form of FORTRAN comment statements.

The last eight columns of each card contain the subprogram name, e.g. XR1, and the line

number, e.g. 210, as in Peavy et al. (1970). This permits easy referencing and updating.

This outline should help to solve most of the problems which arise in implementing OMNITAB

II on a computer. This documentation was prepared as a preliminary guide rather than a definitive manual. A systems programmer's manual is being considered to complement this

guide. Readers who experience any difficulties in the implementation of OMNITAB are asked

to report their results to us so that improvements may be made and passed on to others.

OMNITAB II, Version 5.0 was designed for use in the batch processing mode. However, the

very nature of OMNITAB makes it adaptable for remote batch processing or time-sharing.

Walter J. Gilbert adapted an earlier version of OMNITAB for time-sharing which is being used

at the University of Maryland and the University of Rome, Italy. One of the changes that

has to be made in a time-sharing version is in the subroutine ERROR, see Peavy et al.

(1970). The statement on line ERR 210 needs to be changed so that NERROR=0 instead of

NERR0R=NERR0R+1. If a terminal having only 72 characters per line is used, major changes

would have to be made in the subprograms which execute the instructions which have a comprehensive automatic printing. See Hogben et al. (1970) for a listing. Also, a change
would have to be made in the command PLOT, although the command PAGE PLOT could be modified

and used to replace PLOT. When the necessary hardware and software are added to the NBS

computer configuration, a timesharing version of CMNITAB will be developed.

We thank David Hogben for his valuable assistance and suggestions.

The documentation is grouped as follows

:

XREF Page 5

This gives an outline of the information contained in XREF1, XREF2, XREF3, XREF4, XREF5

and XREF6.

XREF1 Page 6

A complete list of all the labeled common statements and the subprograms in which they

appear is given. A list of subprograms and the referencing subprograms is provided. Also

included are the system library functions (e.g. SIN, COS, etc.) used in the OMNITAB II system and the subprograms which reference them. Enumerated are the subprograms which

contain DOUBLE PRECISION statements along with the variables typed as double precision.

EQUIVALENCE statements are used in a number of subprograms. The variables and their

equivalent variables are listed along with the subprograms in which they appear. Hence, it

should be easier to resolve any difficulties caused by the use of (or changes in the use of)

DOUBLE PRECISION or EQUIVALENCE statements.

XREF2 Page 20

Each CMNITAB command is assigned two unique values for internal use. As the OMNITAB

instruction is scanned, the values for that particular command are stored in the variables

LI and L2. These variables are used as switches by the subprogram XECUTE and others to control the flow of the program in executing the instructions. XREF2 contains a list of all

the command names, appropriate values for LI and L2, and the name of the subprogram used to execute the command.

XREF

3

Page 24

The sizes of the memory spaces allocated in the OMNITAB II system for the worksheet and

scratch areas are 12,500 and 13,500 words respectively. It may be necessary to decrease the

size of these areas if a smaller computer is employed or segmentation is not used. In other

instances, due to large amounts of input data, larger worksheet and scratch areas may be

needed. Because of the limitations of ANSI FORTRAN, a number of changes are necessary if

the worksheet and scratch areas are to be redimensioned. A complete set of instructions is

included in this section which spell out the labeled common statements that must be

modified, the subprograms which contain these labeled common statements, and a number of

parameters that must be redefined in the subprogram SETUP.

The OMNITAB II program allows the user to provide headings for any fifty columns. Also,

formats may be specified for user controlled input-output. If the number of characters per

word is other than six, modifications are necessary before the FORMAT command can be used.

Instructions for handling formats and for changing the number of headings permitted are

given.

An outline of the segmentation used at NBS and the number of memory locations needed by

each segment is described. There are fourteen segments, the last two containing

subsegments. The segments are labeled as PARTI, PART2, ... , PART14 and contain the names of all the subprograms and labeled common areas of each segment. The first segment is

resident in the computer at all times while only one of the other segments, PART2 through

PART14, is in memory storage at a particular moment. The size of each segment is an approximation and dependent on the fact that a floating-point number occupies one memory
location. The outline given may not necessarily be the ideal structure for overlay on other

computer configurations. The sizes given are based on the use of a 36 binary bit word.

Four physical input-output units are used by CMNITAB II. These units are pre-assigned

logical units in the subprogram SETUP which is called only once when CMNITAB first enters

the operating system. The logical units may be reassigned, if necessary, for different

installations.

XREF4 Page 34

Bounds are imposed by CMNITAB on the subprograms evaluating trigonometric, exponential,

square root, and logarithmic functions. These bounds are defined in DATA statements that

may need modification due to a different word length of other computers. Furthermore, the

section indicates variables which must be altered if the number of significant digits for a real number is not eight or if the largest integer variable allowed is less than 2 to the

52nd power.

XREF5 Page 37

Four of the subprograms in the CMNITAB II system are BLOCK DATA subroutines. This

section lists the variables and the values assigned to the variables by the DATA statements

in the BLOCK DATA subprograms. Also included here, is the character set recognized and used

by the CMNITAB II system.

XREF6 Page 39

As CMNITAB instructions are scanned by the program, a thorough check is made for errors

and messages are printed whenever an error is detected. Messages are printed for fatal

errors, arithmetic faults and informative diagnostics. A list of all the messages and their

corresponding number which is used when a subprogram calls the ERROR routine are included.

See also Hogben et al. (1970) for further details.

References

AMERICAN STANDARD FORTRAN (1966). American National Standards Institute, New York.

HILSENRATH, J., ZIEGLER, C. G. , MESSINA, C. G. , WALSH, P. J. and HERBOLD, R. J. (1966).

CMNITAB: A computer Program For Statistical And Numerical Analysis. National Bureau Of

Standards Handbook 101, Superintendent of Documents, U.S. Government Printing Office,

Washington, D.C. 20402. Reissued January 1968 with corrections.

HOGBEN, DAVID, PEAVY, S. T. and VARNER, R. N. (1970). CMNITAB II User's Reference Manual.

NBS Technical Note 552. Superintendent of Documents, U.S. Government Printing Office,

Washington, D.C. 20402.

PEAVY, S. T., VARNER, R. N. and HOGBEN, DAVID (1970). Source Listing Of CMNITAB II Program.

NBS Special Publication 339. Superintendent of Documents, U.S. Government Printing

Office, Washington, D.C. 20402.

VARNER, R. N. and PEAVY, S. T. (1970). Test Problems And Results For CMNITAB II. NBS

Technical Note 551. Superintendent of Documents, U.S. Government Printing Office,

Washington, D.C. 20402.
Listing Of XREF

C XREF XRF 10 C VERSION 5.00 XREF 5/15/70 XRF 20 C XRF 30 C COMMENTS AND DOCUMENTATION XRF 40 C XRF 50 C THESE -PROGRAMS -CONTAIN NO EXECUTABLE STATEMENTS. XRF 60

C XRF 70 C THE DOCUMENTATION IS GROUPED AS FOLLOWS XRF 80 C XRF 90 C XREF1 XRF 100

C 1. CROSS REFERENCE TABLE XRF 110

C A. LIST OF ALL LABELED COMMON XRF 120 C B. CROSS REFERENCE BETWEEN LABELED COMMON AND SUBPROGRAM XRF 130

C C. CROSS REFERENCE BETWEEN SUBPROGRAMS XRF 140

C 2. SUBPROGRAMS WITH DOUBLE PRECISION VARIABLES XRF 150

C 3. EQUIVALENCE VARIABLES IN SUBPROGRAMS XRF 160

C XRF 170

C XREF2 XRF 180

C 1. COMMAND NAMES AND THE SUBPROGRAM THAT PERFORMS THE XRF 190

C INSTRUCTION XRF 195

C XRF 200

C XREF3 XRF 210 C 1. INCREASING OR DECREASING WORKSHEET XRF 220 C 2. CHANGING NO. OF HEADINGS PERMITTED AND PACKING FORMATS XRF 230

C 3. OVERLAY OR SEGMENTATION XRF 240 C 4. ASSIGNMENT OF INPUT-OUTPUT UNITS XRF 250 C 5. DIRECTIONS FOR THE USE OF DUMMY A-F XRF 260

C XRF 290

C XREF4 XRF 300

C 1. CHANGES NEEDED FOR SYSTEM FUNCTIONS XRF 310

C 2. CHANGES NEEDED IN BLOCK DATA SUBPROGRAM XRF 320

C WHICH DEFINES PHYSICAL CONSTANTS XRF 330

C 3. CHANGES NEEDED IF NUMBER OF SIGNIFICANT DIGITS IN FLOATING XRF 340

C POINT IS NOT EQUAL TO 8 XRF 345

C 4. CHANGES NEEDED IF INTEGER NUMBERS ARE LESS THAN 2**32 XRF 350

C XRF 355 C XREF5 XRF 360

C 1. DESCRIPTION OF VARIABLES DEFINED IN BLOCK DATA SUBPROGRAMS XRF 370 C A. BLOCK XRF 380

C B. LBCONS XRF 390

C C. LOOKTB XRF 400

C D. PHYSIC XRF 410

C XRF 415

C XREF6 XRF 420

C 1. LISTING OF ALL OMNITAB ERROR MESSAGES XRF 430

C A. FATAL ERROR MESSAGES XRF 440

C B. ARITHMETIC ERROR MESSAGES XRF 450

C C. INFORMATIVE DIAGNOSTIC MESSAGES XRF 460

C XRP 470

END XRF 480
Listing Of XREF1

C XREF1 XR1 10

C VERSION 5.00 XREF1 5/15/70 XR1 20 C XR1 30 C THIS IS A CROSS-REFERENCE TABLE SHOWING WHICH SUBPROGRAMS XR1 40 C REFERENCE PARTICULAR BLOCKS OF COMMON AND PARTICULAR SUBPROGRAMS. XR1 50 C XR1 80 C OMNITAB USES NO BLANK COMMON. XR1 90

C XR1 100

C XR1 110

C XR1 120

C XR1 130

c 1>A ********mis IS LABELED COMMON AS USED IN CMNITAB*********** XR1 140

C XR1 150

C COMMON /ABCDEF/ L(48) XR1 160

C COMMON/BLOCKAAlODE,M,KARD(83),KARG,ARG,ARG2,NEWCD(80),KRDEND XR1 170

C COMMON/BLOCKB/NSTMT,NSTMTX,NSTMTH,NCOM,LCOM,IOVFL,COM(2000) XR1 180

C COMMON/BLOCKC/KIO,INUNIT,ISCRAT,KBDOUT,KRDKNT,LLIST XR1 190 C COMMON/BLOCKD/ IARGS(IOO) ,KIND(100) ,ARGTAB(100) ,NRMAX, XR1 200

C 1 NROW,NCOL,NARGS,VWXYZ(8),NERROR XR1 210

C COMMON/BLOCKE/ NAME(4),L1,L2,ISRFLG XR1 220 C COMMON/BLOCRC/NRC,RC( 12600) XR1 230

C DIMENSION ARGS(100) XR1 240

C EQUIVALENCE (ARGS(l) ,RC(12501)

)

XR1 250 C COMMON/BLOCKX/ INDEX (6, 8) , LEVEL XR1 260

C COMMON/CODE/IALPH(6),NALPH(5),ID(9,3),IR(300,4),IRD(30,6) XR1 270

C COMMON/CONLB2/ER,ISIGD XR1 275

C COMMON/CONSLB/XTRIG,XEXP XR1 280

C COMMON/CONSTS/PI,E,HALFPI,DEG,RAD,XALOG XR1 290

C COMMON/DCONL2/TRRTPI,NBC,NBM XR1 295

C CCMMON/DCONLB/DSNCOS.DXEXP XR1 300

C DOUBLE PRECISION DSNCOS,DXEXP XR1 310 C COMMON/FMAT/IFMTX( 6) ,IOSWT,IFMTS( 6) ,LHEAD(96) XR1 320 C COMMON/HEADER/NOCARD(80),ITLE(60,6) ,LNCNT,IPRINT,NPAGE,IPUNCH XR1 330 C COMMON/ICODE/NIR,NID,NIRD,LIR,LID,LIRD XR1 340

C COMMON/KFMT/KFMT(100) XR1 350

C CCMMQN/NOTE/NOTEC120) XR1 355

C COMMON/PCONST/JPC,P(40),N(40) XR1 360

C COMMON/PKSWT/IHCNT,IHTP XR1 370

C CCMMON/QRS/NDROW,IFLAG,J,NNARG XR1 380 C COMMON/SCRAT/NS,NS2,A( 13500) XR1 390

C COMMON/SPRV/NERCON,NERR,ISWERR XR1 400

C***THIS LABELED COMMON IS USED ONLY BY THE TAPE OPERATION SUBPROGRAMS* XR1 410

C XR1 420 C COMMON/CODETP/ITP(10,4) XR1 430

C COMMON/ICODTP/NITP,LITP XR1 440

C C0NW0N/TAPE/NAME4(2),NTPCT,IPUNCP,INUNIP,L1TP XR1 450 C XR1 460

C************THIS LABELED COMMON IS USED ONLY BY BESSEL SUBPROGRAMS**** XR1 470 C XR1 480 C CONMON/ABEKI/X,Y,P,Q,S,T XR1 490 C DOUBLE PRECISION X,Y,P,Q,S,T XR1 500 C XR1 510
c l.B * ** CROSS REFERENCE BETWEEN LABELED COMMON AND SUBPROGRAMS*** XR1 515

c XR1 520

c LABELED COMMON SUBP)'i!<;i:W!N XR1 525

c XR1 530

c XR1 535

c ABCDEF XR1 540 c APRINT BLOCK FIXFLO HEADS HISTGM XR1 545

c OMCONV ONEWAY ORTHO ORTPLT PACK PREPAK XR1 550 c RFORMT RPRINT TAPOP TWOWAY XHEAD XCMNIl XR1 560 c ABEKI XR1 570 c BESSEL CBEI CBEK XR1 580

c BLOCKA XR1 590 c AARGS ASTER BEGIN INPUT NNAME NONBLA XR1 600

c NOTEPR OMNIT OUTPUT PHYCON PREPAK READQ XR1 610

c READX SET XR1 620

c SETQ SETUP STMT STORE TAPOP VARCON XR1 630

c XFORMT XHEAD XOMNIT XR1 640

c BLOCKB XR1 650

c BEGIN INPUT LOCATE OMNIT OUTPUT REPINC XR1 660

c STORE XOMNIT XR1 670

c BLOCKC XR1 680

c AERR ERROR INFERR INPUT INTERP INVERT XR1 690

c LIST OMNIT OUTPUT PREPAK READQ READX XR1 700

c RNDOWN SET SETQ SETUP TAPOP2 XOMNIT XR1 710

c XSTOP XR1 720

c BLOCKD XR1 730

c ABRIDG ADRESS ALLSUB APRINT ARITH ARYVEC XR1 740

c BEGIN BESSEL CHANGE CHKCOL CKIND CMSEPA XR1 750 c COALES COMPLX CORREL DEFINE DIMENS XR1 760 c ERASE ERROR EXCHNG XR1 770

c EXPAND EXPCON EXTREM FIXFLO FLIP XR1 780

c FNEC FNEIC FNKC FPROB XR1 790 c FRDIST FUNCT GENER GQUAD HISTGM XR1 800

c IFS INTERP INVERT ITERAT LIST MATRIX XR1 810

c MDAMAD MEIGEN MISC2 MIST XR1 820

c MKRON MMULT MOP XR1 830

c MOVE MPROP MRAISE MSCROW MTRIAN MTXCHK XR1 840

c MXTX QANOVA OCOEFF OCOVAR OMNIT ONEWAY XR1 850

c OPONE ORTHO ORTPLT OUTPUT PDMOTE PLOT XR1 860 c PRINTX PROROW PUNCH RANKS READQ READX XR1 870

c REPINC RESET RPRINT XR1 880 c SELECT SET SETQ SETUP SORDER XR1 890 c SPACE STATIS STORE TAPOP2 THERMO TRANSF XR1 900

c TWOWAY VECTOR XHEAD XOMNIT XR1 910

c XPND XSTOP XR1 920

c BLOCKE XR1 930

c ABRIDG ALLSUB APRINT ARITH ARYVEC BEGIN XR1 940

c BESSEL CMSEPA COALES COMPLX CORREL XR1 950

c EXPAND EXPCON EXTREM FIXFLO FNEC XR1 960

c FNEIC FNKC FRDIST FUNCT HISTGM IFS XR1 970

c INVERT ITERAT LOOKUP MATRIX MDAMAD MISC2 XR1 980

c MOP MPROP MSCROW MXTX OANOVA OCOEFF XR1 990

c OCOVAR OMNIT ONEWAY OPONE ORTHO XR1 1000

c PDMOTE PLOT PRINTX PROROW PUNCH RANKS XR1 1010

c READX REPINC RESET RPRINT SELECT SET XR1 1020

c SORDER STATIS STORE TAPOP TAPOP

2

THERMO XR1 1030

c TRANSF XECUTE XR1 1040

c BLOCRC XR1 1050

c ABRIDG ADRESS ALLSUB APRINT ARITH ARYVEC XR1 1060

c BEGIN BESSEL CHANGE CHKCOL CKIND CMSEPA XR1 1070
c COALES COMPLX CORREL DEFINE DIMENS ERASE

c ERROR EXCHNG EXPAND EXPCON EXTREM FIXFLO

c FLIP FNEC FNEIC FNKC FPROB

c FRDIST FUNCT GENER GQUAD HISTGM

c IFS INTERP INVERT ITERAT LIST MATRIX

c MDAMAD ME I GEN MISC2 MKRON MMULT MOP

c MOVE MPROP MRAISE MSCROW MTRIAN

c MXTX OMNIT ONBVAY OPONE

c ORTHO ORTPLT OUTPUT PDMOTE PLOT

c PRINTX PROROW PUNCH RANKS READQ READX

c REPINC RESET RPRINT

c SELECT SET SETQ SETUP SORDER

c SPACE STATIS STORE THERMO TRANSF TWOWAY

c VECTOR XOMNIT XPND XSTOP

c BLOCKX

c AERR ERROR IFS REPINC RNDOWN SETUP

c XECUTE XOMNIT

c CODE

c LOOKTB LOOKUP TAPOP

c CODETP

c LOOKTB LOOKUP TAPOP

c CONLB2

c CSPINV FIXFLO INFERR LBCONS

c CONSLB

c BESSEL FCOS FEXP FSIN FTANH LBCONS

c CONSTS

c AARGS COMPLX FUNCT LBCONS ORTPLT SETUP

c DCONL2

c ERRINT LBCONS

c DCONLB

c BESSEL DBEJ FDCOS FDEXP FDSIN LBCONS

c FMAT

c ABRIDG APRINT BLOCK FIXFLO HEADS

c OPONE ORTHO ORTPLT PLOT

c PRINTX PUNCH RPRINT XOMNIT

c HEADER

c ABRIDG APRINT BLOCK CORREL HEADS HISTGM

c MIST MPROP NOTEPR OANOVA OCOEFF OCOVAR

c OMNIT ONEWAY OPONE ORTHO ORTPLT

c PAGE PLOT PRINTX PUNCH RPRINT SETUP

c SPACE STATIS TAPOP2 TWOWAY XOMNIT XSTOP

c I CODE

c LOOKTB LOOKUP SETUP

c ICODTP

c LOOKTB LOOKUP SETUP

c KBIT

c ABRIDG APRINT OPONE ORTHO

c PRINTX PUNCH READQ READX TWOWAY

c NOTE

c NOTEPR

c PCONST

c PHYCON PHYSIC SETUP XOMNIT

c PKSWT

c PREPAK SETUP XOMNIT

c QRS

c READQ READX SET SETQ

c SCRAT

c ALLSUB APRINT ARYVEC BESSEL CBEK CMSEPA

c COALES COMPLX CORREL DBEJ EXPCON

c FNEC FNEIC FNKC FUNCT

XR1 1080

XR1 1090

XR1 1100

XR1 1110

XR1 1120

XR1 1130

XR1 1140

XR1 1150

XR1 1160

XR1 1170

XR1 1180

XR1 1190

XR1 1200

XR1 1210

XR1 1230

XR1 1240

XR1 1250

XR1 1260

XR1 1270

XR1 1280

XR1 1290

XR1 1300

XR1 1310

XR1 1320

XR1 1330

XR1 1340

XR1 1350

XR1 1360

XR1 1370

XR1 1380

XR1 1390

XR1 1400

XR1 1410

XR1 1420

XR1 1430

XR1 1440

XR1 1450

XR1 1460

XR1 1470

XR1 1480

XR1 1490

XR1 1500

XR1 1510

XR1 1520

XR1 1530

XR1 1540

XR1 1550

XR1 1555

XR1 1560

XR1 1565

XR1 1570

XR1 1580

XR1 1590

XR1 1600

XR1 1610

XR1 1620

XR1 1630

XR1 1640

XR1 1650

XR1 1660
C HISTGM INTERP INVERT ITERAT XR1 1670

C MATRIX MDAMAD MEIGEN MISC2 MKRON XR1 1680

C MMULT MOP MFROP MRAISE MTRIAN MXTX XR1 1690

C OANOVA OCOEFF OCOVAR ONEWAY OPONE ORTHO XR1 1700

C ORTPLT PLOT PREPAK PROROW RANKS XR1 1710

C RPRINT SELECT SETUP SORDER SORTSM STATIS XR1 1720

C TAPOP2 THERMO TRANSF TWOWAY XSTOP XR1 1730

C SPRV XR1 1740

C AERR ERROR XOMNIT XR1 1750

C TAPE XR1 1760

C LOOKUP OMNIT PUNCH READQ READX XR1 1770

C SET SETQ SETUP TAPOP TAPOP2 XR1 1780

C XR1 1790

C XR1 1800

C XR1 1810

C XR1 1820

C XR1 1830

pAAAAA***************************************************************** XR1 1840

C XR1 1850

c i m Q *************CR0SS REFERENCE BETWEEN SUBPROGRAMS************ XR1 1860

C XR1 1870

C SUBPROGRAM REFERENCING SUBPROGRAMS XR1 1890

C XR1 1900

C AARGS XR1 1920

C ASTER OMNIT XHEAD XR1 1930

C ABRIDG XR1 1940

C XECUTE XR1 1950

C ACCDIG XR1 1960

C ARITH XR1 1970

C ADRESS XR1 1980

C ARITH ALLSUB BESSEL CHANGE CHKCOL XR1 1985

C COMPLX DEFINE EXCHNG FNEC XR1 1990

C FNEIC FNKC FPROB FRDIST FUNCT XR1 2000

C GENER GQUAD HISTGM IFS INVERT MEIGEN XR1 2010

C MISC2 MOP MOVE MPROP MSCROW ORTHO XR1 2020

C PROROW RANKS READX RPRINT SELECT SET XR1 20 30 C THERMO TWOWAY XPND XR1 2040

C AERR XR1 2050

C ERROR OMNIT SETUP XECUTE XR1 2060

C ALLSUB XR1 2070

C XECUTE XR1 2080

C APRINT XR1 2090

C XECUTE XR1 2100

C ARITH XR1 2110

C XECUTE XR1 2120

C ARYVEC XR1 2130

C XECUTE XR1 2140

C ASTER XR1 2150

C OMNIT XR1 2160

C BEGIN XR1 2170

C XECUTE XR1 2180

C BEJN XR1 2190

C BESSEL BINTJ0 STRUVE XR1 2200

C BESSEL XR1 2210

C XECUTE XR1 2220

C BEZERO XR1 2230

C BESSEL XR1 2240

C BEZONE XR1 2250

C BESSEL XR1 2260

C BINTJ0 XR1 2270
c BESSEL

c BJORCK

c CORREL

c CBEI

c BESSEL

c CBEK

c BESSEL

c CHANGE

c XECUTE

c CHKCOL

c ABRIDG CMSEPA

c CORREL ERASE EXTREM FLIP INTERP ITERAT

c MISC2 ONEWAY PDMOTE PLOT PRINTX PROROW

c PUNCH READX RPRINT SELECT SORDER STATIS

c TIVOWAY c CKIND

c APRINT ARYVEC CMSEPA

c COALES EXPCON INTERP

c INVERT MATRIX MDAMAD MEIGEN MISC2 MKRON

c MMULT MOP MPROP MRAISE MXTX

c STATIS TAPOP2 THERMO TRANSF

c CMPARA

c CMSEPA

c CMSEPA

c XECUTE

c COALES

c XECUTE

c COMELL

c BESSEL

c CQMPLX

c XECUTE

c CORREL

c XECUTE

c CSPINV

c INVCOR

c DBEJ

c BESSEL BINTJO STRUVE

c DEFINE

c XECUTE

c DETRNK

c MPROP

c DHRND

c FNEIC

c DIMENS

c XECUTE

c DIM A-F

c XECUTE

c ERASE

c XECUTE

c ERRINT

c FNEC

c ERROR

c AARGS ABRIDG ALLSUB APRINT ARITH ARYVEC

c BEGIN BEJN BESSEL CBEK CHANGE CMSEPA

c COALES CONELL COMPLX CORREL DEFINE DIMENS

c ERASE EXCHNG EXPAND EXPCON EXTREM

c FCOS FDCOS FDEXP FDLOG FDSIN FDSQRT

c FEXP FIXFLO FLIP FLOG FLOG10 FNEC

c FNEIC FNKC FPROB FRDIST FSIN FSQRT

c FUNCT GENER GQUAD HISTGM

XR1 2280

XR1 2290

XR1 2300

XR1 2310

XR1 2320

XR1 2330

XR1 2340

XR1 2350

XR1 2360

XR1 2370

XR1 2380

XR1 2390

XR1 2400

XR1 2410

XR1 2420

XR1 2430

XR1 2440

XR1 2450

XR1 2460

XR1 2470

XR1 2480

XR1 2490

XR1 2500

XR1 2510

XR1 2520

XR1 2530

XR1 2540

XR1 2550

XR1 2560

XR1 2570

XR1 2580

XR1 2590

XR1 2600

XR1 2610

XR1 2620

XR1 2630

XR1 2640

XR1 2650

XR1 2660

XR1 2670

XR1 2680

XR1 2690

XR1 2700

XR1 2710

XR1 2720

XR1 2730

XR1 2740

XR1 2750

XR1 2760

XR1 2770

XR1 2780

XR1 2790

XR1 2795

XR1 2800

XR1 2810

XR1 2820

XR1 2830

XR1 2840

XR1 2850

XR1 2860

10
C IFS INTERP INVERT HERAT MATRIX

C MDAMAD MEIGEN MISC2 MKRON MMULT MOP

C MOVE MPROP MRAISE MSCROW MTRIAN MXTX

C OMNIT ONHVAY ORTHO PDMOTE PLOT PREPAK

C PRINTX PROB PROROW PUNCH RANKS

C READQ READX REPINC RESET RPRINT

C SELECT SET SETQ SORDER SPACE

C STATIS STORE TAPOP2 THERMO

C TRANSF TWOWAY XFORMT XHEAD

C EXCHNG

C XECUTE

C EXPAND

C OMNIT REPINC

C EXPCON

C XECUTE

C EXTREM

C XECUTE

C FCOS

C FUNCT

C FDCOS

C BESSEL CBEI CBEK COMPLX DBEJ FOURIA

C PROB

C FDEXP

C BESSEL CBEI CBEK DBEJ ERRINT PROB

C THERMO

C FDLOG

C CBEK COMELL DBEJ PROB THERMO

C FDPCON

C COMPLX DHRND FNEC GQUAD

C FDSIN

C CBEI CBEK COMPLX DBEJ FOURIA PROB

C FDSQRT

C BEZERO BEZONE CBEI CBEK COMELL

C COMPLX DBEJ ORTHO TWOWAY

C FEXP

C FEXP2 FPPT FUNCT

C FEXP2

C AARGS ARITH MATRIX MISC2

C FIXFLO

C XECUTE

C FLIP

C XECUTE

C FLOG

C CORREL FEXP2 FUNCT ONHVAY THERMO

C FLOG10

C ACCDIG DHRND FREQCY FUNCT ORTHO RFORMT

C RPRINT

C FNEC

C XECUTE

C FNEIC

C XECUTE

C FNKC

C XECUTE

C FOURIA

C BESSEL

C FPPT

C ONHVAY

C FPROB

C XECUTE

C FRDIST

11

XR1 2880

XR1 2890

XR1 2900

XR1 2910

XR1 2920

XR1 2930

XR1 2940

XR1 2950

XR1 2960

XR1 2970

XR1 2980

XR1 2990

XR1 3000

XR1 3010

XR1 3020

XR1 3030

XR1 3040

XR1 3050

XR1 3060

XR1 3070

XR1 3080

XR1 3090

XR1 3100

XR1 3110

XR1 3120

XR1 3130

XR1 3140

XR1 3150

XR1 3160

XR1 3165

XR1 3170

XR1 3180

XR1 3190

XR1 3200

XR1 3210

XR1 3220

XR1 3230

XR1 3240

XR1 3250

XR1 3260

XR1 3270

XR1 3280

XR1 3290

XR1 3300

XR1 3310

XR1 3320

XR1 3330

XR1 3340

XR1 3350

XR1 3360

XR1 3370

XR1 3380

XR1 3390

XR1 3400

XR1 3410

XR1 3420

XR1 3430

XR1 3440

XR1 3450

XR1 3460
c

cc

c

cc

c

c

c

cc

c

c

cc

c

cc

c

c

cc

c

cc

cc

cc

c

cc

cc

c

c

cc

c

c

c

c

cc

c

c

cc

c

c

c

c

ccc

c

c

c

c

c

FREQCY

FSIX

FSQRT

FTANH

FUNCT

GENER

GQUAD

HDIAG

HEADS

HISTGM

IFS

IXFERR

INPUT

IXTERP

INTRP

INVCHK

INVCOR

INVERT

ITERAT

LIST

LOCATE

LOOKUP

MATRIX

mPMAD

MEIGEN

MISC2

MIST

MKRON

NMULT

XECUTE

FRDIST

FUNCT

BJORCK CORREL FPPT FUNCT HDIAG IXTERP

INVCHK INVCOR MSCROW MTRIAN OCOEFF ONBVAY

OPONE ORTHO ORTPLT PROB STATIS TWOWAY

CORREL FUNCT

XECUTE

XECUTE

XECUTE

MEIGEN

PLOT PRINTX RPRINT

XECUTE

XECUTE

ERROR

OMNIT

XECUTE

IXTERP

INVERT MPROP

CORREL

XECUTE

XECUTE

XECUTE

REPINC STORE

OMNIT

XECUTE

XECUTE

XECUTE

XECUTE

CORREL

XECUTE

XR1 3470

XR1 3480

XR1 3490

XR1 3500

XR1 3510

XR1 3520

XR1 3530

XR1 3540

XR1 3550

XR1 3560

XR1 3570

XR1 3580

XR1 3590

XR1 3600

XR1 3610

XR1 3620

XR1 3630

XR1 3640

XR1 3650

XR1 3660

XR1 3670

XR1 3680

XR1 3690

XR1 3700

XR1 3710

XR1 3720

XR1 3730

XR1 3740

XR1 3750

XR1 3760

XR1 3770

XR1 3780

XR1 3790

XR1 3800

XR1 3810

XR1 3820

XR1 3830

XR1 3840

XR1 3850

XR1 3860

XR1 3870

XR1 3880

XR1 3890

XR1 3900

XR1 3910

XR1 39 20 XR1 3930

XR1 3940

XR1 3950

XR1 3960

XR1 39 70 XR1 3980

XR1 3990

XR1 4000

XR1 4010

XR1 4030

XR1 4040

XR1 4050

XR1 4060

XR1 4070

12
C XECUTE XR1 4080

C MOP XR1 4090

C XECUTE XR1 4100

C MOVE XR1 4110

C MISC2 XECUTE XR1 4120

C MPROP XR1 4130

C XECUTE XR1 4140

C MRAISE XR1 4150

C XECUTE XR1 4160

C MSCROW XR1 4170

C XECUTE XR1 4180

C MTRIAN XR1 4190

C XECUTE XR1 4200

C MTXCHK XR1 4210

C APRINT ARYVEC EXPCON INVERT MATRIX MDAMAD XR1 4220

C ME I GEN MKRON MMULT MOP MPROP MRAISE XR1 4230

C MTRIAN MXTX TRANSF XR1 4240

C MXTX XR1 4250

C XECUTE XR1 4260

C MXTXP XR1 4270

C MXTX ORTHRV XR1 4280

C NNAME XR1 4290

C ASTER OMNIT TAPOP XR1 4300

C NONBLA XR1 4310

C ASTER XR1 4320

C NOTEPR XR1 4330

C OMNIT XECUTE XOMNIT XR1 4340

C OANOVA XR1 4350

C ORTHO XR1 4360

C OCOEFF XR1 4370

C ORTHO XR1 4380

C OCOVAR XR1 4390

C ORTHO XR1 4400

C OMCONV XR1 4410

C INPUT XR1 4420

C OMNIT XR1 4430

C MNITAB OMNSYM XR1 4440

C ONEWAY XR1 4450

C XECUTE XR1 4460

C OPONE XR1 4470

C ORTHO XR1 4480

C ORTHO XR1 4490

C XECUTE XR1 4500

C ORTHRV XR1 4510

C MPROP PROCHK XR1 4520

C ORTPLT XR1 4530

C ORTHO XR1 4540

C OUTPUT XR1 4550

C OMNIT XR1 4560

C PACK XR1 4570

C PREPAK XR1 4580

C PAGE XR1 4590

C ABRIDG APRINT CORREL HISTGM MIST MPROP XR1 4600

C NOTEPR OMNIT ONEWAY ORTHO ORTPLT PLOT XR1 4610

C PRINTX RPRINT STATIS TWOWAY XECUTE XSTOP XR1 4620

C PDMOTE XR1 4630

C XECUTE XR1 4640

C PHYCON XR1 4670

C ASTER XECUTE XR1 4680

C PLOT XR1 4690

13
c XECUTE

c PREPAK

c ABRIDG APRINT HEADS ORTHO PRINTX PUNCH

c READX XFORMT XHEAD XOMNIT

c PRINTX

c XECUTE

c PROB

c CORREL FPPT FPROB OANOVA ONEWAY

c STATIS TWOWAY

c PROCHK

c MPROP

c PROROW

c XECUTE

c PUNCH

c TAPOP2 XECUTE

c PVTRI

c MPROP

c RANKO

c ONEWAY

c RANKS

c XECUTE

c RANKX

c CORREL RANKS

c RCSUM

c MPROP

c READQ

c OMNIT

c READX

c TAPOP2 XECUTE

c REPINC

c XECUTE

c RESET

c XECUTE

c RFORMT

c APRINT HISTGM OANOVA OCOEFF OCOVAR OPONE

c ORTHO RPRINT TWOWAY

c RNDOWN

c AERR ERROR

c RNJBK

c FNKC

c RPRINT

c ABRIDG PRINTX

c SELECT

c XECUTE

c SET

c TAPOP2 XECUTE

c SETQ

c OMNIT

c SETUP

c OMNIT

c SKSYMV

c PROCI IK c SORDER

c XECUTE

c SORTSM

c ARYVEC MMULT MRAISE MT'RIAN MXTXP TRANSF

c SPACE

c XECUTE

c SPINV

c INVQ1K

XR1 4700

XR1 4710

XR1 4720

XR1 4730

XR1 4740

XR1 4750

XR1 4760

XR1 4770

XR1 4780

XR1 4790

XR1 4800

XR1 4810

XR1 4820

XR1 4830

XR1 4840

XR1 4850

XR1 4860

XR1 4870

XR1 4880

XR1 4890

XR1 4900

XR1 4910

XR1 4920

XR1 4930

XR1 4940

XR1 4950

XR1 4960

XR1 4970

XR1 4980

XR1 4990

XR1 5000

XR1 5010

XR1 5020

XR1 5030

XR1 5040

XR1 5050

XR1 5060

XR1 5070

XR1 5080

XR1 5090

XR1 5100

XR1 5110

XR1 5120

XR1 5130

XR1 5140

XR1 5150

XR1 5160

XR1 5170

XR1 5180

XR1 5190

XR1 5200

XR1 5210

XR1 5220

XR1 5230

XR1 5240

XR1 5250

XR1 5260

XR1 5270

; XR1 5280

XR1 5290

14
c STATIS

c XECUTE

c STMT

c OMNIT

c STORE

c OMNIT

c STORMT

c MXTX

c STRUVE

c BESSEL

c SYMV

c MTRIAN PROCHK

c TAPOP

c OMNIT

c TAPOP2

c XECUTE

c THERMO

c XECUTE

c TPCTPT

c FPPT ONEWAY

c TRANSF

c MXTX

c TWOWAY

c XECUTE

c VARCON

c ASTER

c VECTOR

c DEFINE ERASE EXTREM FUNCT

c XECUTE

c OMNIT

c XFORMT

c OMNIT

c XHEAD

c OMNIT

c XOMNIT

c OMNIT

c XPND

c EXPAND REPINC

c XSTOP

c OMNIT XOMNIT

cc

***ft**********ft**CROss REFERENCING OF ! C

C FUNCTION REFERENCED BY

C

C *ALOG

C FLOG

c *ALOG10

c FLOG 10 c *ATAN

c FUNCT PROB

c *COS

c FCOS

c *DATAN2

c COMPLX

c *DCOS

c FDCOS

c *DEXP

c FDEXP

c *DLOG

XR1 5300

XR1 5310

XR1 5320

XR1 5330

XR1 5340

XR1 5350

XR1 5360

XR1 5370

XR1 5380

XR1 5390

XR1 5400

XR1 5410

XR1 5420

XR1 5430

XR1 5440

XR1 5450

XR1 5460

XR1 5470

XR1 5480

XR1 5490

XR1 5500

XR1 5510 XRI 5520

XR1 5530

XRI 5540

XRI 5550 XRI 5560

MISC2 MSCROW XRI 5570

XRI 5580

XRI 5590

XRI 5600 XRI 5610 XRI 5620

XRI 5630

XRI 5640 XRI 5650

XRI 5660 XRI 5670

XRI 5680

XRI 5690 XRI 5700

fSTEM FUNCTIONS**************** XRI 5705

XRI 5710

XRI 5715

XRI 5720

XRI 5730 XRI 5740 XRI 5750

XRI 5760 XRI 5770

XRI 5780

XRI 5790 XRI 5800

XRI 5810 XRI 5820

XRI 5830

XRI 5840

XRI 5850

XRI 5860 XRI 5870

15
c FDLOG XR1 5880

c *DSIN XR1 5890

c FDSIN XR1 5900

c *DSQRT XR1 5910

c FDSQRT XR1 5920

c *EXP XR1 5930

c FEXP XR1 5940

c *SIN XR1 5950

c FSIN XR1 5960

c *SQRT XR1 5970

c FSQRT XR1 5980

c *TANH XR1 5990

c FTANH XR1 6000

c XR1 6010

c 2 **********£)Qljg]j: PRECISIC DEFINED IN SUBPROGRAMS** XR1 6020

c XR1 6030

c SUBPROGRAM VARIABLE XR1 6040

c XR1 6050

c ARYVEC XR1 6060

c SUM X XR1 6070

c BEJN XR1 6080

c A B C D E F XR1 6090

c G P Q R X Y XR1 6100

c Z XR1 6105

c BESSEL XR1 6110

c AA B BINTJO COMELL DBEJ DSNCOS XR1 6120

c DXEX DXEXP E FDCOS FDEXP P XR1 6130

c Q S T W X XEX XR1 6140

c Y Z XR1 6145

c BEZERO XR1 6150

c A AA AB AC B FDSQRT XR1 6160

c X Y XR1 6165

c BEZONE XR1 6170

c A B FDSQRT R S T XR1 6180

c X Y XR1 6190

c BINTJO XR1 6200

c A B BINTJO C DBEJ X XR1 6210

c Z XR1 6220

c BJORCK XR1 6223

c c D FDSQRT R Y XR1 6225

c CBEI XR1 6230

c A AA B C D E XR1 6240

c F FDCOS FDEXP FDSIN FDSQRT G XR1 6250

c H P Q R S T XR1 6260

c U V W X Y Z XR1 6270

c CBEK XR1 6280

c A AA AB AC AD AE XR1 6290

c B C D E F FDCOS XR1 6300

c FDEXP FDLOG FDSIN FDSQRT G H XR1 6310

c P Q R S T U XR1 6320

c V W X Y Z XR1 6330

c COMELL XR1 6340

c A B C COMELL D E XR1 6350

c FDLOG FDSQRT P Q X Z XR1 6360

c COMPLX XR1 6370

c D FDCOS FDSIN FDSQRT X Y XR1 6380

c DBEJ XR1 6390

c A B C D DBEJ DSNCOS XR1 6400

c DXEXP E FDCOS FDEXP FDLOG FDSIN XR1 6410

c FDSQRT H S T X Y XR1 6420

16
c DHRND

c Z c ERRINT

c AN BN CI CONS DN ERF

c ERFC F FDEXP FN FNM1 FNM2

c FOUR GN GNM1 GNM2 ONE P

c PREV RNBC SCF SUM TN TOLER

c TRRTPI TWO ULCF ULCPS WN X

c Y YSQ

c FDCOS

c DCOS DSNCOS DXEXP FDCOS X

c FDEXP

c DEXP DSNCOS DXEXP FDEXP X

c FDLOG

c DLOG FDLOG X

c FDPCON

c D X

c FDSIN

c DSIN DSNCOS DXEXP FDSIN X

c FDSQRT

c DSQRT FDSQRT X

c FOURIA

c A AA AB AC AD BA

c BB FDCOS FDSIN R Y

c FNEC

c Y Z c GQUAD

c B BMA BPA C DELGQ STORE1

c STORE

2

c LBCONS

c DSNCOS DXEXP TRRTPI

c MMULT

c SUM X

c MRAISE

c SUM X

c MTRIAN

c SUM X

c MXTX

c AP

c MXTXP

c SUM XP

c OANOVA

c YSUM

c OCOEFF

c YSUM

c ORTHO

c DK2 FDSQRT SUM YSUM

c ORTHRV

c XP

c PROB

c A B C FDCOS FDEXP FDLOG

c FDSIN G ONE TA TB W

c X

c RFORMT

c Z c SORTSM

c SAVE SUM X

c STRUVE

c A B C DBEJ P Q

c R S X Z

XRl 6430

XR1 6440

XRl 6450

XRl 6460

XRl 6470

XRl 6480

XRl 6490

XRl 6500

XRl 6510

XRl 6520

XRl 6530

XRl 6550

XRl 6560

XRl 6570

XRl 6580

XRl 6590

XRl 6600

XRl 6610

XRl 6620

XRl 6630

XRl 6640

XRl 6650

XRl 6660

XRl 6670

XRl 6680

XRl 6690

XRl 6693

XRl 6695

XRl 6697

XRl 6700

XRl 6710

XRl 6720

XRl 6730

XRl 6740

XRl 6750

XRl 6760

XRl 6770

XRl 6780

XRl 6790

XRl 6800

XRl 6810

XRl 6820

XRl 6830

XRl 6840

XRl 6850

XRl 6860

XRl 6870

XRl 6873

XRl 6875

XRl 6880

XRl 6890

XRl 6900

XRl 6910

XRl 6920

XRl 6930

XRl 6940

XRl 6950

XRl 6960

XRl 6970

XRl 6980

17
c THERMO XR1 6990

c EXDIF EXX FDEXP FDLOG G QO XR1 7000

c Ql Q2 QQ x XR1 7010

c TRANSF XR1 7020

c SUM X XR1 7030

c TWOWAY XR1 7040

c DK2 FDSQRT SUM XR1 7050

c XR1 7060

c XR1 7070

c 3 ** ^VARIABLES EQUIVALEK EACH OTHER********************** XR1 7080

c XR1 7090

C - SUBPROGRAM VARIABLE EQUIVALENCE VARIABLE EQUIVALENCE XR1 7100

c XR1 7105

c ALLSUB XR1 7110

c SCRA A Lll LL1 XR1 7120

c L22 LL2 XR1 7130

c ARITH XR1 7140

c 11 11(1) 12 11(2) XR1 7150

c 13 11(3) 14 11(4) XR1 7160

c 15 11(5) XR1 7170

c ARYVEC XR1 7180

c X A XR1 7190

c BESSEL XR1 7200

c A(l) AA A(2001) B XR1 7210

c A(4001) W R RC XR1 7220

c IA IARGS KI KIND XR1 7230

c NR NRMAX XR1 7240

c CBEK XR1 7250

c SCRAT(1700) AA SCRAT(1800) AB XR1 7260

c COMPLX XR1 7270

c 11 IARGS (1) 12 IARGS (2) XR1 7280

c 13 IARGS (3) 14 IARGS (4) XR1 7290

c 15 IARGS (5) 16 IARGS (6) XR1 7300

c CORREL XR1 7310

c A(13301) AVG(l) A(13401) SD(1) XR1 7320

c A(13401) T(D XR1 7330

c DBEJ XR1 7340

c SCRAT(1201) S SCRAT(1451) T XR1 7350

c FLIP XR1 7360

c I IARGS (100) J IARGS (99) XR1 7370

c K IARGS (9 8) KK IARGS (9 7) XR1 7380

c M IARGS (96) MM IARGS (9 5) XR1 7390

c MMM IARGS (9 4) N IARGS (9 3) XR1 7400

c NN IARGS (9 2) A ARGS(l) XR1 7410

c FUNCT XR1 7420

c 11 11(1) 12 11(2) XR1 7430

c 13 H(3) 14 11(4) XR1 7440

c IFS XR1 7450

c 11 11(1) 12 11(2) XR1 7460

c 13 H(3) XR1 7470

c MMULT XR1 7480

c X A XR1 7490

c MRAISE XR1 7500

c X A XR1 7510

c MTRIAN XR1 7520

c X A XR1 7530

c MXTX XR1 7540

c A AP XR1 7550

c ONEWAY XR1 7560

c BLANK L(45) SLO L(22) XRI 7570

18
c HIGH L(18) A2C1) A(2701) XR1 7580

c A3 C.I) A(5401) A4(l) A(8101) XR1 7590

c A5(l) A(10801) Bid) ACD XR1 7600

c B2(l) A(541) B3(l) A(1081 XR1 7610

c B4(l) A(1621) B5(l) AC2161) XR1 7620

c B6(l) A(2701) B7fl) A(3241) XR1 7630

c B8(l) A(3781) B9(l) A(4321) XR1 7640

c B10(l) A(4861) XR1 7650

c OPONE XR1 7660

c IIRGS(l) KFMT(l) XR1 7670

c ORTHO XR1 7680

c IIRGS(l) KFMT(l) BCD IB XR1 7690

c ORTPLT XR1 7700

c IU A XR1 7710

c PLOT XR1 7720

c TIT ITLE(1,6) TITX ITLE(1,5) XR1 7730

c RC(1) X(l) PRINT A XR1 7740

c XO XMIN XI XMAX XR1 7750

c YO YMIN Yl YMAX XR1 7760

c LHEAD IH IPR A(200) XR1 7770

c PREPAK XR1 7780

c A IAA XR1 7790

c RFORMT XR1 7800

c C(D LCD BLANK L(45) XR1 7810

c PERIOD L(38) CPLUS L(40) XR1 7820

c CMINUS L(39) CASTER L(41) XR1 7830

c RPRINT XR1 7840

c NWIDTH(l) A(1001) NDECS(l) A(1101) XR1 7850

c NBLANK(l) AC1201) IRGS(l) A(1301) XR1 7860

c NCOUNT(l) A(1401) XR1 7870

c NWMAX(l) A(1601) NSIGDS(l) A(1701) XR1 7880

c AL(1) LCD NF(1) A(1801) XR1 7890

c NWM(l) AC1901) XR1 7900

c SORTSM XR1 7910

c X A XR1 7920

c STATIS XR1 7930

c A(101) ISA A(3226) SA XR1 7940

c THERMO XR1 7950

c A(l) QQ1 XR1 7960

c TRANSF XR1 79 70 c X A XR1 7980

c TWOWAY XR1 7990

c KFMT IIRGS ND1 KIND(IOO) XR1 8000

c ND2 KIND(99) ND3 KIND(98) XR1 8010

c ND4 KIND (9 7) ND5 KIND(96) XR1 8020

c ND6 KIND(95) ND7 KIND(94) XR1 80 30 c ND8 KIND (9 3) ND9 KIND(92) XR1 8040

c ND10 KIND(91) ND11 KIND(90) XR1 8050

c ND12 KIND(89) ND13 KIND(88) XR1 8060

c ND14 KIND(87) ND16 KIND(86) XR1 8070

c ND17 KIND(85) ND18 KIND(84) XR1 8080

c ND19 KIND (8 3) XR1 8090

c XSTOP XR1 8100

c END

ITEMP(l) A(l) XR1 8110

XR1 8120

19
Listing Of XREF2

c XREF2 XR2 10 c VERSION 5. 00 XREF2 5/15/70 XR2 20 c XR2 30 c l.**COMMAND NAMES , CORRESPONDING LI ,L2 VALUES AND SUBPROGRAM THAT XR2 40

c EXECUTES» THE INSTRUCTION************************************** XR2 50

c XR2 60

p*** FOLLOWING COMMANDS ARE EXECUTED IN OMNIT AND LI ,L2 ARE NOT DEFINED XR2 70 c XR2 80

c FINISH FORMAT A FORMAT B FORMAT C XR2 90 c FORMAT D FORMAT E FORMAT F HEAD XR2 100

c NOTE NOTE1 NOTE2 OMNITA XR2 110

c STOP TITLE 1 TITLE 2 TITLE 3 XR2 120 c TITLE 4 TITLE X TITLE Y XR2 130 c XR2 140

c**** FORMAT CALLS XFORMT XR2 150

C XR2 160

C XR2 170

p*AA COMMAND LI L2 CALLING COMMAND LI L2 CALLING XR2 180

c NAMES SUB NAMES SUB XR2 190

c XR2 200

c AADD 18 4 MATRIX AAVERA 18 10 COALES XR2 210

c ABRIDG 6 1 ABRIDG ABRIDG A 6 2 ABRIDG XR2 220

c ABRIDG B 6 3 ABRIDG ABRIDG C 6 4 ABRIDG XR2 230 c ABRIDG D 6 5 ABRIDG ABRIDG E 6 6 ABRIDG XR2 240

c ABRIDG F 6 7 ABRIDG ABS 12 31 FUNCT XR2 250 c ABSOLU 12 31 FUNCT ACCURA 11 6 ARITH XR2 260

c ACOALES 18 9 COALES ACOS 12 6 FUNCT XR2 270

c ACOSD 12 14 FUNCT ACOSH 12 28 FUNCT XR2 280 c ACOT 12 8 FUNCT ACOTD 12 16 FUNCT XR2 290

c ACOTH 12 30 FUNCT ADD 11 1 ARITH XR2 300

c ADEFIN 15 1 MOP ADIV 18 7 MATRIX XR2 310

c ADIVID 18 7 MATRIX AERASE 15 2 MOP XR2 320 c AMOVE 23 6 MOVE AMULT 18 6 MATRIX XR2 330

c AMULTI 18 6 MATRIX ANTILO 12 22 FUNCT XR2 340 c APRINT 4 1 APRINT APRINT A 4 2 APRINT XR2 350

c APRINT B 4 3 APRINT APRINT C 4 4 APRINT XR2 360

c APRINT D 4 5 APRINT APRINT E 4 6 APRINT XR2 370

c APRINT F 4 7 APRINT APROPE 27 2 MPROP XR2 380

c ARAISE 18 8 MATRIX ASIN 12 5 FUNCT XR2 390 c ASIND 12 13 FUNCT ASINH 12 27 FUNCT XR2 400 c ASUB 18 5 MATRIX ASUBTR 18 5 MATRIX XR2 410

c ATAN 12 7 FUNCT ATAND 12 15 FUNCT XR2 420 c ATANII 12 29 FUNCT ATOMIC 31 3 THERMO XR2 430 c ATRANS 18 3 MATRIX AVERAG 20 4 MSCROW XR2 440 c AZERO 15 2 MOP BACKSPACE TAPE50 1 TAPOP2 XR2 450 c BEGIN 14 1 BEGIN BESIN 30 38 BESSEL XR2 460

c BESJN 30 32 BESSEL BESKN 30 39 BESSEL XR2 470 c BIONE 30 6 BESSEL B I ZERO 30 5 BESSEL XR2 480

c BJONE 30 2 BESSEL BJZERO 30 1 BESSEL XR2 490 c BKONE 30 8 BESSEL BKZERO 30 7 BESSEL XR2 500

c BOLDIS 31 9 THERMO BYONE 30 4 BESSEL XR2 510

c BYZERO 30 3 BESSEL CADD 32 1 COMPLX XR2 520

20
C*** COMMAND LI L2 CALLING COMMAND LI L2 CALLING

C

C

C

NAMES SUB NAMES SUB

CDIVID 32 4 COMPLX CEIONE 30 26 BESSEL XR2 530 C CEIZER 30 25 BESSEL CEKONE 30 28 BESSEL XR2 540

C CEKZER 30 27 BESSEL CENSOR 25 3 SELECT XR2 550

C CERF 21 19 FNEC CGS 13 10 PHYCON XR2 560

C CHANGE 21 13 CHANGE CIONE 30 22 BESSEL XR2 570

C CI ZERO 30 21 BESSEL CKONE 30 24 BESSEL XR2 580 C CKZERO 30 23 BESSEL CLOSE UP 23 1 MISC2 XR2 590

C CMULTI 32 3 COMPLX COMPAR 14 15 IFS XR2 600

C CORREL 24 11 CORREL COS 12 2 FUNCT XR2 610

C COSD 12 10 FUNCT COSH 12 24 FUNCT XR2 620

c COT 12 4 FUNCT COTD 12 12 FUNCT XR2 630

c COTH 12 26 FUNCT COUNT 23 2 MISC2 XR2 640

c CPOLAR 32 6 COMPLX CREAD TAPE 46 1 TAPOP2 XR2 650

c CREAD TAPE A 46 2 TAPOP2 CREAD TAPE B 46 3 TAPOP2 XR2 660

c CREAD TAPE C 46 4 TAPOP2 CREAD TAPE D 46 5 TAPOP2 XR2 670

c CREAD TAPE E 46 6 TAPOP2 CREAD TAPE F 46 7 TAPOP2 XR2 680

c CRECTA 32 5 COMPLX CSET TAPE 49 1 TAPOP2 XR2 690

c CSUBTR 32 2 COMPLX CTOF 31 1 THERMO XR2 700

c DEFINE 21 3 DEFINE DEMOTE 23 11 PDMOTE XR2 710

c DIM 23 12 DIMENS DIMENS 23 12 DIMENS XR2 720 c DIV 11 4 ARITH DIVIDE 11 4 ARITH XR2 730 c DUMMY A 54 2 DUMMYA DUMMY B 54 3 DUMMYB XR2 740 c DUMMY C 54 4 DUMMYC DUMMY D 54 5 DUMMYD XR2 750 c DUMMY E 54 6 DUMMYE DUMMY F 54 7 DUMMYF XR2 760

c DUPLIC 23 5 MISC2 EINSTE 31 5 THERMO XR2 770 c ELLIPT FIRST 30 30 BESSEL ELLIPT SECOND 30 31 BESSEL XR2 780

c ENDFIL TAPE 50 1 TAPOP2 ERASE 21 10 ERASE XR2 790 c ERROR 21 18 FNEC EXGHAN 21 11 EXCHNG XR2 800

c EXECUT 14 3 REPINC EXIONE 30 10 BESSEL XR2 810 c EXIZER 30 9 BESSEL EXKONE 30 12 BESSEL XR2 820

c EXKZER 30 11 BESSEL EXP 12 18 FUNCT XR2 830

c EXPAND 23 4 MISC2 EXPONE 12 18 FUNCT XR2 840 c EXTREM 29 4 CMSEPA F PROBAB 24 5 FPROB XR2 850 c FIT 22 3 ORTHO FIXED 13 3 FIXFLO XR2 860 c FLEXIB 13 12 FIXFLO FLIP 21 12 FLIP XR2 870 c FLOATI 13 4 FIXFLO FRACTI 12 33 FUNCT XR2 880 c FREQUE 24 10 FRDIST FTOC 31 2 THERMO XR2 890 c GAUSS QUADRA 24 4 GQUAD GENERA 13 1 GENER XR2 900

c HARMON 30 37 BESSEL HERMIT 19 3 ALLSUB XR2 910

c HIERAR 21 14 SORDER HISTOG 24 8 HISTGM XR2 920

c IFEQ 14 10 IFS IFGE 14 12 IFS XR2 930

c IFGT 14 11 IFS IFLE 14 14 IFS XR2 940

c IFLT 14 9 IFS IFNE 14 13 IFS XR2 950

c INCREM 14 6 REPINC INSERT 29 3 CMSEPA XR2 960

c INTEGE 12 32 FUNCT INTERP 25 4 INTERP XR2 970

c INTJO 30 29 BESSEL INVERT 16 1 INVERT XR2 980

c ISETUP 28 2 ITERAT ISOLAT 28 3 ITERAT XR2 990

c ITERAT 28 1 ITERAT KBIONE 30 14 BESSEL XR2 1000

c KBIZER 30 13 BESSEL KBKONE 30 16 BESSEL XR2 1010

c KBKZER 30 15 BESSEL KEXION 30 18 BESSEL XR2 1020

c KEXIZE 30 17 BESSEL KEXKON 30 20 BESSEL XR2 1030

c KEXKZE 30 19 BESSEL LAGUER 19 2 ALLSUB XR2 1040

c LEGEND 19 5 ALLSUB LIST 21 15 LIST XR2 1050

c LOG 12 20 FUNCT LOGE 12 20 FUNCT XR2 1060

c LOGTEN 12 21 FUNCT MADD 18 1 MATRIX XR2 1070

c MATCH 25 5 SELECT MAX 21 5 EXTREM XR2 1080

c MAXIMU 21 5 EXTREM MAXMIN 29 4 CMSEPA XR2 1090

21
r*A* COMMAND LI L2 CALLING COMMANI) LI L2 CALLING

c

cc

NAMES SUB NAMES SUB

MDEFIN 15 1 MOP MDIAGO 15 4 MOP XR2 1100

c ME I GEN 17 5 MEIGEN MERASE 15 2 MOP XR2 1110

c MI DENT 15 3 MOP MIN 21 6 EXTREM XR2 1120

c MINIMU 21 6 EXTREM MINVER 16 1 INVERT XR2 1130

c MKRON 17 3 MKRON MMATVE 26 3 EXPCON XR2 1140

c MMOVE 23 6 MOVE MMULT 17 1 MMULT XR2 1150

c MMULTI 17 1 MMULT MOLWT 31 4 THERMO XR2 1160

c MORTHO 22 5 ORTHO MOVE 23 6 MOVE XR2 1170

c MPRINT 7 1 APRINT MPRINT A 7 2 APRINT XR2 1180

c MPRINT B 7 3 APRINT MPRINT C 7 4 APRINT XR2 1190

c MPRINT D 7 5 APRINT MPRINT E 7 6 APRINT XR2 1200

c MPRINT F 7 7 APRINT MPROPE 27 1 MPROP XR2 1210

c MRAISE 17 2 MRAISE MSCALA 18 6 MATRIX XR2 1220

c MSUB 18 2 MATRIX MSUBTR 18 2 MATRIX XR2 1230

c MTRANS 18 3 MATRIX MTRIAN 17 4 MTRIAN XR2 1240

c MULT 11 3 ARITH MULTIP 11 3 ARITH XR2 1250

c MVECDI 26 1 EXPCON MVECMA 26 2 EXPCON XR2 1260

c MZERO 15 2 MOP M(AD) 52 1 MDAMAD XR2 1270

c M(AV) 53 1 ARYVEC M(DA) 52 2 MDAMAD XR2 1280

c M(V'A) 53 2 ARYVEC M(XAX') 51 3 MXTX XR2 1290

c MQC'AX] 51 "j MXTX M(XX') 51 1 MXTX XR2 1300

c M(X'X) 51 2 MXTX NEGEXP 12 19 FUNCT XR2 1310

c NEW PAGE 13 8 PAGE NHISTO 24 9 HISTGM XR2 1320

c NO LIST 21 16 LIST NORMLA 19 1 ALLSUB XR2 1330

c NPRINT 8 1 PRINTX NPRINT A 8 2 PRINTX XR2 1340

c NPRINT B 8 3 PRINTX NPRINT C 8 4 PRINTX XR2 1350

c NPRINT D 8 5 PRINTX NPRINT E 8 6 PRINTX XR2 1360

c NPRINT F 8 7 PRINTX NULL 21 17 XECUTE XR2 1370

c ONEWAY 24 13 ONEWAY ORDER 21 9 SORDER XR2 1380

c PAGE PLOT 13 6 PLOT PARPRO 20 2 MSCROW XR2 1390

c PARSUM 20 1 MSCROW PARTFU 31 8 THERMO XR2 1400

c PERFOR 14 3 REPINC PFATOM 31 7 THERMO XR2 1410

c PFTRAN 31 6 THERMO PLOT 13 5 PLOT XR2 1420

c POLYFI 22 1 ORTHO PRINT 2 1 PRINTX XR2 1430

c PRINT A 2 2 PRINTX PRINT B 2 3 PRINTX XR2 1440

c PRINT C 2 4 PRINTX PRINT D 2 5 PRINTX XR2 1450

c PRINT E 2 6 PRINTX PRINT F 2 7 PRINTX XR2 1460

c PRINT NOTE 13 13 NOTEPR PRODUC 21 2 PROROW XR2 1470

c PROMOT 23 10 PDMOTE PUNCH 3 1 PRINTX XR2 1480

c PUNCH A 3 2 PRINTX PUNCH B 3 3 PRINTX XR2 1490

c PUNCH C 3 4 PRINTX PUNCH D 3 5 PRINTX XR2 1500

c PUNCH E 3 6 PRINTX PUNCH F 3 7 PRINTX XR2 1510

c RAISE 11 5 ARITH RANKS 24 3 RANKS XR2 1520

c READ 5 1 READX READ A 5 2 READX XR2 1530

c READ B 5 3 READX READ C 5 4 READX XR2 1540

c READ D 5 5 READX READ E 5 6 READX XR2 1550

c READ X 5 7 READX READ TAPE 45 1 TAPOP2 XR2 1560 c READ TAPE A 45 2 TAPOP2 READ TAPE B 45 3 TAPOP2 XR2 1570

c READ TAPE C 45 4 TAPOP2 READ TAPE D 45 5 TAPOP2 XR2 1580

c READ TAPE E 45 6 TAPOP2 READ TAPE F 45 7 TAPOP2 XR2 1590

c REPEAT 14 3 REPINC RESET 1 1 RESET XR2 1600

c RESET V 1 3 RESET RESET W 1 4 RESET XR2 1610

c RESET X 1 5 RESET RESET Y 1 6 RESET XR2 1620

c RESET z 1 7 RESET RESTOR 14 8 REPINC XR2 1630

c REWIND TAPE 50 1 TAPOP2 RMS 20 3 MSCROW XR2 1640

c ROUND 13 14 FNEIC ROW SUM 21 1 PROROW XR2 1650

c ROWSUM 21 1 PROROW SAPROP 27 4 MPROP XR2 1660

22
p&rt* COMMAND LI L2 CALLING COMMAND LI L2 CALLING

cc

c

NAMES SUB NAMES SUB

SCAN 14 2 BEGIN SCORRE 24 12 CORREL XR2 1670

c SEARCH 25 2 SELECT SELECT 25 1 SELECT XR2 1680

c SEPARA 29 2 CMSEPA SET 13 2 SET XR2 1690

c SET TAPE 48 1 TAPOP2 SFIT 22 4 ORTHO XR2 1700

c SHORTE 23 3 MISC2 SI 13 11 PHYCON XR2 1710

c SIN 12 1 FUNCT SIND 12 9 FUNCT XR2 1720

c SINH 12 23 FUNCT SKIP TAPE 50 1 TAPOP2 XR2 1730

c SMPROP 27 3 MPROP SOLVE 16 2 INVERT XR2 1740

c SONEWA 24 14 ONEWAY SORT 21 8 SORDER XR2 1750

c SPACE 13 9 SPACE SPOLYF 22 2 ORTHO XR2 1760

c SQRT 12 17 FUNCT SQUARE 12 34 FUNCT XR2 1770

c SSTATI 24 2 STATIS STATIS 24 1 STATIS XR2 1780

c STRUVE ONE 30 36 BESSEL STRUVE ZERO 30 35 BESSEL XR2 1790

c STWOWA 24 7 TWOWAY SUB 11 2 ARITH XR2 1800

c SUBTRA 11 2 ARITH SUM 20 5 MSCROW XR2 1810

c TAN 12 3 FUNCT TAND 12 11 FUNCT XR2 1820

c TANH 12 25 FUNCT TCHEBY 19 6 ALLSUB XR2 1830

c TWOWAY 24 6 TWOWAY UCHEBY 19 4 ALLSUB XR2 1840

c UNIFOR RANDOM 24 15 FNKC WRITE TAPE 47 1 TAPOP2 XR2 1850

c MITE TAPE A 47 2 TAPOP2 WRITE TAPE B 47 3 TAPOP2 XR2 1860

c WRITE TAPE C 47 4 TAPOP2 WRITE TAPE D 47 5 TAPOP2 XR2 1870

c WRITE TAPE E 47 6 TAPOP2 WRITE TAPE F 47 7 TAPOP2 XR2 1880

c ZEROS BJONE 30 34 BESSEL ZEROS BJZERO 30 33 BESSEL XR2 1890

c END

XR2

XR2

1900

1910

23
Listing Of XREF3

C XREF3

C VERSION 5.00 XREF3 5/15/70

C

C THE COMMENTS CONTAINED IN THIS SECTION ARE FOR SYSTEM

C IMPLEMENTATION.

C l.**CHANGES NEEDED IF WORK SHEET AND SCRATCH AREAS ARE TO BE

C MODIFIED****************************************************** XR3

C

C OMNITAB HAS A WORK SHEET OF 12500 MACHINE WORDS. IF THIS IS TO

C BE INCREASED OR DECREASED, THE FOLLOWING CORRECTIONS WILL BE

C NECESSARY.

C

C LET NSIZRC BE THE SIZE OF WORK SHEET REQUIRED. (NSIZRC

C MUST BE A CONSTANT AND NOT A VARIABLE .

)

C

C THEN LET

C

C NSIZR2 = INTEGRAL PART OF (NSIZRC + 1000) /2

CC NSIZR4 = INTEGRAL PART OF NSIZRC/4

C

C NSIZR5 = INTEGRAL PART OF (NSIZRC + 1000) /5

C

C NSIZR6 = INTEGRAL PART OF NSIZRC/5

C

C

C I THE STATEMENTS

C

C COMMON/BLOCRC/NRC,RC(12600)

C EQUIVALENCE (ARGS(l) ,RC( 12 501))

C

C MUST BE CHANGED TO

CC CQMMON/BLOCRC/NRC,RC(NSIZRC+100)

C EQUIVALENCE (ARGS(l) ,RC(NSIZRC+1))

C

C IN THE FOLLOWING SUBPROGRAMS:

C

C ABRIDG ADRESS ALLSUB APRINT ARITH ARYVEC BEGIN BESSEL CHANGE

C QIKCOL CKIND CMSEPA COALES COMPLX CORREL DEFINE DIMENS ERASE

C ERROR EXCHNG EXPAND EXPCON EXTREM FIXFLO FLIP FPROB FNEC

C FNEIC FNKC FRDIST FUNCT GENER GQUAD HISTGM IFS INTERP

C INVERT ITERAT LIST MATRIX MDAMAD MEIGEN MISC2 MKRON MMULT

C MOP MOVE MPROP MRAISE MSCROW MTRIAN MXTX OMNIT ONEWAY

C OPONE ORTHO ORTPLT OUTPUT PDMOTE PLOT PRINTX PROROW PUNCH

C RANKS READQ READX REPINC RESET RPRINT SELECT SET SETQ

C SETUP SORDER SPACE STATIS STORE THERMO TRANSF TWOWAY VECTOR

C XCMNIT XPND XSTOP

C

C II THE STATEMENT

CC COMMON/SCRAT/NS,NS2,A(13500)

24

XR3 10 XR3 20 XR3 30 XR3 40

XR3 50 XR3 60 XR3 70 XR3 80 XR3 90

XR3 100

XR3 110

XR3 120

XR3 130

XR3 135

XR3 140

XR3 141

XR3 142

XR3 143

XR3 144

XR3 145

XR3 146

XR3 147

XR3 148

XR3 149

XR3 150

XR3 160

XR3 170

XR3 180

XR3 190

XR3 200

XR3 210

XR3 220

XR3 230

XR3 240

XR3 250

XR3 260

XR3 270

XR3 280

XR3 290

XR3 300

XR3 310

XR3 320

XR3 330

XR3 340

XR3 350

XR3 360

XR3 370

XR3 380

XR3 390

XR3 400

XR3 410

XR3 420
C XR3 430

C MUST BE CHANGED TO XR3 440

C XR3 450 C COMMON/SCRAT/NS,NS2,A(NSIZRC+1000) XR3 460

C XR3 470

C IN THE FOLLOWING SUBPROGRAMS: XR3 480

C XR3 490

C ALLSUB APRINT ARYVEC BESSEL CBEK CMSEPA COALES CCMPLX CORREL XR3 500

C DBEJ EXPCON FNEC FNEIC FNKC FUNCT HISTGM INTERP INVERT XR3 510

C ITERAT MATRIX MDAMAD ME I GEN MISC2 MKRON MMULT MOP MPROP XR3 520

C MRAISE MTRIAN MXTX OANOVA OCOEFF OCOVAR ONEWAY OPONE ORTHO XR3 530

C ORTPLT PLOT PREPAK PROROW RANKS RPRINT SELECT SETUP SORDER XR3 540

C SORTSM STATIS TAPOP2 THERMO TRANSF TWOWAY XSTOP XR3 550

C XR3 560

C XR3 570 C THE STATEMENT XR3 580

C XR3 590 C COMMON/SCRAT/ NS,NS2,SCRAT( 13500) XR3 600

C XR3 610

C MUST BE CHANGED TO XR3 620

C XR3 630

C CCMMON/SCRAT/NS,NS2,SCRAT(NSIZRC+1000) XR3 640

C XR3 650

C IN THE FOLLOWING SUBPROGRAMS: XR3 660

C XR3 670

C CBEK DBEJ XR3 680

C XR3 690

C ( THE VARIABLE A WAS CHANGED TO SCRAT BECAUSE A WAS USED XR3 700 C IN THE SUBPROGRAMS) XR3 710 C XR3 720

C III THE STATEMENTS XR3 730 C XR3 735 C NS=13500 XR3 740

C NRC=12500 XR3 750 C XR3 760

C MUST BE CHANGED TO XR3 770 C XR3 780

C NS=NSIZRC+1000 XR3 790 C NRC=NSIZRC XR3 800

C XR3 810

C IN THE SUBPROGRAM SETUP. XR3 820

C XR3 830 C IV THE STATEMENTS XR3 840

C XR3 850 C DIMENSION SAC3125, 3), ISA(3125) XR3 860 C EQUIVALENCE (A(101), ISA), (A(3226),SA),(NRMAX,NARMAX) XR3 870 C XR3 880 C MUST BE CHANGED TO XR3 890 C XR3 900

C DIMENSION SA(NSIZ4, 3), ISA(NSIZR4) XR3 910

C EQUIVALENCE (A(101) ,ISA) , (A(NSIZR4+101) ,SA) , (NRMAX,NARMAX) XR3 920

C XR3 930

C IN THE SUBPROGRAM STATIS. XR3 940

C XR3 950

C V THE STATEMENTS XR3 960

C XR3 970

C NROW=201 XR3 980

C NCOL=62 XR3 990

C XR3 1000

C MUST BE CHANGED TO XR3 1010

25
C XR3 1020

C NROW=KROW XR3 1030

C NCOL=KCOL XR3 1040

C WHERE KCOL*KROW IS LESS THAN OR EQUAL TO NSIZRC XR3 1050

C XR3 1060

C IN THE SUBPROGRAM XOMNIT XR3 1070

C XR3 1075

C VI THE STATEMENT XR3 1080

C XR3 1090

C EQUIVALENCE (A(13301) ,AVG(1)) , (A(13401) ,SD(1) ,T (1)) XR3 1100

C XR3 1110

C MUST BE CHANGED TO XR3 1120

C XR3 1130

C EQUIVALENCE (A(NSIZRC+801) ,AVG(1)) , (A(NSIZRC+901) ,SD(1) ,T(1)) XR3 1140

C XR3 1150

C IN THE SUBPROGRAM CORREL XR3 1160

C XR3 1165

CVII THE STATEMENT XR3 1170

C XR3 1180

C DIMENSION QQ(6750) XR3 1190

C XR3 1200

C MUST BE CHANGED TO XR3 1210

C XR3 1220

C DIMENSION QQ(NSIZR2) XR3 1230

C XR3 1240

C IN THE SUBPROGRAM THERMO XR3 1250

C XR3 1260

CVIII THE STATEMENTS XR3 1270

C XR3 1280

C DIMENSION A2 (2700),A3(2700) ,A4(2700) ,A5 (2700) XR3 1290

C EQUIVALENCE (A2(l) ,A(2701)

)

XR3 1300

C EQUIVALENCE (A3(l) ,A(5401)) , (A4(l) ,A(8101) ) , (A5(l) ,A(10801)

)

XR3 1310

C DIMENSION B1(540),B2(540),B3(540),B4(540),B5(540),B6(540), XR3 1320

C 1B7(540),B8(540),B9(540),B10(540) XR3 1330

C EQUIVQLENCE (B1(1),A(1)),(B2(1) ,A(541)),(B3(1),A(1081)), XR3 1340

C 1(B4(1),A(1621)),(B5(1),A(2161)),(B6(1),A(2701)),(B7(1),A(3241)), XR3 1350

C 2(B8(1),A(3781)),(B9(1),A(4321)),(B10(1),A(4861)) XR3 1360

C XR3 1370

C MUST BE CHANGED TO XR3 1380

C XR3 1395

C DIMENSION A2(NSIZR5),A3(NSIZR5),A4(NSIZR5),A5(NSIZR5) XR3 1400

C EQUIVALENCE (A2(l) ,A(NSIZR5+1)) XR3 1410

C EQUIVALENCE (A3(l) ,A(2*NSIZR5+1)) ,(A4(1) ,A(3*NSIZR5+1)) , (A5(l) , XR3 1420

C 1A(4*NSIZR5+1)) XR3 1430

C DIMENSION B1(NSIZR6),B2(NSIZR6),B3(NSIZR6),B4(NSIZR6),B5(NSIZR6), XR3 1440

C 1B6(NSIZR6),B7(NSIZR6),B8(NSIZR6),B9(NSIZR6),B10(NSIZR6) XR3 1450

C EQUIVALENCE (Bl(l) ,A(1) ) , (B2(l) ,A(NSIZR6+1) ) , (B3(l) ,A(2*NSIZR6+1))XR3 1460

C 1,(B4(1),A(3*NSIZR6+1)),(B5(1),A(4*NSIZR6+1)),(B6(1),A(5*NSIZR6+1))XR3 1470

C 2,(B7(1),A(6*NSIZR6+1)),(B8(1),A(7*NSIZR6+1)),(B9(1),A(8*NSIZR6+1))XR3 1480

C 3,(B10(1),A(9*NSIZR6+1)) XR3 1490

C XR3 1500

C IN THE SUBPROGRAM ONEWAY XR3 1510

C XR3 1520

C 2.**CHANGING THE NUMBER OF HEADINGS PERMITTED AND PACKING OF XR3 1530

FORMATS  ** *************************************** ********** xr3 1540

C XR3 1550

2. A INCREASING OR DECREASING THE NUMBER OF HEADINGS FOR COLUMNS. XR3 1560

C XR3 1570

2.A.1 THE VERSION DISTRIBUTED ALLOWS HEADINGS FOR 50 COLUMNS. ANY XR3 1580

50 COLUMNS MAY BE HEADED. IF MORE THAN 50 COLUMNS ARE XR3 1590

26
C HEADED, THE HEADINGS OF THE COLUMNS HEADED EARLIEST WILL BE XR3 1600

C DELETED. (I.E. HEADS ARE STORED IN A PUSH DOWN TABLE WITH XR3 1610

C END OF TABLE PUSHED OUT. THIS APPLIES NO MATTER HOW MANY XR3 1620

C HEADINGS AN INSTALLATION PERMITS.) XR3 1630

C XR3 1640

C 2. A. 2 LET NHEAD= NUMBER OF COLUMN HEADINGS PERMITTED. THEN THE XR3 1650

C FOLLOWING CHANGES HAVE TO BE MADE IN THE SUBPROGRAM PREPAK. XR3 1660

C XR3 1665

C I DIMENSION IFMT(12, 6), IHEAD(5, 50) XR3 1670

C XR3 1675

C MUST BE CHANGED TO XR3 1680

C XR3 1685

C DIMENSION IFMT(12, 6), IHEAD(5,NHEAD) XR3 1690

C XR3 1695

C II DATA I 1/12/, LA/50/ XR3 1700

C XR3 1705

C MUST BE CHANGED TO XR3 1710

C XR3 1715

C DATA II/12/,LA/NHEAD/ XR3 1720

C XR3 1730

C (NOTE: NHEAD STANDS FOR A INTEGER CONSTANT) XR3 1740

C XR3 1750

C 2.B FORMATS MODIFICATION FOR COMPUTER SYSTEMS WHICH DO NOT PACK XR3 1760

C SIX CHARACTERS PER MACHINE WORD XR3 1770

C XR3 1780

C 2.B.1 FORMATS MUST BE PACKED TO THE MAXIMUM NUMBER OF CHARACTERS XR3 1790

C PER WORD IN ORDER FOR THE NH CONVERSION TO WORK XR3 1800

C XR3 1810

C 2.B.2 LET NCHAR=NUMBER OF CHARACTERS THAT CAN BE PACKED PER WORD XR3 1820

C FOR A PARTICULAR CONFIGURATION. THEN THE FOLLOWING CHANGES XR3 1830

C HAVE TO BE MADE IN THE SUBPROGRAM PREPAK, WHERE K IS XR3 1840

C DEFINED AS. XR3 1850

C XR3 1860

C K=72/NCHAR+M M=0 IF 72/NCHAR HAS NO REMAINDER XR3 1870

C M=l IF 72/NCHAR HAS A REMAINDER XR3 1880

C XR3 1885

C I DIMENSION IFMT(12,6),IHEAD(5,50) XR3 1890

C XR3 1895

C MUST BE CHANGED TO XR3 1900

C XR3 1905

C DIMENSION IFMT(K, 6) ,IHEAD(5,50) XR3 1910

C XR3 1915

C II DATA 11/12/, LA/50/ XR3 1920

C XR3 1925

C MUST BE CHANGED TO XR3 1930

C XR3 1935

C DATA II/K/, LA/50/ XR3 1940

C XR3 1945

C III 230 FORMAT (12A6) XR3 1950

C XR3 1955

C MUST BE CHANGED TO XR3 1960

C XR3 1965

C 230 FORMAT (KA=NCHAR) XR3 1970

C XR3 1980

C (NOTE: NHEAD AND K STAND FOR AN INTEGER CONSTANT) XR3 1990

C XR3 2000

C 2.B.3 BECAUSE A FORMAT COMMAND IS TREATED THE WAY ALL OTHER XR3 2010

C OMNITAB COMMANDS ARE TREATED, THE FORMAT COMMAND IS XR3 2020

C DUMPED ON THE SCRATCH TAPE (OR OTHER I/O) AND REREAD XR3 2030

C IN THE PACK MODE USING THE FORMAT 230. THIS IS XR3 2040

27
C NECESSARY IN ORDER TO MAINTAIN MACHINE INDEPENDENCE. XR3 2050

C XR3 2060

C **OVERLAY OR SEGMENTATION*************************************** XR3 2070

C XR3 2080

C 3. A SEGMENTATION AS USED ON THE NBS UNIVAC 1108. XR3 2090

C XR3 2100

C I PART 1 IS RESIDENT IN MEMORY AT ALL TIMES, XR3 2110

C CONTAINS MOST OFTEN USED ROUTINES, SUBPROGRAMS XR3 2120

C NEEDED BY MORE THAN ONE OF THE OTHER PARTS, AND XR3 2130

C ALL OF LABELED COMMON (EXCLUDING LABELED COMMON XR3 2140

C FOR BESSEL SUBPROGRAMS.) XR3 2150

C XR3 2160

C II ONLY ONE OF THE OTHER PARTS (PART 2 THRU PART 14) XR3 2170

C IS IN MEMORY AT A PARTICULAR MOMENT. XR3 2180

C XR3 2190

C 3.B OUTLINE OF SEGMENTATION XR3 2200

C XR3 2210

C PART 1 XR3 2220

C XR3 2230

c ************************************** XR3 2240

C * AARGS ACCDIG ADRESS AERR ARITH * XR3 2250

C * ASTER BEGIN CHKCOL CKIND DIMENS * XR3 2260

C * ERROR EXPAND FCOS FDCOS FDEXP * XR3 2270

C * FDLOG FDPCON FDSIN FDSQRT FEXP * XR3 2280

C * FEXP2 * XR3 2290

C * FLOG FLOG10 FSIN FSQRT FTANH ** XR3 2300

C * FUNCT GENER HEADS IFS INFERR * XR3 2310

C * INPUT LIST LOCATE LOOKUP * XR3 2320

C * MTXCHK MXTXP NNAME NONBLA * XR3 2330

C * OMCONV OMNIT OMNSYM OUTPUT PACK * XR3 2340

C * PAGE PHYCON PREPAK PROB READQ * XR3 2350

C * READX REPINC RESET RFORMT RNDOWN * XR3 2360

C * RPRINT SET SORTSM SPACE SYMV * XR3 2370

C * VARCON VECTOR XECUTE XPND XSTOP * XR3 2380

C * * XR3 2390

C * PLUS LIBRARY FUNCTIONS (SIN, ETC) * XR3 2400

C * * XR3 2410

C * * XR3 2420

C * ALL OF LABELED COMMON * XR3 2430

C * * XR3 2440

C * ABCDEF BLOCKA BLOCKB BLOCKC BLOCKD * XR3 2450

C * BLOCKE BLOCKX BLOCRC CODE CODETP * XR3 2460

C * CONLB2 CONSLB CONSTS DCONL2 DCONLB * XR3 2470

C * FMAT HEADER ICODE ICODTP KFMT * XR3 2480

C * PCONST PKSWT QRS SCRAT SPRV * XR3 2490

C * TAPE * XR3 2500

C * * XR3 2510

C * * XR3 2520

C * THERE ARE 4 BLOCK DATA PROCEDURES. * XR3 2530

* THE COMMENT OR 2ND CARD HAS THE * XR3 2540

C * LABEL NAME. THE 4 PROCEDURES ARE * XR3 2550

C * * XR3 2560

C * BLOCK LBCONS LOOKTB PHYSIC * XR3 2570

C * * XR3 2580

Q ************************************** y^7 2590

C * XR3 2600

C * XR3 2610

C * XR3 2620

C * XR3 2630

C * XR3 2640

28
Q A********* A ************************************************** XR3

C 1*1 1 1 1 1 1 * 1 * XR3

CPART * 2 * 1P4RT * 4 XR3

C * P **********

C * ARYVEC *

* EXPCON *

* MATRIX *

* MDAMAD *

* MKRON *

* MMULT

* MOP

* MRAISE

* MTRIAN

* MXTX

* STORMT

* TRANSF

**********

1

11

1

1

1

1

1

11

1

1

1 *

CPART 1 5 p **********

C * ABRIDG * * APRINT *

* FIXFLO *

* NOTEPR *

* PRINTX * * PUNCH *

* SETQ

* SETUP

STMT

STORE

TAPOP

TAPOP2

* XFORMT

C * XHEAD

C * XOMNIT

p **********

C

cc

cCPART * 8 C * **********

* COALES *

* pppT *

* HDIAG *

* MEIGEN *

* ONEWAY * RANKO *

* TpPTPT * **********

1

11

1

11

11

11

I1

1 1*1

*

*

*

**

**

*** **********

*

C

C

C

C

CC

C

C

c

C

C

C

1PART * 3

1 *

i **********

1 * ALLSUB *

1 * CHANGE *

1 * CMSEPA * 1 CMPARA * EXCHNG * FLIP *

PART

1

1 9

1 **********

* TWOWAY * **********

INTERP * 1 INTRP * 1 ITERAT * 1 PROROW * 1 SELECT * 1

* SORDER * 1 ********** i

PART 1 6 **********

* COMPLX

* DEFINE

* ERASE

* EXTREM

* GQUAD

* MISC2

* MSCROW

* MOVE

* PDMOTE

********** i1

11

1

1

1

11

* 1 XR3

*

] *

* 1 **********Y]?'Z

*

1 * DETRNK *XR3

*

1 * INVCHK *XR3

* 1 * INVERT *XR3

* 1 * MPROP *XR3

* 1 * ORTHRV *XR3

*

1 * PROCHK *XR3

* 1 * PVTRI *XR3

* RCSUM *XR3

1 * SPINV *XR3

* 1 * SKSYMV *XR3

* i **********vp2

*

1 * !

* 1 *

!

* _. **

*

**

**

*

*

*

*

*

*

*

**

*

*

*

*

**

**********

*

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

PART 1 7 XR3

******** **vp3

* DUMMYA *XR3

* DUMMYB *XR3

* DUMMYC *XR3

* DUMMYD *XR3

* DUMMYE *XR3

* DUMMYF *XR3 ********* *XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

******** **XR3

* DHRND *XR3

* ERRINT *XR3

* FNEC *XR3

* FNEIC *XR3

* FNKC *XR3

* PLOT *XR3

* RNJBK *XR3

* THERMO *XR3

**********^P3

XR3

XR3

PART * 10

*

2650

2660

2670

2680

2690

2700

2710

2720

2730

2740

2750

2760

2770

2780

2790

2800

2810

2820

2830

2840

2850

2860

2870

2880

2890

2900

2910

2920

2930

2940

2950

2960

2970

2980

2990

3000

3010

30 20 3030

3040

3050

3060

3070

3080

3090

3100

3110

3120

3130

3140

3150

3160

3170

3180

3190

3200

3210

3220

3230

3240

29
c

C 1 — 11 1

11 XR3

XR3

3250

3260

CPART 1 11 1 PART 1 12 1 XR3 3270

C 1 1 1 1 XR3 3280

r ********** 1 ********** 1 XR3 3290

C * FPROB * 1 * BJORCK * 1 XR3 3300

C * FRDIST * 1 * CORREL * 1 XR3 3310

C * FREQCY * 1 * CSPINV * 1 XR3 3320

C * HISTGM * 1 * INVCOR * 1 XR3 3330

C * STATIS * 1

* MIST * 1 XR3 3340

c ********** 1 * RANKS *

1 XR3 3350

C 1 * RANKX * 1 XR3 3360

c 1

•kkkkkk-k'k'k'k 1 XR3 3370

c 1 1 XR3 3380

c 1 1 XR3 3390

c 1 1 XR3 3400

c PART 1 13 PART 1 14 XR3 3410

c 1 1 XR3 3420

c ********** ********** XR3 3430

c * ABEKI * * ORTHO * XR3 3440

c * BESSEL * ********** XR3 3450

c ********** * XR3 3460

c 1 1 XR3 3470

c 1 1 XR3 3480

r XR3

XR3

3490

C 1 1 l 1 1 1 1 3500

C 1 1 l 1 1 1 1 XR3 3510

C 1 1 l 1 1 1 1 XR3 3520

C 1 1 l 1 1 1 1 XR3 3530

r ********** •kkk-kk'kk'k'k'k ********** 1 ********** 1 ********* *XR3 3540

C * BEJN * * BEZERO * * OANOVA *1 * OCOEFF *1 * OCOVAR *XR3 3550

C * BINTJO * * BEZONE * ********** 1 ********** 1 ********* *XR3 3560

C * DBEJ * * CBEI *

1 1 XR3 3570

C * STRUVE * * CBEK *

1 1 XR3 3580

c ********** * COMELL * 1 1 XR3 3590

C * FOURIA *

1 1 XR3 3600

c ********** ********** ********** XR3 3610

c * OPONE * * ORTPLT * XR3 3620

c ********** *^:******** XR3 3630

c XR3 3640

c XR3 3650

C 3.C APPROXIMATE SIZE OF EACH PART ( 3N THE NBS COMPUTER XR3 3660

C THE SIZE AS INDICATED BELOW ASSUMES THAT A MACHINE WORD XR3 3670

C IS 36 BITS LONG, A FLOATING POINT NUMBER USES ONE MEMORY XR3 3680

C WORD, AND EACH MACHINE INSTRUCTION IS ONE WORD LONG. XR3 3690

C XR3 3700

C I PART 1 44000 WORDS XR3 3710

C XR3 3720

C A. ROUTINES 13000 XR3 3730

c B. LABELED COMMON 31000 (THIS INCLUDES A WORK XR3 3740

c SHEET OF 1250C WORDS) XR3 3750

c C. NOT COUNTED BUT MUST BE ADDED ARE XR3 3760

c LIBRARY FUNCTION ] ROUTINES (I. E. SIN, COS, ETC), XR3 3770

c AND INPUT OUTPUT ROUTINES. (FOR THE NBS UNTVAC XR3 3780 c 1108 SYSTEM THIS IS ABOUT 4500 WORDS] XR3 3790 c XR3 3800

C II PART 2 4100 WORDS XR3 3810

c XR3 3820

C III PART 3 4300 WORDS XR3 3830 c XR3 3840

30
c

c

c

c

c

cc

c

c

cc

cc

c

c

cc

cc

c

c

c

c

c

c

cc

c

cc

c

c

c

c

cc

c

c

c

c

c

c

c

c

cc

c

cc

c

c

c

c

cc

c

c

c

cc

IV PART 4 4000 WORDS XR3 3850

XR3 3860

V PART 5 2900 WORDS XR3 3870

XR3 3880

VI PART 6 3000 WORDS XR3 3890

XR3 3900

VII PART 7 4400 WORDS XR3 3910

(THIS SEGMENT IS WHERE THE USER MAY XR3 3920

ADD HIS SUBROUTINES) XR3 3930

XR3 3940

VIII PART 8 3700 WORDS XR3 3950

XR3 3960

IX PART 9 4000 WORDS XR3 3970

XR3 3980

X PART 10 3200 WORDS XR3 3990

XR3 4000

XI PART 11 3400 WORDS XR3 4010

XR3 4020 XII PART 12 3800 WORDS XR3 4030

XR3 4040 XIII PART 13 3300 WORDS XR3 4050

XR3 4060

XIV PART 14 4400 WORDS XR3 4070

XR3 4080

4.**ASSIGNMENT OF INPUT AND OUTPUT UNITS************************** XR3 4090

XR3 4100

4. A I/O UNITS FOR OMNITAB SYSTEM XR3 4110

XR3 4120

FOUR PHYSICAL I/O UNITS ARE USED BY OMNITAB. THROUGH OUT ALLXR3 4130

THE ROUTINES THESE UNITS ARE REFERRED BY THE FOLLOWING XR3 4140

VARIABLES XR3 4150

THIS UNIT IS USED FOR READING THE OMNITAB XR3 4170

COMMANDS AND DATA. (USUALLY CARD READER) XR3 4180

THE OUTPUT UNIT-RESULTS AND COMMANDS XR3 4190

(USUALLY ON LINE PRINTER) XR3 4200

(62 LINES/PAGE, 120 CHAR/LINE) XR3 4205

3. IPUNCH - THE OUTPUT UNIT FOR GENERATING PUNCH CARDS XR3 4210

4. ISCRAT - A MAGNETIC TAPE, DRUM, OR DISC, WHICHEVER XR3 4220

IS AVAILABLE AND FASTER. (THIS I/O UNIT IS XR3 4230

USED BY OMNITAB.) XR3 4240

XR3 4250

THE I/O VARIABLES ARE ASSIGNED LOGICAL UNITS IN THE XR3 4260

SUBPROGRAM CALLED SETUP. XR3 4270

THEY MAY NEED TO BE REASSIGNED FOR DIFFERENT XR3 4280

CONFIGURATIONS. XR3 4290

XR3 4300

THE PRESENT ASSIGNMENT IS AS FOLLOWS XR3 4310

XR3 4320

INUNIT=5 XR3 4330

IPRINT=6 XR3 4340

IPUNCH=3 XR3 4350

ISCRAT=45 XR3 4360

XR3 4370

4.B PHYSICAL TAPE ASSIGNMENTS FOR OMNITAB TAPE COMMANDS XR3 4380 IN ADDITION USERS MAY REFERENCE 6 TAPE UNITS THROUGH OMNITABXR3 4400

COMMANDS (E.G. AS READ TAPE 'L' , WRITE TAPE 'L' ,ETC, WHERE XR3 4410

'L' IS EITHER A,B,C,D,E OR F) . THE LOGICAL UNITS ASSIGNED TOXR3 4420

THESE TYPES ARE 7 THRU 12, RESPECTIVELY. IF THESE LOGICAL XR3 4430

UNITS HAVE TO BE CHANGED, THE FOLLOWING CORRECTIONS MUST BE XR3 4440

MADE IN THE SUBPROGRAM TAPOP XR3 4450

INUNIT

I PRINT

31
c

c

c

c

c

c

c

c

c

cc

c

c

cc

cc

c

c

c

c

cc

cc

c

c

cc

c

ccc

c

c

c

c

c

c

c

cc

c

c

cc

cc

c

cc

c

cc

c

cc

c

c

c

ITAPE=I+6

MUST BE CHANGED TO

ITAPE=I+NEWTP (WHERE NEWTP IS THE FIRST LOGICAL I/O UNIT

AVAILABLE MINUS 1)

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

5.**DIRECTIONS FOR THE USE OF DUMMY A-F*************************** XR3

XR3

THE OMNI TAB COMMANDS DUMMY A THROUGH DUMMY F ENABLE THE USERXR3

TO EXECUTE HIS OWN SUBPROGRAM(S) THROUGH OMNITAB. THE FIRST XR3

STATEMENT OF THE SUBPROGRAM MUST BE ONE OF THE FOLLOWING XR3

SUBPROGRAM DUMMYA

SUBPROGRAM DUMMYB

SUBPROGRAM DUMMYC

SUBPROGRAM DUMMYD

SUBPROGRAM DUMMYE

SUBPROGRAM DUMMYF

IN WRITING THE SUBPROGRAM(S) THE FOLLOWING INFORMATION

SHOULD BE FOUND USEFUL.

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

THE FOLLOWING IABELED COMMON, DIMENSION AND EQUIVALENCE

STATEMENTS SHOULD BE USED. ******************* XR3

COMMON/BLOCRC/NRC,RC(NSIZRC+100j **SEE XREF3 PARTI** XR3

DIMENSION ARGS(IOO) **FOR DEFINITION ** XR3

EQUIVALENCE ARGS(l) ,RC(NSIZRC+1) **OF NSIZRC ** XR3

COMMON/SCRAT/NR,NS2,A(NSIZRC+1000) ******************* XR3

COMMON/BLOCKD/IARGS(100) ,KIND(100) ,ARGTAB(100) ,NRMAX,

1 NROW,NCOL,NARGS,VWXYZ(8),NERROR

**VMIABLE**

RC

NRC

NROW

NCOL

NRMAX

A

NS

NS2

IARGS

MGS

NARGS

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

XR3

ALL INTEGER ARGUMENTS OF THE DUMMYA, . .

. ,OR XR3

DUMMYF INSTRUCTION XR3

XR3

ALL FLOATING ARGUMENTS OF THE DUMMYA, . .

. , XR3

OR DUMMYF INSTRUCTION XR3

XR3

THE NUMBER OF ARGUMENTS OF THE DUMMYA, . .

. , XR3

OR DUMMYF INSTRUCTION XR3

XR3

**MEANING OF VARIABLE**

THE NAME OF THE WORKSHEET AREA

THE SIZE OF THE WORKSHEET AREA

THE NUMBER OF ROWS IN THE WORKSHEET

THE NUMBER OF COLUMNS IN THE WORKSHEET

THE PRESENT NO. OF ROWS BEING USED

AVAILABLE SCRATCH AREA

SIZE OF THE SCRATCH AREA

ONE-HALF THE SIZE OF THE SCRATCH AREA

4460

4470

4480

4490

4500

4510

4520

4530

5090

5100

5110

5120

5130

5150

5160

5170

5180

5190

5200

5230

5240

5250

5260

5270

5280

5290

5300

5310

5320

5330

5340

5350

5360

5370

5380

5390

5400

5410

5420

5430

5440

5450

5460

5470

5480

5490

5500

5510

5520

5530

5540

5550

5560

5570

5580

5590

5600

5610

5620

5630

32
C KIND KIND(I)=0 IF ITU ARGUMENT IS AN INTEGER XR3 5640

C KIND (I) =1 IF ITH ARGUMENT IS FLOATING XR3 5650

C XR3 5660

C NERROR NERROR=0 IF THERE HAVE BEEN NO FATAL ERRORS XR3 5670

C NERROR= NO. OF PREVIOUS FATAL ERRORS XR3 5680

C XR3 5690

C THE FOLLOWING STATEMENTS SHOULD BE USED FOR CHECKING XR3 5700

C PURPOSES AND TO COMPUTE ADDRESSES OF COLUMNS IN THE XR3 5710

C WORKSHEET. XR3 5720

C XR3 5 7 30 C CALL CKIND(J) XR3 5740

C XR3 5750

C CKIND IS A SUBPROGRAM WHICH CHECKS THE KIND OF ARGUMENTS XR3 5760

C (INTEGER OR FLOATING) .J VALUES ARE CHECKED. XR3 5770

C UPON RETURN FROM THE SUBPROGRAM J IS SET TO THE FOLLOWINGXR3 5780

C J=0 IF ALL ARGUMENTS ARE INTEGER XR3 5790

C J=l IF ALL ARGUMENTS ARE FLOATING XR3 5800

C J=2 IF SOME ARE INTEGER AND SOME FLOATING XR3 5810

C XR3 5820

C CALL ADRESS(I,J) XR3 5830

C XR3 5840

C ADRESS IS A SUBPROGRAM WHICH COMPUTES THE STARTING XR3 5850

C ' ADDRESS IN THE WORKSHEET OF IARGS(I) . XR3 5860

C UPON RETURN FROM THE SUBPROGRAM J IS SET TO THE FOLLOWINGXR3 5870

C J=-(I+NRQ IF ARGUMENT IS A FLOATING POINT NUMBER XR3 5880

C J=0 COLUMN NUMBER IS ILLEGAL XR3 5890

C J= ADDRESS OF A LEGAL COLUMN XR3 5900

C XR3 5910

C CALL CHKCOL(J) XR3 5920

C XR3 5930

C CHKCOL IS A SUBPROGRAM WHICH CHECKS THAT ALL NARGS XR3 5940

C ARGUMENTS ARE LEGAL COLUMN NUMBERS AND CONVERTS THESE XR3 5950

C COLUMN NUMBERS TO THEIR STARTING ADRESSES IN THE WORK- XR3 5960

C SHEET AND STORES THE ADDRESSES IN IARGS(I) , 1=1,... NARGS XR3 5970

C UPON RETURN FROM THE SUBPROGRAM J IS SET TO THE FOLLOWINGXR3 5980

C J=l IF THERE ARE ANY ILLEGAL ARGUMENTS XR3 5990

C J=0 IF ALL ARGUMENTS ARE LEGAL COLUMN NUMBERS XR3 6000

C XR3 6010

C CALL ERROR(I) XR3 6020

C XR3 6030

C IF THERE ARE ERRORS IN THE USER'S CHECKING OF ARGUMENTS, XR3 6040

C ERROR MESSAGES ARE PROVIDED BY CALLING SUBPROGRAM ERROR. XR3 6050

C SEE XREF6 FOR EXISTING ERROR MESSAGES. XR3 6060

C XR3 6070

C IF(NERROR.NE.O) RETURN XR3 6080

C XR3 6090

C THE USER SHOULD NOT PROCEED WITH HIS SUBPROGRAM IF NERRORXR3 6100

C IS NOT ZERO, BECAUSE INFORMATION IN THE WORKSHEET MAY NOTXR3 6110

C BE CORRECT. XR3 6120

C XR3 6130

Q ***********************WApMT|jr**********************ft****** XR3 6140

C XR3 6150

C THE USE OF DUMMYA, . . .DUMMYF IS NOT RECOMMENDED FOR NOVICEXR3 6160

C USERS. THE INFORMATION GIVEN ABOVE IS PRIMARILY FOR THE XR3 6170

C CONVENIENCE OF THE SYSTEM PROGRAMMER. XR3 6080

C XR3 6190

q *********************************************************** XR3 6200

C XR3 6210

END XR3 6220

33
Listing Of XREF4

C XREF4 XR4 10 C VERSION 5.00 XREF4 5/15/70 XR4 20 C XR4 30 C l.**CHANGES THAT MAY BE NEEDED BECAUSE SYSTEM FUNCTIONS ARE USED  XR4 40

C XR4 50 C DIFFERENT COMPUTERS HAVE DIFFERENT SIZE WORDS WHICH AFFECT THE XR4 60 C SIZE OF REAL AND DOUBLE PRECISION NUMBERS. BECAUSE OF THIS, AN XR4 70 C UPPER LIMIT HAS TO BE PLACED FOR COMPUTING THE EXPONENTIAL, SINE,XR4 80

C AND COSINE FUNCTIONS. THESE LIMITS MAY HAVE TO BE CHANGED FOR XR4 90

C SOME COMPUTERS. XR4 100

C XR4 110

C IN ORDER TO MAKE OMNITAB AS SYSTEM FREE AS POSSIBLE, NONE OF THE XR4 120

C PROGRAMS MAKE DIRECT REFERENCE TO THE FOLLOWING FUNCTIONS: XR4 130

C ALOG,ALOG10,COS,DCOS,DEXP,DLOG,DSIN,DSQRT,EXP,SIN,SQRT,TANH, XR4 140

C AND RAISING A VARIABLE TO ANOTHER VARIABLE. XR4 150

C INSTEAD THE OMNITAB SUBPROGRAMS CALL THE FOLLOWING FUNCTIONS: XR4 160

C FLOG,FLOG10,FCOS,FDCOS,FDEXP,FDLOG,FDSIN,FDSQRT,FEXP,FEXP2, XR4 170

C FSIN,FSQRT,FTANH, XR4 180

C WHICH CHECK TO SEE IF ARGUMENT OF THE FUNCTION IS ILLEGAL OR OUT XR4 190 C OF BOUNDS. IF THE FUNCTION CAN NOT BE EVALUATED INFORMATIVE XR4 200

C DIAGNOSTIC IS PRINTED AND VALUE OF ZERO IS RETURNED. IF THE XR4 210

C ARGUMENT IS VALID THE SYSTEM FUNCTIONS ARE USED. XR4 220

C XR4 230 C THE LIMITS ARE DEFINED IN BLOCK DATA SUBPROGRAM WHICH HAS THE XR4 240 C FOLLOWING CARD; XR4 250

C BLOCK DATA LBCONS XR4 260

C THESE VALUES MAY HAVE TO BE REDEFINED IN SOME CASES. XR4 270

C XR4 280

C XTRIG IS NOW SET = 3.3E7 UPPER ABSOLUTE LIMIT FOR SIN, COS SINGLEXR4 290

C PRECISION REAL NUMBER XR4 300

C XR4 310

C XEXP IS NOW SET =88.0 UPPER LIMIT FOR EXP SINGLE PREC. REAL NO. XR4 320 C XR4 330

C DSNCOS IS NOW SET =3.5D16 UPPER ABSOLUTE LIMIT FOR SIN, COS XR4 340

C DOUBLE PRECISION REAL NO. XR4 350 C XR4 360

C DXEXP IS NOW SET = 70 4. DO UPPER LIMIT FOR EXP DOUBLE PREC. REAL NOXR4 370

C XR4 380

C ER IS NOW SET=l.E-8 COMPUTER ZERO. XR4 390

C XR4 400 C NBC IS NOW SET=11 NUMBER OF BINARY BITS IN CHARACTERISTIC OF A XR4 410 C DOUBLE PRECISION NUMBER XR4 420 C XR4 430

C NBM IS NOW SET=60 NUMBER OF BINARY BITS IN MANTISSA OF A DOUBLE XR4 440

C PRECISION NUMBER XR4 450 C XR4 460 C TRRTPI IS NOW SET=1. 12837916 7095512574D0 THE VALUE OF 2.0/SQRT(PI)XR4 470

C TO 19 DECIMAL PLACES XR4 480

C XR4 490

C XR4 500 C XR4 510

XR4 520

34
C I FLOG,FLOG10: IF X IS GREATER THAN ZERO, LOGE OR LOG10 IS COMPUTED XR4 530 C IF X IS ZERO OR NEGATIVE , INFORMATIVE DIAGNOSTIC IS XR4 540 C PRINTED AND ZERO VALUE RETURNED XR4 550 C II FDLOG: SAME AS I EXCEPT RESULT IS DOUBLE PRECISION VALUE XR4 560 C III FCOS,FSIN: IF X IN ABSOLUTE VALUE IS LESS THAN OR EQUAL TO 3.3E7XR4 570

C COS OR SIN IS EVALUATED. OTHERWISE AN INFORMATIVE XR4 580

C DIAGNOSTIC IS PRINTED AND ZERO VALUE RETURNED XR4 590 C IV FDCOS,FDSIN: IF X IN ABSOLUTE VALUE IS LESS THAN OR EQUAL TO XR4 600

C 3.5D16, COS OR SIN IS EVALUATED IN DOUBLE PREC. XR4 610

C OTHERWISE AN INFORMATIVE DIAGNOSTIC IS PRINTED AND XR4 620

C ZERO VALUE RETURNED. XR4 630

C V FEXP: IF X IS LESS OR EQUAL TO 88.0 THE EXP(X) IS COMPUTED XR4 640

C IF VALUE IS NOT EVALUATED, THEN AN INFORMATIVE XR4 650

C DIAGNOSTIC IS PRINTED AND ZERO VALUE RETURNED XR4 660

C VI FDEXP: EXP(X) IS COMPUTED IN DOUBLE PREC, IF X IS NOT XR4 670

C GREATER THAN 704. DO. OTHERWISE ZERO VALUE IS RETURNEDXR4 680

C AND AN INFORMATIVE DIAGNOSTIC IS PRINTED XR4 690

C VII FEXP2: EVALUATES X** A. IF ABSOLUTE VALUE OF A IS XR4 700

C GREATER THAN 60 OR NON INTEGER, FEXP2 CALLS FLOG. XR4 710

CVIII FSQRT,FDSQRT: COMPUTES SQUARE ROOT OF SINGLE AND DOUBLE PREC. XR4 720

C NO. RESPECTIVELY IF X IS EQUAL OR GREATER THAN ZERO.XR4 730 C OTHERWISE AN INFORMATIVE DIAGNOSTIC IS PRINTED AND XR4 740

C ZERO VALUE RETURNED XR4 750 C IX FTANH: IF ABSOLUTE VALUE OF 2*X IS LESS THAN OR EQUAL TO 88 XR4 760

C THE TANH IS COMPUTED. IF NOT, INFORMATIVE DIAGNOSTIC XR4 770 C IS PRINTED AND ZERO RETURNED. XR4 780

C XR4 785 C 2.**CHANGES NEEDED IN BLOCK DATA SUBPROGRAM WHICH DEFINES THE XR4 790

C PHYSICAL CONSTANTS********************************************** XR4 800 C XR4 810 C DUE TO THE DIFFERENCE IN THE SIZE OF WORDS IN DIFFERENT COMPUTERS XR4 820

C THE VALUES OF THE PHYSICAL CONSTANTS MAY NEED TO BE CHANGED. XR4 830

C THESE VALUES ARE DEFINED IN THE BLOCK DATA SUBPROGRAM WHICH HAS XR4 840

C THE FOLLOWING CARD; XR4 850

C BLOCK DATA PHYSIC XR4 860 C XR4 870 C 3.**NUMBER OF SIGNIFICANT DIGITS OF FLOATING POINT NOT EIGHT******** XR4 880 C XR4 890

C CHANGES MUST BE MADE IN SUBPROGRAMS ACCDIG AND LBCONS, IF XR4 900

C NUMBER OF SIGNIFICANT DIGITS OF A FLOATING POINT NUMBER IS XR4 910

C NOT 8. XR4 920

C XR4 930

C LET SIGDIG= NUMBER OF SIGNIFICANT DIGITS FOR FLOATING POINTXR4 940

C NUMBER. XR4 950

C XR4 955

C THEN THE FOLLOWING CHANGE MUST BE MADE IN ACCDIG XR4 960

C XR4 970

C DATA ADMAC/8.0/ XR4 975

C XR4 980

C MUST BE CHANGED TO XR4 985

C XR4 990

C DATA ADMAX/SIGDIG/ (WHERE SIGDIG IS A REAL CONSTANT) XR4 1000

C XR4 1010

C LET JSIGD= NUMBER OF SIGNIFICANT DIGITS FOR FLOATING POINT XR4 1020

C NUMBER XR4 1030

C XR4 1035

C THEN THE FOLLOWING CHANGE MUST BE MADE IN LBCONS. XR4 1040

C XR4 1050

C DATA ISIGD/8/ XR4 1055

C XR4 1060

35
C MUST BE CHANGED TO XR4 1065

C XR4 1070

C DATA ISIGD/JSIGD/ (WHERE JSIGD IS AN INTEGER CONSTANT) XR4 1080

C XR4 1090

C 4.**WARNING FOR COMPUTERS WHERE INTEGER NUMBERS MUST BE LESS THAN XR4 1100

C 10 E9 OR ?**32 ************************************************* XR4 1110

C XR4 1120

C CHANGES MUST BE MADE IN SUBPROGRAMS DHRND AND RPRINT, IF XR4 1130

C INTEGER NUMBERS ARE LESS THAN 2**32. XR4 1140

C XR4 1150

C CHANGES INVOLVE SCALING INTEGERS. XR4 1160

C XR4 1170

END XR4 1180

36
Listing Of XREF5

C XREF5 XR5 10 C VERSION 5.00 XREF5 5/15/70 XR5 20

C XR5 30 C l.A VARIABLES DEFINED IN SUBPROGRAM BLOCK********************** XR5 40 C XR5 50 C L IS DEFINED IN LABELED COMMON ABCDEF. XR5 60

C IT CONTAINS THE ANSI FORTRAN CHARACTER SET: XR5 70 C THE LETTERS A-Z XR5 80 C THE NUMBERS 0-9 XR5 90

C THE CHARACTERS =+-*/(),.$ AND BLANK XR5 100

C THE CHARACTER ' IS ALSO INCLUDED XR5 110 C XR5 120 C NOCARD IS DEFINED IN LABELED COMMON FMAT. XR5 130

C IT INITIALLY CONTAINS THE WORD OMNITAB XR5 140

C XR5 150

C ITS CONTENTS ARE CHANGED WHEN AN OMNITAB CARD IS XR5 160

C ENCOUNTERED DURING EXECUTION. XR5 170 C XR5 180

C IFMTS IS DEFINED IN LABELED COMMON FMAT. XR5 190

C IT CONTAINS THE STANDARD OMNITAB OUTPUT FORMAT. XR5 200

C IFMTS(1)=1H( IFMTS(2)=2H1P IFMTS(3)=2H8E XR5 210 C IFMTS(4)=3H15. IFMTS(5)=1H6 IFMTS(6)=1H) XR5 220 C XR5 230 C XR5 240

C l.B *VARIABLES DEFINED IN SUBPROGRAM LBG3NS********************* XR5 250

C XR5 260 C XTRIG IS DEFINED IN LABELED COMMON CONSLB. XR5 270

C XTRIG=3.3E7 AND IS USED BY FSIN AND FCOS SUBPROGRAMS. XR5 280

C SEE XREF4 FOR ADDITIONAL NOTES. XR5 290

C XR5 300

C XEXP IS DEFINED IN LABELED COMMON CONSLB. XR5 310

C XEXP=88.0 AND IS USED BY FEXP SUBPROGRAM. XR5 320

C SEE XREF4 FOR ADDITIONAL NOTES. XR5 330

C XR5 340

C DSNCOS IS DEFINED IN LABELED COMMON DCONLB. XR5 350

C DSNCOS=3.5D16 AND IS USED BY FDSIN AND FDCOS SUBPROG XR5 360

C SEE XREF4 FOR ADDITIONAL NOTES XR5 370

C XR5 380

C DXEXP IS DEFINED IN LABELED COMMON DCONLB. XR5 390

C DXEXP=704.0D0 AND IS USED BY FDEXP. XR5 400

C SEE XREF4 FOR ADDITIONAL NOTES. XR5 410

C XR5 420

C ER IS DEFINED IN LABELED COMMON CONLB2. XR5 430

C ER=l.E-8 AND IS USED BY CSPINV. XR5 440 C SEE XREF4 FOR ADDITIONAL NOTES XR5 450

C XR5 460

C NBC IS DEFINED IN LABELED COMMON DCONL2. XR5 470 C NBOll AND IS USED BY ERRINT. XR5 480

C SEE XREF4 FOR ADDITIONAL NOTES XR5 490 C XR5 500

C NBM IS DEFINED IN LABELED COMMON DCONL2. XR5 510

C NBM=60 AND IS USED BY ERRINT. XR5 520

37
c

c

c

c

c

c

cc

cc

cc

c

c

c

cc

c

c

c

c

c

c

cc

c

c

c

c

c

c

c

c

c

c

c

c

c

cc

c

c

cc

cc

cc

c

cc

cc

c

SEE XREF4 FOR ADDITIONAL NOTES XR5

XR5

TRRTPI IS DEFINED IN LABELED COMMON DCONL2. XR5

TRRTPI=1.128379167095512574D0 AND IS USED BY ERRINT XR5

SEE XREF4 FOR ADDITIONAL NOTES XR5

XR5

ISIGD IS DEFINED IN LABELED COMMON CONLB2. XR5

ISIGD=8 AND IS USED BY FIXFLO AND INFERR XR5

SEE XREF4 FOR ADDITIONAL NOTES XR5

XR5

THE FOLLOWING VARIABLES ARE DEFINED IN LABELED COMMON CONSTS.XR5

XR5

PI

E

HALFPI

DEG

RAD

XALOG

PI=3. 14159265 (THE VALUE OF PI)

E=2. 71821818 (BASE OF NATURAL LOGS)

HALFPI=1. 5707963 (VALUE OF PI/2)

DEG=57. 2957795 (NUMBER OF DEGREES IN ONE RADIAN)

RAD=. 0174532925 (NUMBER OF RADIANS IN ONE DEGREE)

XALOG=38. (EXPONENT BOUND)

***NOTE*** PI AND E ARE ALSO DEFINED IN PHYSIC SUBPROGRAM

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

XR5

l.C VARIABLES DEFINED IN SUBPROGRAM LOOKTB********************* XR5

XR5

**THE FOLLOWING VARIABLES ARE DEFINED BY LABELED COMMON CODE** XR5

XR5

IALPHA XR5

CONTAINS THE OMNITAB NUMERIC REPRESENTATION FOR THE XR5

CHARACTERS A,B,C,D,E, AND F XR5

THESE CHARACTERS ARE USED BY THE INPUT OUTPUT COMMANDXR5

XR5

NALPHA XR5

CONTAINS THE OMNITAB NUMERIC REPRESENTATION FOR THE XR5

CHARACTERS V,W,X,Y, AND Z XR5

THESE CHARACTERS ARE USED BY THE RESET COMMAND XR5

XR5

XR5

l.D *VARIABLES DEFINED IN SUBPROGRAM PHYSIC********************* XR5

XR5

**THE LABELED COMMON PCONST DEFINES THE FOLLOWING VARIBLES** XR5

P (I) XR5

CONTAINS THE NUMERICAL VALUES FOR THE PHYSICAL XR5

CONSTANTS IN SI AND CGS UNITS. XR5

N (I) XR5

CONTAINS THE OMNITAB NUMERIC REPRESENTATION OF THE XR5

PHYSICAL CONSTANT NAMES USED IN OMNITAB. XR5

XR5

***NOTE*** SEE SUBPROGRAM PHYSIC FOR ADDITIONAL DETAILS. XR5

END XR5

530

540

550

560

570

580

590

600

610

620

630

640

650

660

670

680

690

700

710

720

730

740

750

760

770

780

790

800

810

820

830

840

850

860

870

880

890

900

910

920

925

930

935

940

950

960

970

980

990

1000

1010

1020

1025

1030

1040

38
Listing Of XREF6

c XREF6 XR6 10 c VERSION 5.00 XREF6 5/15/70 XR6 20

c XR6 30 c THIS SBC Nov n\)Vi\!\s \]j ommt.vk ERROR MESSAGE XR6 40 c XR6 50 c XR6 60

c 1 A *A*AAAAAAA*pA'pAT ERROR fvfESSAGES***************************** XR6 70 c XR6 80 c ERROR (1) XR6 90

c NAME NOT FOUND IN LIBRARY XR6 100

c ERROR (2) XR6 110

c ILLEGAL STATEMENT NUMBER XR6 120 c ERROR (3) XR6 130 c ILLEGAL ARGUMENT ON CARD XR6 140

c ERROR (5) XR6 150

c COMMAND NOT ALLOWED IN REPEAT MODE XR6 160

c ERROR (6) XR6 170

c STATEMENT NUMBER MAY NOT BEGIN ANY CARD BETWEEN BEGIN XR6 180

c AND FINISH CARDS XR6 190 c ERROR (7) XR6 200

c ILLEGAL *STATEMENT* XR6 210

c ERROR (9) XR6 240

c NRMAX=0 XR6 250

c ERROR (10) XR6 260

c I4,34H IS AN ILLEGAL NUMBER OF ARGUMENTS XR6 270

c ERROR (11) XR6 280 c COLUMN NUMBER TOO BIG OR LESS THAN 1 XR6 290

c ERROR (12) XR6 300

c COMMAND STORAGE AREA OVERFLOW XR6 310

c ERROR (13) XR6 320

c STATEMENT NUMBER NOT FOUND XR6 330

c ERROR (15) XR6 360

c DIMENSIONED AREA EXCEEDS LIMIT XR6 370

c ERROR (16) XR6 380

c ILLEGAL SIZE ROW NUMBER XR6 390

c ERROR (17) XR6 400

c DEFINED MATRIX OVERFLOWS WORKSHEET XR6 410

c ERROR (18) XR6 420

c INTEGER ARGUMENT LESS THAN -8191 XR6 430

c ERROR (19) XR6 440

c STORED PERFORM STATEMENT WILL EXECUTE ITSELF XR6 450

c ERROR (20) XR6 460

c IMPROPER TYPE OF ARGUMENT XR6 470

c ERROR (21) XR6 480

c COMMAND MUST BE STORED XR6 490

c ERROR (22) XR6 500

c MATRIX IS (NEARLY) SINGULAR XR6 510 c ERROR (23) XR6 520

c INSUFFICIENT SCRATCH AREA XR6 530 c ERROR (24) XR6 540

c DEGREE IS LARGER THAN NO. OF NON-ZERO WEIGHTS XR6 550

c ERROR (25) XR6 560

39
C NEGATIVE WEIGHTS MAY NOT BE USED XR6 570

C ERROR (26) XR6 580 C NUMBER OF COLUMNS IS GREATER THAN NUMBER OF ROWS XR6 590

C ERROR (27) XR6 600

C FORMAT NOT FOUND XR6 610

C ERROR (28) XR6 620

C INCORRECT TAPE UNIT. COMMAND IS NOT EXECUTED. XR6 630

C ERROR (29) XR6 640

C NUMBER OF ARGUMENTS SHOULD BE, 1

2

XR6 650

C ERROR (30) XR6 660

C AN INCREMENT COMMAND CAN NOT INCREMENT ITSELF. XR6 670

C XR6 680

c 1#B ***********aritHMETIC ERROR MESSAGES************************ XR6 690

C XR6 700

C ERROR (101) XR6 710

C NEGATIVE ARGUMENT TO SORT, LOG OR RAISE XR6 720

C ERROR (102) XR6 730

C EVALUATION OF EXPONENT PRODUCES OVERFLOW XR6 740 C ERROR (103) XR6 750

C ARGUMENT OUT OF BOUNDS TO INVERSE FUNCTION XR6 760

C ERROR (104) XR6 770

C ARGUMENT TOO LARGE FOR SIN OR COS, ZERO. RETURNED , I 4 , 6H XR6 780

C TIMES XR6 790

C ERROR (105) XR6 800

C BESSEL ARGUMENTS SCALED TO AVOID OVER/UNDER FLOW. XR6 810

C RETURNED, I 4, 6H TIMES XR6 820

C ERROR (106) XR6 830

C DIVISION BY ZERO, RESULT SET=0,I4,6H TIMES XR6 840

C ERROR (107) XR6 850

C TRIG FUNCTION NOT DEFINED RESULTS SET=0,I4,6H TIMES XR6 855

C ERROR (108) XR6 860

C ONE OF THE VALUES COMPARED IS ZERO, ABSOLUTE TOLERANCE XR6 865

C WAS USED, 14, 6H TIMES XR6 870

C ERROR (109) XR6 873

C X FOR ELLIPTICAL INTEGRALS IS = 1.0 OR GREATER. RESULT ISXR6 875

C SET TO 0.0.,I4,6H TIMES XR6 877

C XR6 879

C i.c *********** INpoRMATIVE DIAGNOSTIC MESSAGES****************** XR6 880

C XR6 890

C ERROR (201) XR6 900

C TOO MUCH DATA IN SET, READ OR GENERATE, SPILL LOST XR6 910

C ERROR (202) XR6 920

C COMMAND NOT ALLOWED IN REPEAT MODE. EXECUTED BUT NOT XR6 930

C STORED XR6 940

C ERROR (203) XR6 950

C VALUE REQUESTED IN SHORTEN, ACOALESCE OR AAVERAGE NOT XR6 960

C FOUND. XR6 9 70 C ERROR (204) XR6 980

C BAD HEAD. COLUMN GT 50 OR NO / XR6 990

C ERROR (205) XR6 1000

C THIS COMMAND WAS NOT EXECUTED BECAUSE ITS MEANING WAS XR6 1010

C QUESTIONABLE XR6 1020

C ERROR (206) XR6 1030

C F LESS THAN , SET = XR6 1040

C ERROR (207) XR6 1050

C NU1 OR NU2 LESS THAN 1 XR6 1060

C ERROR (208) XR6 1070

C NU1 OR NU2 TRUNCATED TO INTEGER - XR6 1080

C ERROR (209) XR6 1090

C IMPROPER TITLE NUMBER, ASSUMED 1 XR6 1100

40
c ERROR ( c

f-1 ERROR ( cc ERROR ( c

c ERROR ( c

c ERROR ( c

c ERROR ( c

c ERROR ( cc

c ERROR ( c

c ERROR ( c

c

c ERROR ( c

c ERROR ( cc

c ERROR ( c

c ERROR ( c

c ERROR f cc ERROR ( c

c ERROR ( c

cc ERROR ( c

c

c ERROR ( c

cc ERROR ( c

cc ERROR ( c

c ERROR ( cc

c ERROR ( cc ERROR ( c

c

c ERROR ( c

c ERROR ( c

c

210) XR6 1110

NO OF ROWS NOT = TO COLS. MATRIX USED LARGEST SQUARE XR6 1120

211) XR6 1130

ASTERISK STRING IMPLYING 'THRU' INCORRECT, IGNORED XR6 1140

212) XR6 1150

UNNECESSARY ARGUMENTS IN COMMAND IGNORED XR6 1160

213) XR6 1170

PARTIAL STORAGE OF MATRIX XR6 1175

214) XR6 1180

INSUFFICIENT SCRATCH AREA XR6 1185

215) XR6 1190

NRMAX IS NOT LARGE ENOUGH TO ALLOW ITERATION XR6 1200

216) XR6 1210

1ST COLUMN OF ISETUP OR ISOLATE IS NOT MONOTONIC OR IS XR6 1220

CONSTANT. XR6 1230

217) XR6 1240

ITERATION HAS FOUND NO VALUES. XR6 1250

218) XR6 1260

WORKSHEET IS TOO SHORT TO ACCOMODATE ALL THE VALUES XR6 1270

GENERATED BY THIS COMMAND. XR6 1280

219) XR6 1290

MAXMIN HAS FOUND NO EXTREMA. XR6 1300

220) XR6 1310

MAXMIN HAS FOUND AND IGNORED A TRIAD OF X'S WITH AT LEASTXR6 1320

TWO IDENTICAL VALUES. XR6 1330

221) XR6 1340

MORE THAN ONE ARGUMENT IN COMMAND. ONLY FIRST ONE IS USEDXR6 1350

222) XR6 1360

FORMAT NOT FOUND. READABLE FORMAT IS USED XR6 1370

223) XR6 1380

ONE, SOME OR ALL WEIGHTS ARE NEGATIVE XR6 1390

224) XR6 1400

ALL WEIGHTS ARE ZERO. COMMAND IS NOT EXECUTED XR6 1410

225) XR6 1420

ARG FOR BESIN,BESJN,BESKN GIVES A RESULT TOO LARGE/SMALL . XR6 1430

COMMAND NOT EXECUTED. XR6 1440

226) XR6 1450

COLUMN NOT LONG ENOUGH TO STORE ALL ELEMENTS. ONLY NROW XR6 1460

WILL BE STORED. XR6 1470

227) XR6 1480

NOT ENOUGH DATA ON COL TO RESTORE MATRIX/ARRAY. DATA XR6 1490

AVAILABLE WILL BE USED. XR6 1500

228) XR6 1510

SUM OF SQRS DO NOT ADD UP-ABS. VALUE OF (TOTAL -ROW- COL- XR6 1520

RES.) /TOTAL EXCEEDS 5.E-7 ) XR6 1530

229) XR6 1540

MORE THAN 50 HEAD COLUMN COMMANDS HAVE BEEN USED. XR6 1550

230) XR6 1560

ATTEMPT TO PROMOTE FROM BELOW NRMAX. FIRST ARGUMENT IS XR6 1570

RESET TO NRMAX. XR6 1580

231) XR6 1590

ATTEMPT TO DEMOTE OFF THE WORKSHEET. SPILL IS LOST. XR6 1600

233) XR6 1640

NEGATIVE VALUE (S) WERE ENCOUNTERED BY PARTITION FUNCTION. XR6 1650

ZEROES STORED. XR6 1660

234) XR6 1670

NEGATIVE ABSOLUTE TEMPERATURES CONVERTED. XR6 1680

235) XR6 1690

CAUTION, USE EXPERIMENTALLY ONLY. NOT OPTIMUM IN ORDER XR6 1700

TO MAKE IT MACHINE INDEPENDENT. REFERENCES - J.B.KRUSKAL,XR6 1710

41
C ACM, 12, 92. AND J.H. HALTON,SIAM REV., 12,1. XR6 1720

C ERROR (236) XR6 1730

C COMMAND IGNORED - S BEFORE COMMAND NAME MEANINGLESS IF XR6 1740

C NO STORAGE REQUESTED. XR6 1750

C ERROR (237) XR6 1760

C NUMBER OF SIGNIFICANT DIGITS AFTER DECIMAL PT HAS BEEN XR6 1770

C SET TO, 13. XR6 1780

END XR6 1790

42
Latest developments in the subject area of this publication, as well as

in other areas where the National Bureau of Standards is active, are

reported in the I\IBS Technical News Bulletin. See following page.
HOW TO KEEP ABREAST OF MBS ACTIVITIES

Your purchase of this publication indicates an interest in the research,

development, technology, or service activities of the National Bureau of Standards.

The best source of current awareness in your specific area, as well as in other NBS programs of possible interest, is the TECHNICAL NEWS BULLETIN,

a monthly magazine designed for engineers, chemists, physicists, research

and product development managers, librarians, and company executives.

If you do not now receive the TECHNICAL NEWS BULLETIN and would like to subscribe, and/or to review some recent issues, please fill out and return the

form below.

Mail to: Office of Technical Information and Publications

National Bureau of Standards

Washington, D. C. 20234

Name

Affiliation

Address

City State Zip

| | Please send complimentary past issues of the Technical News Bulletin.

I | Please enter my 1 — yr subscription. Enclosed is my check or money

order for $3.00 (additional $1.00 for foreign mailing).

Check is made payable to: SUPERINTENDENT OF DOCUMENTS.

TO 550


NBS TECHNICAL PUBLICATIONS

PERIODICALS NONPERIODICALS

JOURNAL OF RESEARCH reports National

Bureau of Standards research and development in physics, mathematics, chemistry, and engineering.

Comprehensive scientific papers give complete details of the work, including laboratory data, experimental

procedures, and theoretical and mathematical analyses. Illustrated with photographs, drawings, and

charts.

Published in three sections, available separately: • Physics and Chemistry

Papers of interest primarily to scientists working in these fields. This section covers a broad range of physical and chemical research, with major emphasis on standards of physical measurement, fundamental

constants, and properties of matter. Issued six times a year. Annual subscription: Domestic, $9.50; for- eign, $11.75*. • Mathematical Sciences

Studies and compilations designed mainly for the mathematician and theoretical physicist. Topics in mathematical statistics, theory of experiment design,

numerical analysis, theoretical physics and chemistry, logical design and programming of computers

and computer systems. Short numerical tables. Issued quarterly. Annual subscription: Domestic,

$5.00; foreign, $6.25*. • Engineering and Instrumentation

Reporting results of interest chiefly to the engineer

and the applied scientist. This section includes many

of the new developments in instrumentation resulting

from the Bureau's work in physical measurement,

data processing, and development of test methods.

It will also cover some of the work in acoustics,

applied mechanics, building research, and cryogenic

engineering. Issued quarterly. Annual subscription:

Domestic, $5.00; foreign, $6.25*.

TECHNICAL NEWS BULLETIN

The best single source of information concerning the Bureau's research, developmental, cooperative and

publication activities, this monthly publication is designed for the industry-oriented individual whose

daily work involves intimate contact with science and

technology

for engineers, chemists, physicists, re- search managers, product-development managers, and

company executives. Annual subscription: Domestic,

$3.00; foreign, $4.00*.

• Difference in price is due to extra cost of foreign mailing.

Applied Mathematics Series. Mathematical tables, manuals, and studies. Building Science Series. Research results, test methods, and performance criteria of building materials, components, systems, and structures. Handbooks. Recommended codes of engineering

and industrial practice (including safety codes) de- veloped in cooperation with interested industries,

professional organizations, and regulatory bodies. Special Publications. Proceedings of NBS conferences, bibliographies, annual reports, wall charts,

pamphlets, etc. Monographs. Major contributions to the technical

literature on various subjects related to the Bureau's

scientific and technical activities. National Standard Reference Data Series.

NSRDS provides quantitative data on the physical

and chemical properties of materials, compiled from

the world's literature and critically evaluated. Product Standards. Provide requirements for sizes, types, quality and methods for testing various industrial products. These standards are developed cooperatively with interested Government and industry groups

and provide the basis for common understanding of product characteristics for both buyers and sellers. Their use is voluntary.

Technical Notes. This series consists of communications and reports (covering both other agency and

NBS-sponsored work) of limited or transitory interest. Federal Information Processing Standards Publications. This series is the official publication within

the Federal Government for information on standards

adopted and promulgated under the Public Law

89-306, and Bureau of the Budget Circular A-86

entitled, Standardization of Data Elements and Codes

in Data Systems.

Order NBS publications from: Superintendent of Documents

Government Printing Office

Washington, D.C. 20402
U.S. DEPARTMENT OF COMMERCE

WASHINGTON, D.C. 20230

OFFICIAL BUSINESS POSTAGE AND FEES PAID

U.S. DEPARTMENT OF COMMERCE

PENALTY FOR PRIVATE USE, $300

1 H