4 курс_1 / Мат.мод 313 / Лаб6

function u=func1(m,n,tau)

u=zeros(m,n);

for i=1:m

u(i,1)=1;

end

for j=1:n

u(1,j)=1+tau*(j-1);

end

for j=2:n

for i=2:m-1

u(i,j)=0.5*(u(i+1,j)+u(i-1,j));

end

u(m,j)=1;

end

u=u';

>> u=func1(5,7,0.16)

u =

1.0000 1.0000 1.0000 1.0000 1.0000

1.1600 0.5800 0.2900 0.1450 1.0000

1.3200 0.6600 0.3300 0.1650 1.0000

1.4800 0.7400 0.3700 0.1850 1.0000

1.6400 0.8200 0.4100 0.2050 1.0000

1.8000 0.9000 0.4500 0.2250 1.0000

1.9600 0.9800 0.4900 0.2450 1.0000

>> surf(u)

>> u=func1(11,9,0.16)

u =

1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000

1.1600 0.5800 0.2900 0.1450 0.0725 0.0362 0.0181 0.0091 0.0045 0.0023 1.0000

1.3200 0.6600 0.3300 0.1650 0.0825 0.0413 0.0206 0.0103 0.0052 0.0026 1.0000

1.4800 0.7400 0.3700 0.1850 0.0925 0.0462 0.0231 0.0116 0.0058 0.0029 1.0000

1.6400 0.8200 0.4100 0.2050 0.1025 0.0513 0.0256 0.0128 0.0064 0.0032 1.0000

1.8000 0.9000 0.4500 0.2250 0.1125 0.0563 0.0281 0.0141 0.0070 0.0035 1.0000

1.9600 0.9800 0.4900 0.2450 0.1225 0.0612 0.0306 0.0153 0.0077 0.0038 1.0000

2.1200 1.0600 0.5300 0.2650 0.1325 0.0663 0.0331 0.0166 0.0083 0.0041 1.0000

2.2800 1.1400 0.5700 0.2850 0.1425 0.0713 0.0356 0.0178 0.0089 0.0045 1.0000

>> surf(u)

>>