Pedagogs international research journal

“PEDAGOGS” 
 international research journal ISSN: 
2181-4027
_SJIF: 
4.995
 
www.pedagoglar.uz
 
Volume-8, Issue-4, April - 2022
 
253 
2. M.file da f1 ni saqlab olinadi. 
function
y=f1(x,t) 
y=x-t; 
end 
3. M.file da u0 ni saqlab olinadi. 
function
y=u0(x) 
y=x-x^2; 
end 
4. M.file da g1 ni saqlab olinadi. 
function
z=g1(t) 
z=t; 
end
5. M.file da g2 ni saqlab olinadi. 
function
y=g2(t) 
y=t; 
end
2. M.file da pdt_Oshkormas_teskari_matrissa
ni saqlab olinadi. Shu faylni saqlagandan keyin 
dasturni ishga tushiramiz. Iloji bo`lsa, bitta papkaga hammasini saqlashimiz lozim. Dasturni ishga 
tushirishda aniq harakat qilish kerak.
function
U=pdt_Oshkormas_teskari_matrissa(u0,f,c1,c2,a,b,n,m) 
h=a/(n-1); tau=b/(m-1); r=tau/h^2; U=zeros(n,m); A=zeros(n,n); U(1,1:m)=c1; U(n,1:m)=c2; 
for
j=1:m 
for
i=1:n 
f1(i,j)=feval(f,h*(i-1),tau*(j-1)); 
end
end
for
i=1:n
U(i,1)=feval(
'u0'
,h*(i-1)); 
end
A(1,1)=1;A(n,n)=1; 
for
i=2:n-1 
A(i,i-1)=-r;A(i,i)=1+2*r;A(i,i+1)=-r 
end

---

“PEDAGOGS” 
 international research journal ISSN: 
2181-4027
_SJIF: 
4.995
 
www.pedagoglar.uz
 
Volume-8, Issue-4, April - 2022
 
254 
%AU=D ni yaratish va yechish
D(1)=feval(
'g1'
,0);D(n)=feval(
'g2'
,b); 
for
j=2:m 
for
i=2:n-1 
D(i)=U(i-1,j)+tau*f1(i,j); 
end
V=inv(A);X=V*D';U(1:n,j)=X'; 
end
U=U'; 
% Command Window oynasiga quyidagi kodni yozish lozim
% [U]=pdt_Oshkormas_teskari_matrissa('u0','f',0,0,1,0.05,10,10)
% Grafigi mesh(U)
3. Natija chiqarish bo`limi. 
Command Window oynasiga quyidagilarni yozish lozim
>> [U]=pdt_Oshkormas_teskari_matrissa('u0','f',0,0,1,0.05,10,10) 
A = 
1.0000 0 0 0 0 0 0 0 0 0 
-0.4500 1.9000 -0.4500 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 1.0000 
A = 
1.0000 0 0 0 0 0 0 0 0 0 
-0.4500 1.9000 -0.4500 0 0 0 0 0 0 0 
0 -0.4500 1.9000 -0.4500 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 0 

---