|
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: »tspln3dif«
└─⟦00964e8f7⟧ Bits:30007478 RC8000 Dump tape fra HCØ.
└─⟦b2ec5d50f⟧
└─⟦09b4e9619⟧ »thcømat«
└─⟦this⟧
message spln3dif
spln3dif=algol message.no
cubic spline 1. derivative
external
real procedure spln3dif(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>***spln3dif 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;
spln3dif:=(Y(i+1)-Y(i)-xp*xm*(M(i+1)-M(i))/2)/(xm+xp)
+(M(i+1)*(xm+xm-xp)+M(i)*(xm-xp-xp))/6
end spln3dif
; end
▶EOF◀