Var Rom, pR,pO,kap,tok,u,Ut,co,cr,ir,V,h,er,eo,l,dx :real;

k,n,NN: integer;

{----------------------------------------------------------------------------------------------------}

Procedure kappa; {локальна електропровідність електроліту kap}

BEGIN

kap:=····{функція концентрацій f3(C^R,C^O);

End;

{---------------------------------- ------------------ -------------------------------------- --------}

Procedure AMPER; { локальна густина струму AMP ,A/сm2}

Label 3;

BEGIN

pR :=0; pO:=0; Rom:=h*2/kap;

tl :=0; tp := U*kap/(h*2); tpa:=0; tpk:=0;

for n:=1 to 50 do begin {алгоритм половинного ділення}

AMP:=(tl+tp)/2;

рR:=… f(I, C^R,C^O); {поляризаційна характеристика аноду}

pO:=…. f(I, C^R,C^O); { поляризаційна характеристика катоду }

ir := AMP*h/kap {падіння напруги в електроліті}

Ut:= E + (ir + pR + pO );

if Ut >U then tp:=t; if Ut < U then tl:=t;

Ut := E - (ir pO+ pR+);

if ABS(Ut - U)*2/(Ut + U)<0.001 then goto 3;

end; {n}

3: END;

{--------------------------------------------------------------------------------}

BEGIN {головна програма}

Read CO,CR,U,ER,EO,V,h,L,NN

Dx:=L/(NN+1);

for k:= 1 to NN do begin {рух по координаті довжини х}

x:= x*k;

KAPPA; AMPER;

CR:=CR– AMP*ER*dx/(V*h);

CO:=CO– AMP*EO*dx/(V*h);

writeln (х:10:0, CR:7:3, CO:7:3, AMP:10:3);

END;