blob: a3d627066ef3d7c7bd5e4bec63443d7b16d417de (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
//Example No. 15_05
//Initial Value Problems
//Pg No. 494
clear ; close ; clc ;
h = 1 ;
k = 2 ;
tau = h^2/(2*k)
for i = 2:4
f(1,i) = 50*( 4 - (i-1) )
end
f(1:7,1) = 0 ;
f(1:7,5) = 0 ;
for j = 1:6
for i = 2:4
f(j+1,i) = ( f(j,i-1) + f(j,i+1) )/2
end
end
disp(f,'The final results are ')
|
