veto14

ALC.Subst1

Oct 10th, 2018
65
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.85 KB | None | 0 0
  1. A = []
  2. with open('somefile.txt', 'r') as f:
  3. for line in f:
  4. A.append(list(map(int,line.split())))
  5. tamanho = len(A[0]) #adquire se o tamanho da matriz
  6. truetam = tamanho-1
  7. satisfaz = 0
  8. x = [0] * truetam;
  9. mat = [0] * truetam;
  10. for i in range(truetam):
  11. mat[i] = [0]*truetam # separamos a matriz Ax do vetor b
  12. for i in range(truetam):
  13. for j in range(truetam):
  14. mat[i][j] = A[i][j]
  15. Vx = [0] * truetam
  16. for i in range(truetam):
  17. Vx[i] = A[i][3] #aqui se cria o vetor b
  18. for i in range(truetam):
  19. existencia = truetam - (i+1)
  20. cont = 0
  21. for j in range(truetam):
  22. if(j>i):
  23. if(mat[i][j] == 0):
  24. cont = cont+1
  25. if(cont == existencia):
  26. satisfaz = satisfaz + 1
  27. if(satisfaz == truetam):
  28. for k in range(truetam):
  29. x[0] = (Vx[0] / mat[0][0])
  30. if(k==0):
  31. k=k+1
  32. cont=1
  33. results = [0] * k
  34. for m in range(k):
  35. if(cont<=(truetam-1)):
  36. if(k==1):
  37. x[cont]= (Vx[cont] - (mat[k][m]*x[m])/mat[k][k])
  38. if(m == k-1):
  39. cont +=cont
  40. k+=1
  41. if(k>1):
  42. while(m!=(k-1)):
  43. soma = Vx[cont] - (mat[k][m]*x[m])
  44. m+=1
  45. if(m==(k-1)):
  46. soma -= (mat[k][m]*x[m])
  47. soma = soma/mat[k][k]
  48. x[cont] = soma
  49. if(m == k-1):
  50. cont +=cont
  51. print("As raízes do seu sistema são:");
  52. for fim in range(truetam):
  53. print(x[fim])
  54. else:
  55. print("Não é hoje que você resolve seu sistema, senhor Tobias, tente no próximo Show do Milhão!")
Add Comment
Please, Sign In to add comment