Here is a plot of the numerical solution to the heat equation (in time and space) on the intervall [0,1] with known border value and initial value conditions using the Crank-Nicolson method.
If you should be interested; the calculations were done in Octave with the following files:
- main.m: main script
- f.m: initial values for t = 0
- g0.m: border values for x = 0
- g1.m: border values for x = 1
- G.m: function to calculate the vector G used in main.m.
As you might notice the scripts are general so that it is possible to change the border and initial values. (I hope.)
Also, some credit should go to my friend Jens who helped me make the scripts, but if there are any mistakes they are obviously his.
Comments (2)
Ble 'linket' hit av reklame på datasal i Realfagbygget. Flere fine bilder, god løsning av varmeledningslikninga, mange gode tegneserier.
hehe. gøy å høre at markedsføringen min bærer frukter, og glad bildene falt i smak :)