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, 1 * ones(1, k), k, k);
  5. Elow = sparse(2 : k, 1 : k-1, -2 *ones(1, k-1), k, k);
  6. Elow2 = sparse(3 : k, 1 : k - 2, ones(1, k - 2), k, k);
  7. S = Elow + D + Elow2;
  8. S(1,1) = 1;
  9. S(1,2) = 0;
  10. S(2,1) = 0; S(2,2) = 1; S(2,3) = 0;
  11. b = h^2 * ones(1,k)'; b(1) = 0; b(2) = 0;
  12. u = S \ b
  13. plot(x, u)
Success #stdin #stdout 0.31s 104960KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
u =

   0.00000
   0.00000
   0.01235
   0.03704
   0.07407
   0.12346
   0.18519
   0.25926
   0.34568
   0.44444
https://ideone.com/L5txJR
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!