Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- double[] pastX = new double[n];
- double l1 = 0;
- double l2 = 0;
- double REL = 0;
- double DIFF = 0;
- do{
- for (int i =0; i < n; i++){
- pastX[i] = currentX[i];
- }
- for (int i = 0; i < n; i++){
- currentX[i]=0;
- for(int j = 0; j < n; j++){
- currentX[i] += A1[i][j]*pastX[j];
- }
- }
- l1=0;
- l2=0;
- for (int i = 0; i < n; i++){
- l1+=currentX[i];
- l2+=pastX[i];
- }
- l1=sqrt(l1);
- l2=sqrt(l2);
- for (int i = 0; i < n; i++){
- currentX[i]/=l1;
- pastX[i]/=l2;
- }
- REL = l1/l2;
- DIFF = abs(currentX[0]-pastX[0]);
- for (int i = 1; i < n; i++){
- if (abs(currentX[i]-pastX[i]) > DIFF){
- DIFF = abs(currentX[i]-pastX[i]);
- }
- }
- } while (DIFF < eps);
- System.out.println(1/REL);
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement