# test_deriv.py  Q&D test of deriv.py, see nderiv.py for full output

from deriv import *  # deriv, rderiv


print "test_deriv.py running"
a,b = deriv(1,3,1)
print "a,b = deriv(1,3,1)"
print "a=",a
print "b=",b

npoints = 5
h = 1.0/(npoints-1);
for order in range(1,npoints):
    for term in range (0,order):
        c = rderiv(order, npoints, term, h)
        print "order=",order,", term=",term,", h=",h,", n=",npoints
        print c

# end test_deriv.py