Mathematica 2.0 for MS-DOS 386/7 (June 21, 1991)
Copyright 1988-91 Wolfram Research, Inc.

In[1]:= <<na4.m
Roman Maurer, uporabna matematika, 29-Aug-1996
4. domača naloga iz Numerične analize:
Naloga 1. c; rešitev na [0,1]x[0,1] PDE  u_t = u_xx,
               u(x,0) = sin(pi*x),
               u(0,t) = u(1,t) = 0.
- z eksplicitno metodo (dx=0.1, lambda=1/4),
- z implicitno metodo (dx=0.1, lambda=1/4),
- s Crank-Nicholsonovo metodo (dx=0.1, lambda=1).
 
metoda                   čas      rel. napaka
eksplicitna:        {1.6 Second, -0.0400328}
implicitna:         {1.75 Second, 0.219892}
Crank-Nicholson:    {1.7 Second, 0.0758169}

In[2]:= u1

Out[2]= {0.0000159833, 0.0000304021, 0.0000418449, 0.0000491917, 
 
>    0.0000517232, 0.0000491917, 0.0000418449, 0.0000304021, 0.0000159833}

In[3]:= ExplicitAt1

Out[3]= {0.0000153435, 0.000029185, 0.0000401698, 0.0000472224, 0.0000496526, 
 
>    0.0000472224, 0.0000401698, 0.000029185, 0.0000153435}

In[4]:= ImplicitAt1

Out[4]= {0.0000194979, 0.0000370873, 0.0000510463, 0.0000600085, 
 
>    0.0000630967, 0.0000600085, 0.0000510463, 0.0000370873, 0.0000194979}

In[5]:= CrankNicholsonAt1

Out[5]= {0.0000171952, 0.0000327071, 0.0000450175, 0.0000529212, 
 
>    0.0000556447, 0.0000529212, 0.0000450175, 0.0000327071, 0.0000171952}