Samarqand davlat universiteti giperbolik tipdagi tenglamali chegaraviy masalalarni sonli

 
2-ilova. O’zgaruvchan koeffisiyentli ko'chirish tenglamasini har xil ayirmali 
sxemalar yordamida sonli yechishning Pascal dasturi
 
Oshkor sxema. Chap ayirmali sxema. p0<0, pN<0 
uses crt; 
const n=15; j0=20; tt=1; l=1; A=0.01; 
a1=1; q=2; 
type m=array[0..n] of real; 
hi=array[0..n] of real; 
var i,j:integer; 
x,h,t,tau,d:hi; 
u,u1,g,u2,u11,u12:m; 
function ut(p,r:real):real; 
begin ut:= A*exp(p+r);end; 
function fi(p,r:real):real; 
begin 
fi:=A*exp(p+r)*(p*(p+1)+r*(r+1)+7); 
end; 
function ro(p,r:real):real; 
begin ro:=sqr(p)+sqr(r)+5;end; 
function p1(p,r:real):real; 
begin p1:=-(p+r+2);end; 
begin 
clrscr; 
writeln ( 'yugiruvchi hisob sxemasi'); 
writeln(' p0<0, pN<0 hol'); 
writeln(' chap ayirmali sxema'); 
readln; 
h[0]:=0; 
h[1]:=a1; 
for i:=2 to n do 
h[i]:=l/n; 
tau[j]:=tt/j0; 
t[j]:=0;j:=1; 
for i:=0 to n do 
begin 
x[i]:=i*h[i];t[j]:=j*tau[j]; 
u[i]:= A*(exp(x[i])); 
end; 
while t[j]<=tt do begin clrscr; 
t[j]:=t[j]+tau[j]; 
u1[0]:=A*exp((t[j])); 
for i:=n-1 downto 0 do 
begin 
g[i]:=tau[j+1]*p1(x[i],t[j+1])/h[i+1]; 
u11[i]:=(-
g[i]*u1[i+1])+(ro(x[i],t[j+1])*u[i]); 
u12[i]:=tau[j+1]*fi(x[i],t[j]); 
u1[i]:=(u11[i]+u12[i])/(ro(x[i],t[j+1])+
g[i]); 

---

70 
end; 
for i:=n-1 downto 0 do 
u[i]:=u1[i]; 
writeln('---------------------------------'); 
write(' ',j,'qatlam'); 
writeln(' '); 
writeln('----------------------------------'); 
writeln('N Taqribiy yechim Aniq 
yechim Xatolik '); 
writeln('---------------------------------'); 
for i:=0 to n do 
begin 
d[i]:=abs(ut(x[i],t[j])-u1[i]); 
write('',I,' ',u1[i]:6:8,' '); 
writeln(ut(x[i],t[j]):6:8,' ',d[i]:6:8,' '); 
end; 
j:=j+1; 
writeln('-----------------------------------'); 
readln; end; 
end.