Search Documents

By J. K. Brooke, F. A. Renninger

During 1966, crushed stone production in the United States totaled just over 811 million tons valued at almost $1.2 billion. This represented in- creases of 4 % in tonnage and 5 % in value over that for 1965. Excluding the crushed stone used for agricultural and certain industrial purposes (cement, lime and flux), crushed stone mineral aggregate production totaled 616.3 million tons valued at almost $910 million, also up roughly 4% in tonnage and 5% in value over the 1965 figures. Limestone was by far the most widely used rock variety, followed by basalt, granite and sandstone in that order. Concrete aggregate and roadstone continued as the main uses for stone, accounting for 63% of all tonnage produced in 1966 and 60% of total stone value. Sand and gravel production in 1966 increased 3% over that for 1965. Total 1966 output was just in excess of 934 million tons valued at almost $985 million as compared with 908 million tons at $957 million for the previous year. Excluding the higher valued industrial sands, 903.2 million tons of sand and gravel valued at $891.6 million were used directly as mineral aggregates, representing a 2.9% increase in tonnage and a 2.4% increase in value over 1965.

Jan 2, 1968

By F. F. Aplan

The past year has been another year of vigorous fundamental research in mineral processing as evidenced by the active publication of results, particularly in Germany, Japan, United Kingdom, United States and USSR. It was also noteworthy for several significant meetings and conferences. The Basic Science sessions at the AIME Annual Meeting continued to attract an increasing audience and the sessions were especially notable for the active discussion of papers. The third Engineering Foundation Conference on Particulate Matter Systems,1 chaired by F. F. Aplan, was held in Milwaukee in August and attracted over 100 experts from five continents. Sessions at this week-long conference covered comminution, particle characterization, separation of fine particles, fluid-solid separations as applied to pollution control, flotation kinetics, particle-particle interactions and particle agglomeration. This popular conference will be continued under John Karpinski of Bethlehem Steel Corp. as 1968 chairman.

Jan 2, 1968

By N. E. Grosvenor

The study of rock mechanics is becoming more important each year and several groups are currently engaged in the study of rock mechanics instrumentation and techniques that can be used to determine rock properties in place. As it becomes necessary to mine lower grade ores, larger equipment and tonnages are required, which in turn call for larger openings underground. In addition, rock mechanics studies are necessary for designing large underground openings for defense and powerhouses, and highway and water tunnels. Meramec Mining Co. has a rock mechanics program in which they are trying to develop a means of detecting ground instability prior to failure. In one case advance warning resulted in valuable lead time for back filling stopes before pillar failure had progressed to a dangerous extent.1

Jan 2, 1968

By J. D. Corbett

Mining geophysics in 1967 continued to expand with greater application in the exploration and development of mineral deposits. Moreover, geophysical activity increased without the particular emphasis of a Pine Point or Timmins discovery to spur interest. This trend is expected to continue through 1968 as geophysics more and more becomes an integral part of a mineral exploration program. This past year saw the "little black box" of the old-fashioned "doodlebugger" pushed further into the background, replaced by increasingly sensitive equipment of sophisticated design. The trend also is away from the one man geophysical operation wherein the same individual is responsible for the development, design and operation of the field gear, the planning and execution of a survey and analyzing and interpreting the data. In today's technology a staff of specialists, electronics engineers, mathematicians, physicists-to name but a few- are required to back up the exploration geophysicist in the development of new methods and techniques, in the reduction of the voluminous amount of field data and in the interpretation of the anomalies which can be enhanced by various high-speed computer techniques. Obviously a full back-up staff of these varied abilities is impractical to most exploration companies. However, the equivalent function is provided to the lone staff geophysicist by the many competent geophysical contractors and individual consultants.

Jan 2, 1968

By T. A. Simpson

The search for ore continued at its relentless pace throughout 1967. Canada, South Africa and Australia plus a few scattered localities on the globe reported minerals finds of significant importance. Many of the reported finds were not high grade ore, but they nonetheless will assist in sup- plying a highly industrialized, mineral hungry world. No longer does inaccessibility and remoteness close an area to mineral exploration and development. The communication and transportation systems of today have closed the gap and, moreover, the demand for essential mineral commodities is causing price shifts which encourage exploitation of deposits heretofore considered uneconomical.

Jan 2, 1968

By Wayne S. Cavender

During the past year, geochemical prospecting appears to have come of age as an exploration method, and its acceptance by the mining industry is widespread. There is a growing recognition that applied geochemistry is a scientific tool and must be performed and interpreted by scientific personnel who are trained in its utilization. At present, most personnel training is done by geochemical departments in government organizations and industry groups. Most schools in the United States and Canada do not offer geochemical training, particularly at graduate levels, due to faculty limitations, lack of budgetary support and the erroneous impression that applied geochemistry is not sufficiently academic. As a result, some companies have found it expedient to send selected personnel to European schools for advanced training. A critical personnel shortage now exists will probably continue into the near future. Indications are that the rate of development of new geochemical techniques has levelled off, and the probability of new breakthroughs is not high. Present research is essentially limited to efforts by government agencies and a few private companies. Progress has been the most rapid in the development of analytical techniques, but the need remains for investigation into sampling technology and the development of more quantitative methods for the interpretation of geochemical results. Geochemists need a better understanding of the numerous physical and chemical factors which affect and control primary and secondary dispersion phenomena in all geologic environments. Detailed studies must be made of the relationship between geologic, geochemical and geophysical factors. With better interpretive techniques, the return on the geochemical exploration dollar will increase.

Jan 2, 1968

By M. E. Prunty

Activity in coal preparation in 1967 was on the up-swing and all indications are that demands will be greater in the future. Feed tonnages to plants are increasing and a significant factor is that the top size being fed to plants is decreasing. This has been brought about by market demands, better cokability in some cases and emphasis on sulfur reduction in others. New plants are being built and some are in operation in the coal fields of southern Illinois. These plants are producing metallurgical coal in circuits that contain heavy medium cyclones and froth flotation on the - ½ mm x 0.

Jan 2, 1968

By E. P. Pfleide, L. W. Gibbs. J. R. Gross

A method using computer techniques is described for the comparative evaluation of truck performance over any given haul road using readily obtainable manufacturers' data. An additional computer program is utilized to calculate individual cycle times for each truck and shovel combination, fleet requirements for any desired production rate, and to estimate comparative costs. The introduction of a large number of variables into the program permits management to determine the critical system factors having the greatest bearing on costs. Thus, the method described can be used not only for analysis of needs for a new surface mine but for cost improvement in existing mines where changes in truck and shovel operations are contemplated. The systems analysis approach to the solution of problems in the mining industry is gaining in industry acceptance. This is a logical development in view of the fact that a mine constitutes one of the most dynamic systems of production. In contrast with a manufacturing plant, situated in a location selected through choice and utilizing a regulated source of quality controlled raw materials, the miner must go where the ore is and he is faced with ever changing physical conditions. He is subjected to the vagaries of nature, and he can never mine the same ore twice. What is systems engineering? There seems to be no concise and universally accepted definition of the term and its meaning is still fairly local. It deals largely with trade-offs of subsystem characteristics to achieve some desired result in an overall system. This concept of suboptimization can be explained by an example. If the objective of a mining company is to produce a certain quality ore at minimum cost, the optimum operation is not achieved by minimizing all costs, since all phases of mining are interrelated. Reduction of blasting costs might well increase loading, hauling, and crushing costs out of proportion to the resulting saving in explosives. Hence tradeoffs are necessary. It is evident then that an optimum system is achieved through a blending of suboptimal components which taken together yield the desired result. OBJECTIVE For our mathematical model we will assume the existence of a hypothetical new organization called Taconite Mining Co. (TMC). TMC has made an analysis of their orebody, their proposed operation, and

Jan 1, 1968

By Robert M. Grogan, James C. Bradbury

The Illinois-Kentucky mining district has, since 1880, accounted for 80 per cent of all U.S. production of fluorspar. The ore deposits are of two types: vein deposits formed by fissure fillings along faults and bedding-replacement deposits formed by the replacement of strata at certain favored stratigraphic positions along minor faults and fractures. Fluorite is the chief valuable mineral in the deposits. Sphalerite and galena commonly are present in small amounts but, in some ore bodies, may form a substantial part of the valuable constituents of the ore. Calcite is the chief gangue mineral; quartz and barite Exposed strata range from Devonian are present in various amounts. through Pennsylvanian in age. Igneous intrusive rocks occur as mafic dikes and sills and intrusive breccia dikes and plugs. Structurally, the district is situated on a generally northwest- trending arch that has been sliced into a series of long, narrow blocks by normal faults that trend northeastward. The structural The localization of the ore deposits by faults and fractures, the association of fluorite mineralization with intrusive breccias, and various lines of petrographic and geochemical evidence indicate that the fluorite-zinc-lead deposits were epigenetic and deposited by solutions emanating from a deep alkalic magma.

Jan 1, 1968

By Arnold S. Epstein

Tlze green electroluminescence occurring at liquid-nitvogen temperature in epitaxial gallium phosphide diodes is exarnined using the donor dopants silicon and sulfur. Zinc is used as the p-type diffusant. Two tt~ain peaks are found, namely til 0.56 and 0.535 p. The intensity of the gYeen peak occurring at 0.56 is found to oary as -4"' wlzere Jf is the forward current density and to first order appears to be dependent on the concet~tvation of ionized donovs. Tile vadintive lifetinze fov this process is -10-'sec. The intensity of the 0.56-p lu~ninescence is nt least several ovdevs of magnitude greate,, at liquid-nitvogen temperature than that noted jor the 0.538 radiation for current densities up to about 35 amp per sq cm. The variation oj- the intensity of. the 0.538-p peak with current density caries according to the donor dopant. With sulfur as the donor one has a J12 dependence, whereas when silicon is the donov a J1.5 dependence is noted. Tl~e above can be fitted to a In molecular process (in first order) such as a band to band recombination. The constant B fir band to band vecombination) is found from our experitt~etlt to be of the order of 2 CZI cm per sec, zcl~ereas for tlze impurity vecombination peak (0.56 p]) tile A cotzstor~t is of the order of. 7 x 10-' cu cm per sec. The structuve o/. the diodes at liquid-nitvog-en tett~peruture is oj tlze double injection type, and by carefi~lly selecting. the thickness of the structure one can regulate the onset of the 0.538-p radiation. The diode clzavaclevistics follow an I-Vx laiv where x is >1. Llregatiue resistance in these structures con be correlated with saturation of tlze 0.56- peak. ThE development of high-purity, large single crystals of epitaxially grown gallium phosphide is the ultimate aim of the studies to be reported. While a number of papers and reviews' have been presented on properties of gallium phosphide grown from gallium solution and from the melt, less attention has been paid to the properties of epitaxial gallium phosphide grown from the vapor. In the studies to be reported we have examined the effects of shallow donor impurities, primarily silicon and sulfur, at liquid-nitrogen temperature, on both the electrical and green electroluminescent properties of carefully characterized single-crystal epitaxial Gap. I) EXPERIMENTAL PROCEDURES A) Material Preparation and Properties. The samules used in this investigation were single-crystal n-type epitaxial gallium phosphide. he single crystals were synthesized from the elements by vapor transport and epitaxial deposition on single-crystal gallium arsenide substrates. Typical transport agents used were chlorides. The growth procedures are described in more detail in Ref. 2. The samples were grown in the (111) orientation and had thicknesses up to 480 . Following the growth, the GaAs substrate was removed from the epitaxial Gap by etching in nitric acid. The Hall effect and resistivity of these n-type samples were measured as a function of temperature and analyzed for impurity content and compared with spectroscopic emission analysis. The results are given in greater detail in Ref. 3 but a summary is noted in Table I. In Table I, AED denotes the donor ionization energy, Nd is the donor impurity concentration, lvA is the acceptor impurity concentration, NI is the total impurity concentration, Rip vs temperature is the measured Hall mobility-temperature relationship from liquid nitrogen to about 400"K, n(77"K) is the measured carrier concentration at liquid nitrogen, and R/p(77'~) is the measured Hall mobility at liquid nitrogen. Good agreement is noted between electrical and emission analysis for the major donor impurity, silicon, in the two materials where it is a major constituent. In the case of the two samples SC90-6 and AP138-3(111) B which were sulfur-doped the ionization energy found from the electrical data is in good agreement with that reported for sulfur.4 The acceptor impurities are not known, although from the presence of magnesium in the samples as indicated by emission analysis and its known role as an acceptor in GaP,4 magnesium could be present as an acceptor. Similarly, since silicon is known to be an amphoteric impurity4 it may also be present as an acceptor in the Gap.* The *A small band has been noted at 0.625 11 in samples SC90-6, SC76-4, and SC77-4, suggesting according to Ref. 5 that sillcon could be present as an acceptor In these samples since these authors have found this peak to represent a donor-silicon emission. presence of calcium deduced from emission analysis is not correlatable with the electrical analysis and its presence as an electrically active impurity should be questioned. B) Sample Preparation and Diode Fabrication. The original slices were of the order of 2 sq cm in area witi the (111)B face up. The material was cleaved and sectioned and the best pieces selected (i.e., those which visually appeared to have the least strain and defects). The samples were polished, lightly etched, and cleaned in trichloroethylene. Zinc was diffused in a closed tube at 800°C to form a p-n junction. Diodes were made from these samples. The depth of junction varied, but did not exceed 10 p depth. The back side of the slice was lapped back and a low-resistance n-type contact was formed (at 600°C in forming gas) to an Au-Mo tab with a hole in the middle for

Jan 1, 1968

By Thomas O. Ziebold, Robert E. Ogilvie

Experimental analysis of diffusion samples covering the entire Cu-Ag-Au system at 725°C has been carried out. Experimental coefficients are computed at twenty-eight points in the lermary system. It is found that the direct coefficienl for copper is nearly independent of the silver, content for 1ow-silver alloys and that the cross coefficient for silver correlates with the ther-nzodynatnic properlies. This lalter observalion leads to the conclusion that the mobililies of copper and gold are nearly equal. Thermodynarmic actiuilies for Cu-Ag-Au were computed from binary data and adjusted to be consistent with neasured tie lines across the two-phase region which extends into the ternary diagram at 725°C. Application of these therknodynarnic calculaliorls to the dijftcsion data terifies the Onsager reciprocal relations to the extent that experimental uncertainlies will allow. It is jound that diffusion neav the critical point of the 1200-phase field in the isothermal section at 725°C causes the composition gradients for copper and silver to approach infinity while the direct and cross coeficients become equal. From the standpoint of experimental Procedures, it has been shown that microanalysis of ternary samples may be carried out using a simple expression for the conversion of X-ray zntensily to composilion. The accuracy of' this method has been demonstrated by the analysis of standard alloys and by the reproducibility of composition determinations using differed X-ray lines. DIFFUSION in a ternary alloy system is basically and significantly different from diffusion in binary alloys. When only two elements are present there is just one independent composition variable, and the diffusion profile of composition vs distance across a sample must change monotonically through all solid solutions bracketed by the terminal alloys. If two-phase regions are stable for compositions between the terminal values, these must appear as planar interfaces in the diffused sample, with the composition profile exhibiting sharp discontinuities. Adding a third element may alter this picture completely. Because of the additional degree of freedom in the composition, we may see composition profiles which do not change monotonically through solid-solution ranges, and stable, nonplanar phase interfaces. Even though two-phase regions may lie between the terminal compositions of a diffusion sample on the constitution diagram, they will not necessarily appear in the diffusion profiles. In binary diffusion only one coefficient is needed to describe the interdiffusion process. From Fick's law this coefficient relates the diffusive flux of one component to its own composition gradient. In ternary diffusion we must allow for cross interaction between the two independent species. By a linear extension of Fick's law we introduce "cross" or "off-diagonal" diffusion coefficients which relate the flux of one component to the composition gradient of the other independent constituent. Thus, we require four interdiffusion coefficients to describe the transport process in a ternary system. In its most general form,1,2 the rate equation for mass transport by isothermal diffusion in a system of s components states that the diffusive flux of component i, Ji, is a linear summation of all driving forces, Xk, each force being multiplied by an appropriate compliance coefficient:

Jan 1, 1968

By Perry N. Halstead

The primary objective of the structural geology portion of the Kennecott Copper Corporation-U. S. Bureau of Mines Pit Slope Stability Study was to determine the attitude, geometry, and spatial distribution of geologic structures of the west wall of the Kimbley Pit for use in designing a maximum stable slope for the pit. A corollary objective was to develop and refine techniques for collecting and analyzing geologic structural data with emphasis on numerical methods that could be integrated with other aspects of pit design. Surface mapping techniques used were a modification of conventional joint set mapping and a detailed linear sampling method. For subsurface investigations both oriented core and a borehole camera were used. It was determined that the major faults in the area were favorably oriented for mining a steep stable slope. The jointing or natural fracturing, although complex, was not random in attitude and consisted of several persistent fracture sets and a number of secondary sets. The attitude of the persistent fracture sets was favorable for steeper slope angles, and these attitudes were used in the design considerations. There was no significant difference between the attitude of the fracture sets determined from the former and new slope face.

Jan 1, 1968

By C. F. Cockrell, K. K. Humphreys, H. E. Shajer

It is quite conceivable that the question of using coal or an alternate fuel in the future by large industrial consumers may be decided on the basis of possible profitable utilization of flyash. Accordingly, the Coal Research Bureau at West Virginia University has been attempting to develop new uses for this waste material. This effort has resulted in the development of a process for producing flyash-based face brick and other structural materials which generally exceed published specifications for competitive clay-based materials. These structural materials, containing up to 97% waste coal ash, can be produced from flyash derived from most coals without regard to the physical or chemical properties of the coal ash. Preliminary production cost estimates indicate that face brick can be manufactured at a rate of 9000 brick per hr or more with a resultant gross profit of over $24.00 per thousand brick. Plant payout time for a plant of this size would be only about 1.2 years assuming a market price of approximately $55.00 per thousand, a well established figure for comparable clay brick. With flyash production expected to more than double in the next 15 years, there is an increasing need for a practical and Profitable method for the utilization of large tonnages of this material (see Table I listing total flyash production of utilities in each state of the U.S.). Present industrial and governmental physical and chemical restrictions, particularly for use as a lightweight filler or pozzolanic material, have limited the use of flyash considerably. For those ashes which have met the required standards, utilization is not always assured because of lack of promotion and poor public acceptance. Less than 6% of the total flyash production in the U.S. was sold in 1962. The remainder was dumped in disposal areas at costs as high as $2.42 per ton.' It is quite conceivable that in the near future the question of using coal or an alternate fuel by large industrial consumers may be decided on the basis of the possible utilization of flyash. To help insure future coal demands, the Coal Research Bureau at

Jan 1, 1968

By Robert L. de Beauchamp

The Bureau of Mines designed a chlorination apparatus to prepare CrCl2 free of CrCl3 by the direct reaction of chlorine and chromium. Equipment and procedures were developed for routine preparation of 25-pound batch quantities of anhydrous CrCl2. Calcium, iron, nickel, and silicon impurities were in the 0.01 to 0.05 percent range, and all other impurities were below 0.01 percent.

Jan 1, 1968

By K. F. Veit, R. L. Marovelli, T. S. Chen

There is an error at the bottom of the right-hand column of p. 11, and the top left-hand column of p. 12. The first paragraph of the section EXPERIMENTAL RESULTS should read: "Thermal shock tests were conducted on circular discs of basalt, quartzite, and taconite. Average grain sizes ranged between 0.10 and 0.25 mm. Individual discs were tested only once, either in a cooling or heating shock. Several discs of the same rock type were shocked at a temperature, and the temperature level then was increased for another test by 10° to 20°F until complete fracture occurred when a specimen was quenched. The minimum temperature difference required to initiate cracking, TmiN was observed and used, together with other property values, to compute the dimensionless stress o*max . The observed ?Tmin for cooling shock in different media is tabulated in Table I. The property values E, St, a, and K are shown in Figs. 6, 8, 10, and 11, respectively. Discs of various rocks shocked in different isothermal baths gave o*max values for a wide range of Biot numbers. These o*max values were then compared with the theoretically predicted values."

Jan 1, 1968

By Y. A. Chang

Heat content values of tantalum monocarbide and tungsten monocarbide have been determined from 325" to 985°K by means of a drop-type diphenyl ether calorinzeter. Based on the values obtained in the present study and those available in the literature, best values lected and computed from 298. 15°K to temperatures close to their respectiue melting points. At temperatures higher than 1700°Kfor TaC and 27000 K for WC, values of Cp and related thermodynamic quantities were extrapolated based on the two Debye parameters: one corresponds to the metal sublattice and the other to the carbon sublattice. The contrzbutzon In the total specific heat due to lattice dilation was evaluated using the Nernst-Lindemann relalionship. The calculated values of Cp were in good agreement with the In view of their refractoriness, tantalum monocar-bide (TaC) and tungsten monocarbide (WC) are of con- experimentally determined ualues. Furthermore, the calculated values of the absolute entropies at 298.15"K were 5.03 and 4.15 in comparison to the experimental values of 5.05 and 4.15, all expressed in the unit of cal per degg-atom. Using the same Debye parameters, values of the linear coefficients of ther mal expansion for TaC and WC were calculated from absolute zero to their respective melting points using the Grüneisen equation. The calculated values are in excellent agreement with the experimental values available in the literature. Finally, the values of the bulk modulus of TaC from room tenzperature to temperature near its melting Point were also calculated using the equation recently derived by Anderson. No cot~ipari-son between the calculated and experimental values can be made due to lack of experimental data. siderable practical interest. In order to completely define these materials, it is necessary to know their thermal properties as a function of temperature. Tantalum rnonocarbide having a NaCl structure (B-1 type) melts congruently at 4250PK,1 while tungsten monoearbide with a simple hexagonal structure (Bh type) melts per itectically at 3049°K into melt and graphite.2 In the present study, the heat content values (from which

Jan 1, 1968

By Louis G. Bailey

This paper .reviews some of the key developments which have been made in the synthesis of the III-V compound semiconductors and the associated progress in obtaining high-quality material for device development. The history of the investigation of the semi-conducting properties of these materials has indicated the need for synthesis techniques capable of producing rnaterial of high chemical purity and physical perfection. The progress in the area is described, and the present status of quality of several compounds is discussed. In addition, recent studies of localized epitaxial deposition of GaAs into semi-insulating GaAs substvates for the development of arrays of devices ave reviewed. In conclusion, some important areas for future research and developrnent are discussed. INTEREST in III-V compounds as semiconductors dates back to the period 1950 to 1952, when Goryunova, Ioffe, and, in particular, Welker determined that 111-V compounds were structural analogs of the group IV elements carbon, silicon, germanium, and tin and that they were also semiconductors.'-3 Since that time interest in the properties of these materials has ex- panded rapidly. During 1953 and 1954, the first reports of rectification, photoconductivity, and optical absorption in these materials came from several investigators. In 1955, recombination radiation from samples of GaAs, InP, and GaSb was reported by ~raunstein .' This radiation was produced by injected minority carriers from point and large-area contacts. From 1955 to 1960, efforts increased on preparing large, pure specimens of these materials for the study of basic physical phenomena and device exploration. By 1960, materials preparation and the understanding of the physics of some of the materials had advanced to where some very interesting and potentially useful phenomena came under study. The first observation of electric-field-induced changes in the optical absorption edge (Franz -Keldysh effect) of 111-V compounds was reported for GaAs by Moss in 1961.' In 1962 the study of recombination radiation advanced rapidly. Keyes and Quist6 and Pankove and Berkey-heiser7 reported high radiative efficientices obtained from recombination, and then the first injection electroluminescent lasers in GaAs were announced by Hall ef al.,' by Nathan et al.,O and by Quist ef al.,'O and in Ga(As, P) alloys by Holonyak and Bevacqua." Hil-sum has reviewed the early work leading to the first observations of laser action in GaAS." Recombination processes in GaAs, Gap, and the other 111-V com-

Jan 1, 1968

By S. F. Bartram, G. A. Chase

COMPUTED X-ray diffraction powder patterns of YNi3, YNi5, and CbBe3 are available on request from either author.

Jan 1, 1968

By J. W. Town

The Bureau of Mines determined the physical properties of southeastern Idaho phosphate fines as related to mineral composition, particle size and liberation, beneficiation, and potential use in structural clay products. The phosphate mineral in the fines was identified by unit-cell dimensions to be essentially the same as the phosphate mineral fluorapatite identified in the coarser sizes. The mean particle frequency as determined by electron-microscopy and pipet-analyses was between 0.5 and 0.2 micron. Liberation of the phosphate could not be definitely determined because of the weathered condition of the minerals and their extreme fineness. No beneficiation scheme was developed that would concentrate the phosphate into a marketable fertilizer product; therefore, forming and firing tests were made to demonstrate that the fines can be used to manufacture structural clay products.

Jan 1, 1968

By Paul Cherin

ThE angle of intersection, 0, between h2k2l2 planes and the (loo), (110), (101), and (111) planes was calculated on an IBM 7040 using the well-known relationship

Jan 1, 1968