|
DataMuseum.dkPresents historical artifacts from the history of: DKUUG/EUUG Conference tapes |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about DKUUG/EUUG Conference tapes Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - downloadIndex: S T
Length: 1004 (0x3ec)
Types: TextFile
Names: »S00099.f«
└─⟦db229ac7e⟧ Bits:30007240 EUUGD20: SSBA 1.2 / AFW Benchmarks
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21E/doduc/S00099.f«
└─⟦this⟧ »EUUGD20/AFUU-ssba1.21/ssba1.21F/doduc/S00099.f«
SUBROUTINE S00099(NP,X,T,C)
IMPLICITDOUBLEPRECISION(A-H,O-Z)
DOUBLEPRECISION X(5),T(5),Y(5),Z(5),C(5)
********************************************************************************
* DO 1000 I=1,5 *
********************************************************************************
DO 1000 I=1,NP
1000 C(I)=0.
IF(NP.GT.3)GOTO 4000
DO 3100 I=2,3
3100 Y(I)=(T(I)-T(1))/(X(I)-X(1))
C(1)=(Y(3)-Y(2))/(X(3)-X(2))
C(2)=Y(2)-(X(1)+X(2))*C(1)
C(3)=T(1)-X(1)*C(2)-X(1)*X(1)*C(1)
RETURN
4000 IF(NP.GT.4)GOTO 5000
DO 4100 I=2,4
4100 Y(I)=(T(I)-T(1))/(X(I)-X(1))
DO 4200 I=3,4
4200 Z(I)=(Y(I)-Y(2))/(X(I)-X(2))
C(1)=(Z(4)-Z(3))/(X(4)-X(3))
C(2)=Z(3)-(X(1)+X(2)+X(3))*C(1)
C(3)=Y(2)-C(1)*(X(2)*X(2)+X(1)*X(2)+X(1)*X(1))-C(2)*(X(1)+X(2))
C(4)=T(1)-C(1)*X(1)*X(1)*X(1)-C(2)*X(1)*X(1)-C(3)*X(1)
RETURN
5000 IF(NP.GT.5)RETURN
RETURN
END