SUBROUTINE SPLICO (X,Y,M,C) 



DIMENSION X(1),Y(1),C(4,53),D(53).P(53),E<53).A(53,3),B(53),Z(53) 



MM=M-1 



DO 2 K=1,MM 



D(K 



P(K 



E(K 



00 



B(K 



A(l 



A(l 



A(2 



A(2 



A(2 



B(2 



=X(K+1)-X(K) 

 =D(K)/6. 



=(Y(K+1)-Y(K))/D(K) 

 3 K=2,MM 

 =E(K)-E(K-1) 

 2)=-l.-D(l)/D(2) 

 3)=D(1)/D(2) 

 3)=P(2)-P(1)*A(1,3) 

 2)=2.*(P(1)+P(2))-P(1)*A(1,2) 

 3)=A(2,3)/A(2,2) 

 =B(2)/A(2,2) 



DO 4 K=3,MM 



A(K,2)=2.*(P(K-1)+P(K))-P(K-1)*A(K-1,3) 



B(K)=B(K)-P(K-1)*B(K-1) 



A(K,3)=P(K)/A(K,2) 



B(K)=B(K)/A(K,2) 



Q=D(M-2)/D(M-l) 



A(M,l)=l.+Q+A(M-2,3) 



A(M,2)=-Q-A(M,1)*A(M-1,3) 



B(M)=B(M-2)-A(M,l)*B(M-I) 



Z(M)=B(M)/A(M,2) 



MN=M-2 



DO 6 I=1,MN 



K=M-I 



Z(K)=B(K)-A(K,3)*Z(K+1) 



Z(1)=-A(1,2)*Z(2)-A(1,3)*Z(3) 



DO 7 K=1,MM 



Q=1./(6.*D(K)) 



C(1,K)=Z(K)*Q 



C(2,K)=Z(K+1)*Q 



C(3,K)=Y(K)/D(K)-Z(K)*P(K) 



C(4,K)=Y(K+1)/D(K)-Z(K+1)*P(K) 



RETURN 

 END 



168 



