|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 768 (0x300)
Types: TextFile
Names: »tspln3dif2«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦09b4e9619⟧ »thcømat«
└─⟦this⟧
message spln3dif2
spln3dif2=algol message.no
cubic spline 2.derivative
external
real procedure spln3dif2(x,X,Y,M,n);
value n,x; real x; integer n; array X,Y,M;
begin
real xm,xp;
own integer i;
if x<X(1) or x>X(n) then alarm(<:<10>***spln3dif2 illegal x = :>,
string exactlay(x,i,x),x);
if i>=n then i:=0;
for xp:=X(i+1)-x while xp<=0 and i<n-1 do i:=i+1;
for xm:=x-X(i) while xm<0 do
begin
xp:=-xm; i:=i-1
end;
spln3dif2:=(M(i+1)*xm+M(i)*xp)/(xm+xp)
end spln3dif2
; end
▶EOF◀