fork download
copy 
  1. k = 10;
  2. x = linspace(0, 1, k);
  3. h = (x(end) - x(1)) / (k - 1);
  4. D = sparse(1 : k, 1 : k, -2 * ones(1, k), k, k);
  5. E = sparse(2 : k, 1 : k - 1, ones(1, k-1), k, k);
  6. S = E + D + E';
  7. S(1,1) = 1;
  8. S(1,2) = 0;
  9. S(2,1) = 0; S(2,2) = 1; S(2,3) = 0;
  10. Sf = full(S);
  11. b = h^2 * ones(1,k)'; b(1) = 0; b(2) = 0;
  12. u = S \ b;
  13.  
  14. u'
  15.  
Success #stdin #stdout 0.21s 103360KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
ans =

 Columns 1 through 8:

   0.00000   0.00000  -0.04938  -0.08642  -0.11111  -0.12346  -0.12346  -0.11111

 Columns 9 and 10:

  -0.08642  -0.04938
https://ideone.com/MOXNxB
Finite Difference method Octave
language:
Octave (octave 4.4.1)
created:
7 years 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!