Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- a = []
- b = []
- c = []
- d = []
- x=[]
- alpha = []
- beta = []
- gamma = []
- element = 0.0
- n=input()
- for i in range(0,n):
- element=input()
- b.append(element)
- for i in range(0,n-1):
- element=input()
- a.append(element)
- for i in range (0,n-1):
- element=input()
- c.append(element)
- for i in range (0,n):
- element=input()
- d.append(element)
- gamma.append(b[0])
- beta.append(d[0]/(gamma[0]*1.0))
- alpha.append(-c[0]/(b[0]*1.0))
- for i in range(1,n-1):
- gamma.append(b[i]+a[i]*alpha[i-1])
- alpha.append(-c[i]/(gamma[i]*1.0))
- beta.append((d[i]-a[i]*beta[i-1])/(gamma[i]*1.0))
- x = range(0,n)
- x[n-1]=(beta[n-2])
- for i in range(n-1,1,-1):
- x[i-1]=(alpha[i]*x[i]+beta[i])
- print x
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement