А2 Программа № 2

program lab_1;

uses printer;

label 1,2;

var k,p,ca0,cb0,kr,t,xp,xp1,h,s,x,xx,r,dr,lk,v,VR,sr,ca,l:real;

var i,j: integer;

function f(x:real):real;

var a,b,d,q:real;

begin

a:=k*p/ca0; b:=(cb0-0.5*ca0*x)/(l-0.5*ca0*x);

d:=(l-x)/(l-0.2*x); q:=(l-l/p*sqr(x/(kr*(l-x))))/b;

f:=l/(a*b*d*q);

end;

begin

dr:=0.1; lk:=0.4; v:=7.85E-03; VR:=3.14E-03; sr:=pi*sqr(dr)/4;

p:=l.; t:=693; ca0:=0.2;

writeln (1st,'ca0=',ca0);

{ for i:=l to 6 do begin

t:=t+40;

ca0:=ca0+0.05}; writeln (1st,'t=',t);

x:=0.1;

cb0:=0.21*(l-ca0);

kr:=exp((4905.5/t-4.6455)*ln(10));

if t<673 then k:=exp(2.3*(21.9-14850/t));

if (t>673) and (t<806) then k:=exp(2,3*(1.66-1230/t));

if t>806 then k:=1.361;

xp:=x; {method simple iteration}

l:xpl:=kr/(kr+sqrt((l-0.5*ca0*xp)/(p*(cb0-0,5*ca0*xp))));

if abs(xp-xp1)>0,0001 then begin

xp:=xpl; goto 1;

end;

writeln(lst,'kr=',kr, 'xp=', xp,'k=',k);

x:=0;s:=0; ca:=ca0;

while (x<xp) do begin

ca:=ca0*(l-x);

h:=0.01*x:{integral};

s:=f(0)+f(x);

r:=l/f(x);

xx=0;

for j:=l to 99 do

begin xx:=xx+h;

s:=s+2*f(xx);

end;

s:=h/2*s; l:=s*v/sr; if (l>lk) then goto 2;

writeln (lst,'x=',x,'tau=',s, 'r=',r);

writeln (lst, 'l=',l,'ca=',ca);

x:=x+0.1;

end;

2: {end;}

end.

А3 Программа № 3

program lab_2;

USES PRINTER;

var A,T,P0,R0,SIG,V0,Kl,K2,X,H_sS,XX,TAU,Q,R:REAL;

L,I,J:INTEGER;

FUNCTION F(X:REAL):REAL;

VARQ,Z:REAL;

BEGIN

Q:=(l-X)/(l/P042*K2*X);

Z:=K1/A*SQRT(Q);

F:=l/Z;

END;

BEGIN

T:=1000; P0:=0.1E-4; R0:=260; SIG:=57E+03; V0:=0.041;

A:=1/(SIG*R0*V0);

K1:=EXP(0.001*(LN(0.325E+06)-48000/T));

K2:=EXP(0.001*(LN(0.53E-10)+58500/T));

WRITELN(LST,' A=',A₅'K1=',K1, 'K2=',K2);

FOR L:=l TO 9 DO BEGIN

WRITELN(LST,' PO=',PO);

X:=0; {8:=Р(Х); WRITELN(Т(Х)=*,8);}

FOR I:=1 TO 10 DO

BEGIN

H:=0.01*X;

S:=F(0)+F(X);

XX:=0;

FOR J:=1 TO 99 DO

BEGIN XX:=XX+H;

S:=S+2*F(XX);

END;

S:=H/2*S*3600;

Q:=(l-X)/(l/P0+2*K2*X);

R:=Kl*SQRT(Q)*3600;

WRITELN(LST,' X=',X, 'TAU=',S, 'R=',R);

X:=X+0.1;

END;

P0:=P0+0.1E-4;

END;

END.

А4 Программа № 4

program lab_3;

const m=2;

var k, tau, dtau, са0, сb0, са, сb, сs, сr, z, ха, а, b, с, d, ca1, са2:геа1;

var i,j,N:integer;

var AM, BM, CM, DM:array[l..M] OF REAL;

VAR CAM, CAM1, CAM2, CBM2, CRM, CSM:array[0..M] OF REAL;

begin

N:=2;

k:=0.08;

tau:=600;

ca0:=0.1/N; cb0:=0.1*(N-l)/N;

writeln ('______________________riv____________________');

writeln;

dtau:=0.01*tau;

ca:=ca0; cb:=cb0; cr:=0; cs:=0; t:=0;

writeln('i ca cb cr');

for i:=l to 100 do

begin

z:=k*ca*cb;

ca:=ca-dtau*z;

cb:=cb-dtau*z;

cr:=cr+dtau*z;

cs:=cs+dtau*z;

{t:=t+dtau;}

if (i=10)or(i=20)or(i=30)or(i=40)or(i=50) then writeln(i,ca,cb,cr);

if (i=60)or(i=70)or(i=80)or(i=90)or(i=100) then writeln(i,ca,cb,cr);

end;

xa:=(ca0-ca)/ca0;

writeln ('xa=',xa,'; tau=i*dtau;',' dtau=',dtau);

writeln ('_____________________ ris_____________________');

writeln;

A:=tau*k; B:=(1+tau*k*(cb0-ca0));

с:=-cа0; D:=sqrt(sqr(B)-4*A*C);

cal:=(-B+D)/(2*A);

ca2:=(-B-D)/(2*A);

if cal<0 then ca:=ca2;

ifca2<0 then ca:=ca1;

cb:=cb0/(1+tau*k*ca);

cr:=tau*k*ca*cb;

cs:=сr; ха:=(са0-cа)/са0;

writeln('tau=',tau,' ca=',ca,' cb=',cb);

writeln('cr=',cr,' cs=',cs,' xa=',xa);

writeln('_____________________ris-m __________________');

writeln;

writeln. ('j ca cb cr cs');

cam[0] :=ca0; cbm[0] :=cb0; crm[0] :=0; csm[0] :=0; t:=tau/m;

for j:=l to m do begin

Am[j]:=t*k; Bm[j]:=(l+t*k*(cbm[j-l]-cam[j-l]));

cm[j] :=-cam[j-l]; Dm[j] :=sqrt(sqr(Bm[j])-4*Am[j] *Cm[j]);

caml [j]:=(-Bm[j]+Dm[j])/(2*Am[j]);

cam2[j]:=(-Bm[j]-Dm[j])/(2*Am[J]);

if caml [j]<0 then cam[j]:=cam2[j];

if cam2[j]<0 then cam[j]:=caml[j];

cbm[j]:=cbm[j-l]/(l+t*k*cam[j]);

crm[j] :=t*k*cam[j]*cbm[j];

csm[j]:=crm[j];

writeln (j, cam[j], cbm[j], crm[j], csm[j]);

end;

xa:=(cam[0]-cam[m])/cam[0];

writeln ('m=',m,' xa=',xa);

end.