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

program lab_3_l;

var a,b,c,d,a1,b1,kl,k2,t,ca0,cb0,xa,xariv,xaris:real;

var tau, dtau, x, y, y1,ca,cb,cr,cs,cq,cd,friv,fris:real;

var i,j,n:integer;

label 1,11,10,100;

function f(x:real):real;

begin

f:=A*exp(3*ln(x))+b*sqr(x)+c*x+d;

end;

begin

tau:=100; t:=tau;

kl:=1.692E-3; k2:=1.708E-3;

ca0:=l; cb0:=l;

writeln ('______________________riv_______________________');

writeln;

dtau:=0.01*tau;

ca:=ca0; cb:=cb0; cr:=0; cq:=0;

writeln('i ca cb cr cq’);

for i:=l to l00 do

begin

ca:=ca-dtau*(kl *ca+k2*ca*cb);

cb:=cb-dtau*k2*ca*cb;

cr:=cr+dtau*k1*ca;

cq:=cq+dtau*k2*ca*cb;

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

if (i=60)or(i=70)or(i=80)or(i=90)or(i= 100) then

writeln(i,ca,cb,cr,cq);

end;

xariv:=(ca0-ca)/ca0;

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

friv:=cr/ca0;

writeln (' xariv=',xariv,' friv=',friv);

writeln ('____________________ris_______________________');

writeln;

a:=sqr(tau)*kl*k2;

b:⁼t*k2+tau*kl+tau*k2*cb0;

c:=l-tau*k2*ca0;

d:=-ca0;

al:=0.0l; bl:=cb0;

1:x:=(al+bl)/2;

y:=f(x);

if abs(al-bl)<0.001 then goto 10;

yl:=f(al);

if y*yl<0 then bl:=x else al:=x; goto 1;

10:ca:=x; cb:=cb0/(l+tau*k2*ca);

cr:=tau*kl*ca; cs:=cr;

cq:=Tau*k2*ca*cb; cd:=cq;

xaris:=(ca0-ca)/ca0;

fris:=cr/ca0;

writeln('ca=',ca,'cb=',cb,'cr=’,cr,'cq=',cq);

writeln (' xaris=',xaris,' fris=',fris);

xa:=0.1; while xa<xariv do begin tau:=tau+l;

a:=sqr(tau)*kl*k2;

b:=t*k2+tau*k1+tau*k2*cb0;

c:=l-tau*k2*ca0;

d:=-ca0;

al:=0.01; b:=eb0;

ll:x:=(al+bl)/2;

y:=f(x);

if abs(al-bl)<0.001 then goto 100;

yl:=f(al);

if y*yl<0 then bl:=x else al:=x; goto 11;

100:ca;=x; cb:=cb0/(l+tau*k2*ca);

cr:=tau*kl*ca; cs:=cr;

cq:=Tau*k2*ca*cb; cd:=cq;

xa:⁼(ca0-ca)/ca0;

end;

fris:=cr/ca0;

writeln (‘________________xa riv = xa ris__________________');

writeln;

writeln ('ха=’,хаriv,' tau riv=',t,' friv=',friv);

writeln (‘ tau ris=',tau,' fris=',fris);

end.

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

program lab_4;

var bn0,t,P,k,cb0,ca0,e,xp,x,xx,h,s,Vr,r:real;

var i,j:integer;

function f(x:real):real;

begin

f:=exp(3*LN(l+e*x))*exp(3*LN(T))/k/(sqr(l-x)*(Cb0-

0.5*Ca0*x)*exp(3 *LN(P))); end;

begin

bn0:=0.757;

T:=273; P:=5; k:=1.4E+4;

Cb0:=0.08; Ca0:=0.09; e:=-0.45; xp:=0.8; x:=0;

while (x<xp) do begin

x-x+0.1;

{T:=298+(12.15*x-35.25)/(7.005-0.143*x);

k:=1.4E+4*exp(1.7*LN(293/T));}

h:=0.01*x;{integral}

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

xx:=0;

for j:=l to 99 do

begin

xx:=xx+h;

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

end;

s:=h/2*s; VR:=bn0*s;

writeln ('x=',x,'Vr=',VR);

end;

end.

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

program lab_5;

var y,y1,a1,b1,cp,ca,cb,ca0,cb0,cr0,kt,tau,dtau,ro:real;

var t,t0,tt,tt0,ttk,xa,k,k0,ft,e,r,H,V,x,vt,cpt,cpa,cpr:real;

var i,j,m:integer;

label 1,10;

function f(vt:real):real;

begin

t:=t0; tt:=ttk; ca:=ca0; cb:=cb0; dtau:=0.01*tau;

for i:=l to 100 do begin

k:=exp(ln(k0)-E/(8.31*t));

{ writeln(‘k=’,k,’ca=’,ca,’t’t,'tt=',tt);}

r:=k*ca;

ca:=ca-dtau*r;

cb:=cb+dtau*r;

t:=t+dtau*(H*r+FT*kt/V*(tt-t))/(cp*ro);

tt:=tt+dtau*(ft*kt/(tau*vt*cpt)*(tt-t));

end;

f:=tt-284;

end;

begin

t0:=301; ca0:=1.8; cb0:=0; cr0:=10.8;

cpa:=62; cpr:=53; ro:=1000;

k0:=1.84E+5; e:=4.49E+4; h:=-448E+4; tt0:=284;

ttk:=293; kt:=0.320;

tau:=360; V:=0.283; FT:=14.16; cpt:=4.190;

cp:⁼ca0*cpa+cr0*cpr;

al:=0,l; bl:=4.5;

l:x:=(al+bl)/2;

y:=f(x);

{ writeln (‘x=’,x,'y=’,y);}

if abs(y)<0.1 then goto 10;

yl:=f(al);

if y*yl<0 then bl:=x else al:=x; goto 1;

10:

writeln (‘t0=',t0);

vt:=x;

t:=t0; tt:=ttk; ca:=ca0; cb:=cb0; dtau:=0.01*tau;

for i:=l to 100 do begin

k:=exp(ln(k0)-E/(8.31*t));

r:=k*ca;

ca:=ca-dtau*r;

cb:=cb+dtau*r;

t:=t+dtau*(H*r+FT*WV*(tt-t))/(cp*ro);

tt:=tt+dtau*(ft*kt/(tau*vt*cpt)*(tt-t));

if (i=10)or(i=30)or(i=50)or(i=100) then writeln(i,ca,t,tt);

end;

xa:=(ca0-ca)/ca0;

writeln(‘vt=',vt,'xa=',xa);

end.