spline(tab,p)
Func

	Local i,n,mat,h,l,mi,d,m,s,a,b

	colDim(tab)-1 -> n
	diag(newMat(1,n).+2) -> mat
	newList(n-1) -> h

	For i,1,n-1,1
		tab[1,i+1]-tab[1,i] -> h[i]
	EndFor
		
	1->mat[1,2]
	1->mat[n,n-1]
	
	For i,2,n-1,1
		h[i]/(h[i-1]+h[i]) -> l
		1-l -> mi
		mi -> mat[i,i-1]
		l -> mat[i,i+1]
	EndFor

	newList(n) -> d
	6/(h[1])*((tab[2,2]-tab[2,1])/(h[1])-tab[1,n+1]) -> d[1]
	6/(h[n-1)*(tab[2,n+1]-(tab[2,n]-tab[2,n-1])/(h[n-1])) -> d[n]

	For i,2,n-1,1
		6/(h[i-1]+h[i])*((tab[2,i+1]-tab[2,i])/(h[i])-(tab[2,i]-tab[2,i-1])/(h[i-1])) -> d[i]
	EndFor

	jordan(augmen(mat,list>mat(d,1))) -> m

	newList(n-1) -> s

	For i,1,n-1,1
		(tab[2,i+1]-tab[2,i])/(h[i]) - h[i]/6*(m[i+1,1]-m[i,1]) -> a
		tab[2,i] - h[i]^2/6*m[i,1] -> b
		m[i,1]*(tab[1,i+1]-x)^3/(6*h[i]) + m[i+1,1]*(x-tab[1,i])^3/(6*h[i]) + a*(x-tab[1,i]) + b -> s[i]
	EndFor

	Return s[p]

EndFunc