Advertisement
Guest User

Untitled

a guest
Sep 26th, 2017
70
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.86 KB | None | 0 0
  1. # n = 0
  2. a = []
  3. b = []
  4. c = []
  5. d = []
  6. x=[]
  7. alpha = []
  8. beta = []
  9. gamma = []
  10. element = 0.0
  11. n=input()
  12. for i in range(0,n):
  13. element=input()
  14. b.append(element)
  15. for i in range(0,n-1):
  16. element=input()
  17. a.append(element)
  18. for i in range (0,n-1):
  19. element=input()
  20. c.append(element)
  21. for i in range (0,n):
  22. element=input()
  23. d.append(element)
  24. gamma.append(b[0])
  25. beta.append(d[0]/(gamma[0]*1.0))
  26. alpha.append(-c[0]/(b[0]*1.0))
  27. for i in range(1,n-1):
  28. gamma.append(b[i]+a[i-1]*alpha[i-1])
  29. alpha.append(-c[i]/(gamma[i]*1.0))
  30. beta.append((d[i]-a[i-1]*beta[i-1])/(gamma[i]*1.0))
  31. gamma.append (b[n-1]+a[n-2]*alpha[n-2])
  32. beta.append ((d[n-1]-a[n-2]*beta[n-2])/(gamma[n-1]*1.0))
  33.  
  34. x = range(0,n)
  35. x[n-1]=(beta[n-1])
  36. for i in range(n-1,1,-1):
  37. x[i-1]=(alpha[i-1]*x[i]+beta[i-1])
  38. print x
  39. x[0]=alpha[0]*x[1]+beta[0]
  40. print x
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement