Advertisement
Guest User

Untitled

a guest
Dec 6th, 2016
55
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.08 KB | None | 0 0
  1. double[] pastX = new double[n];
  2. double l1 = 0;
  3. double l2 = 0;
  4. double REL = 0;
  5. double DIFF = 0;
  6. do{
  7. for (int i =0; i < n; i++){
  8. pastX[i] = currentX[i];
  9. }
  10. for (int i = 0; i < n; i++){
  11. currentX[i]=0;
  12. for(int j = 0; j < n; j++){
  13. currentX[i] += A1[i][j]*pastX[j];
  14. }
  15. }
  16. l1=0;
  17. l2=0;
  18. for (int i = 0; i < n; i++){
  19. l1+=currentX[i];
  20. l2+=pastX[i];
  21. }
  22. l1=sqrt(l1);
  23. l2=sqrt(l2);
  24. for (int i = 0; i < n; i++){
  25. currentX[i]/=l1;
  26. pastX[i]/=l2;
  27. }
  28. REL = l1/l2;
  29. DIFF = abs(currentX[0]-pastX[0]);
  30. for (int i = 1; i < n; i++){
  31. if (abs(currentX[i]-pastX[i]) > DIFF){
  32. DIFF = abs(currentX[i]-pastX[i]);
  33. }
  34. }
  35. } while (DIFF < eps);
  36. System.out.println(1/REL);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement