example9_6.py

## example9_6
from numarray import ones
from inversePower5 import *

def Bv(v):          # Compute {z} = [B]{v}
      n = len(v)
      z = zeros((n),type=Float64)
      z[0] = 2.0*v[0] - v[1]
      for i in range(1,n-1):
          z[i] = -v[i-1] + 2.0*v[i] - v[i+1]
      z[n-1] = -v[n-2] + 2.0*v[n-1]
      return z

n = 100                 # Number of interior nodes
d = ones((n))*6.0       # Specify diagonals of [A] = [f\e\d\e\f]
d[0] = 5.0                             
d[n-1] = 7.0
e = ones((n-1))*(-4.0)
f = ones((n-2))*1.0
lam,x = inversePower5(Bv,d,e,f)
print "PL^2/EI =",lam*(n+1)**2
raw_input("\nPress return to exit")