def check_spint_spalde(self):
per = [0, 0, 0]
N = [20, 20, 50]
ia = [0, 0, 0.2*pi]
ib = [0, 0, pi]
a,b = 0,2*pi
dx = 0.2*pi
k = range(1,6)
for i in range(len(per)):
x=a+(b-a)*arange(N[i]+1,dtype=float)/float(N[i])
v=f1(x)
for j in range(len(k)):
uspl = UnivariateSpline(x,v,k=k[j],s=0)
ir = uspl.integral(ia[i],ib[i])
dr = uspl.derivatives(dx)
assert_almost_equal(ir, f1(ib[i],-1)-f1(ia[i],-1), decimal=2)
d=0
for ddr in dr:
if d<k[j]-1:
assert_almost_equal(1, ddr/f1(dx,d), decimal=2)
d=d+1
