# test_laphi.py  test  laphi.py  function phi, phip, ...
from laphi import *    # phi(x,j,n,xg), shape function
                       # phip, first derivative
                       # phipp, second derivative
                       # phippp, third derivative
                       # phipppp, fourth derivative

n = 4
h = 0.25
print "n = ", n
xg = [h*i for i in range(n+1)]
print "xg = ", xg

print "phi, shape function"
for j in range(0,n+1):
  print " "
  print "j = ",j,", at x = ", j*h
  for k in range(0,2*n+1):
    x = (h/2.0)*k
    y = phi(x,j,n,xg)
    print "x= ",x,", y = phi(x,j,n,xg) = ", y

print "phip, first derivative"
for j in range(0,n+1):
  print " "
  print "j = ",j,", at x = ", j*h
  for k in range(0,2*n+1):
    x = (h/2.0)*k
    yp = phip(x,j,n,xg)
    print "x= ",x,", yp = phip(x,j,n,xg) = ", yp

print "phipp, second derivative"
for j in range(0,n+1):
  print " "
  print "j = ",j,", at x = ", j*h
  for k in range(0,2*n+1):
    x = (h/2.0)*k
    ypp = phipp(x,j,n,xg)
    print "x= ",x,", ypp = phipp(x,j,n,xg) = ", ypp

print "phippp, third derivative"
for j in range(0,n+1):
  print " "
  print "j = ",j,", at x = ", j*h
  for k in range(0,2*n+1):
    x = (h/2.0)*k
    yppp = phippp(x,j,n,xg)
    print "x= ",x,", yppp = phippp(x,j,n,xg) = ", yppp

print "phipppp, fourth derivative"
for j in range(0,n+1):
  print " "
  print "j = ",j,", at x = ", j*h
  for k in range(0,2*n+1):
    x = (h/2.0)*k
    ypppp = phipppp(x,j,n,xg)
    print "x= ",x,", ypppp = phipppp(x,j,n,xg) = ", ypppp

# end test_laphi.py