summaryrefslogtreecommitdiff
path: root/1332/CH19/EX19.3/19_3.sce
blob: 7d95b8ddabfca6e9bf90900400260ee5eeb6dd37 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
//Example 19.3
//Simple Explicit Method
//Page no. 658
clc;clear;close;

c=-2;dt=0.07;dx=0.2;
r=abs(c)*dt/dx;
printf('\n  x\ti\t|\tj -->\t')
for i=0:6
    printf('  %i\t',i)
end
printf('\n  |\t|\t|\tt -->\t')
for i=0:6
    printf('%.3f\t',i*dt)
end
printf('\n---------------------------------------------------------------------------------------')
for j=1:6
    printf('\n %.1f\t%i\t|\t\t',(j-1)*dx,j-1)
    for i=1:7
        if i==1 then
            u(j,i)=0;
        elseif j==1 then
            u(j,i)=1
        else
            u(j,i)=(1-r)*u(j,i-1)+r*u(j-1,i-1)
        end
        printf('%.3f\t',u(j,i))
    end
end