fork(1) download
copy 
  1. yy = [-0.729045991406, 0.670017179989, 0.937735542248, -0.557931914035, -0.383665898986, 0.0944411927357, -0.623236047906, 0.985762603836]';
  2. AA =  [-1962    389      0    921    539    113      1      0;...
  3.          389  -1936      0    866    195    486      0      0;...
  4.            0      0   -632    332    300      0      0      1;...
  5.          921    866    332  -3940    950    871      0      0;...
  6.          539    195    300    950  -2951    967      0      0;...
  7.          113    486      0    871    967  -2437      0      0;...
  8.           -1      0      0      0      0      0      0      0;...
  9.            0      0     -1      0      0      0      0      0];
  10.  
  11. x = zeros(1, 8);
  12. x(1) = yy(7) / AA(7, 1);
  13. x(3) = yy(8) / AA(8, 3);
  14. yy = yy - AA(:, 1) * x(1);
  15. yy = yy - AA(:, 3) * x(3);
  16.  
  17. AA(6,:) += AA(5,:)*0.327685530329;
  18. yy(6,:) += yy(5,:)*0.327685530329;
  19. AA(5,:) += AA(6,:)*0.456104517255;
  20. yy(5,:) += yy(6,:)*0.456104517255;
  21.  
  22. A = AA([2, 4, 5, 6],[2, 4, 5, 6])
  23. y = yy([2, 4, 5, 6])
Success #stdin #stdout 0.08s 47024KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
A =

  -1.9360e+03   8.6600e+02   1.9500e+02   4.8600e+02
   8.6600e+02  -3.9400e+03   9.5000e+02   8.7100e+02
   4.4581e+02   1.4893e+03  -2.9510e+03   1.1910e-09
   5.4990e+02   1.1823e+03  -8.7903e-10  -2.1201e+03

y =

  -241.769
  -247.285
   -78.722
   -83.628
https://ideone.com/iP96PI
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!