fork(1) download
copy 
  1. A =  [-578      0    312      0    266      0      0;...
  2.          0   -323      0      0    323      1      0;...
  3.        312      0   -821      0    509      0      0;...
  4.          0      0      0   -726    726      0      1;...
  5.        266    323    509    726  -1824      0      0;...
  6.          0     -1      0      0      0      0      0;...
  7.          0      0      0     -1      0      0      0];
  8. y = [-0.729046, 0.670017, 0.937736, -0.557932, -0.383666, 0.0944412, -0.623236];
  9.  
  10. x = zeros(1, 7);
  11. x(1, 2) = y(1, 6) / A(6, 2);
  12. x(1, 4) = y(1, 7) / A(7, 4);
  13.  
  14. y = y' - A(:, 2) * x(2);
  15. y = y - A(:, 4) * x(4);
  16. AA = A([1, 3, 5], [1, 3, 5]);
  17. yy = y([1, 3, 5]);
  18.  
  19. AA(2,:) += AA(1,:) * 0.539792;
  20. AA(1,:) += AA(2,:) * 0.478099;
  21.  
  22. yy(2) += yy(1) * 0.539792;
  23. yy(1) += yy(2) * 0.478099;
  24.  
  25. coeff1 = (- AA(1, 3)) / AA(1, 1) * AA(3, 1)
  26. coeff2 = (- AA(2, 3)) / AA(2, 2) * AA(3, 2)
  27.  
  28. AA
  29.  
  30. disp(AA^(-1) * yy);
Success #stdin #stdout 0.11s 47652KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
coeff1 =  266.00
coeff2 =  509.00
AA =

   -577.99989291     -0.00018619    578.00007910
      0.00022400   -652.58489600    652.58467200
    266.00000000    509.00000000  -1824.00000000

   0.40323
   0.40159
   0.40242
https://ideone.com/Vby2X0
language:
Octave (octave 4.4.1)
created:
7 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!