Roots of a quadratic function: Difference between revisions
Content added Content deleted
m (Added the Sidef language) |
PatGarrett (talk | contribs) |
||
Line 565: | Line 565: | ||
=={{header|Fortran}}== |
=={{header|Fortran}}== |
||
=={{{header|Fortran 90}}}== |
|||
{{works with|Fortran|90 and later}} |
{{works with|Fortran|90 and later}} |
||
<lang fortran>PROGRAM QUADRATIC |
<lang fortran>PROGRAM QUADRATIC |
||
Line 616: | Line 617: | ||
The roots are real: |
The roots are real: |
||
Root1 = 0.999999999999000E+06 Root2 = 0.100000000000100E-05 |
Root1 = 0.999999999999000E+06 Root2 = 0.100000000000100E-05 |
||
=={{{header|Fortran I}}}== |
|||
Source code written in FORTRAN I (october 1956) for the IBM 704. |
|||
<lang fortran> |
|||
COMPUTE ROOTS OF A QUADRATIC FUNCTION - 1956 |
|||
READ 100,A,B,C |
|||
100 FORMAT(3F8.3) |
|||
PRINT 100,A,B,C |
|||
DISC=B**2-4.*A*C |
|||
IF(DISC),1,2,3 |
|||
1 XR=-B/(2.*A) |
|||
XI=SQRT(-DISC)/(2.*A) |
|||
XJ=-XI |
|||
PRINT 311 |
|||
PRINT 312,XR,XI,XR,XJ |
|||
311 FORMAT(13HCOMPLEX ROOTS) |
|||
312 FORMAT(4HX1=(,2E12.4,6H),X2=(,2E12.4,1H)) |
|||
GO TO 999 |
|||
2 X1=-B/(2.*A) |
|||
X2=X1 |
|||
PRINT 321 |
|||
PRINT 332,X1,X2 |
|||
321 FORMAT(16HEQUAL REAL ROOTS) |
|||
GO TO 999 |
|||
3 X1= (-B+SQRT(DISC)) / (2.*A) |
|||
X2= (-B-SQRT(DISC)) / (2.*A) |
|||
PRINT 331 |
|||
PRINT 332,X1,X2 |
|||
331 FORMAT(10HREAL ROOTS) |
|||
332 FORMAT(3HX1=,E12.5,4H,X2=,E12.5) |
|||
999 STOP |
|||
</lang> |
|||
=={{header|GAP}}== |
=={{header|GAP}}== |