Методы расчета диффузионных структур. Быкадорова Г.В - 14 стр.

                                           14

  const an=7.12e8;
          gn=2.3;
          bn=-3.8;
          cn=0.73;
             mmin=53.24;
            non=1.072e17;
  var a,ml:double;
    begin
      a:=1/(1+exp(bn*ln(293/300))*exp(cn*ln(z/non)));
      ml:=an*exp(-gn*ln(293));
mp:=a*ml+mmin*(1-a)
end;
function fn(z:double):double;
  begin
fn:=cs*erfc(z/dt)-ci
end;
       function rs(a,b:double):double;
    var m,s:double;
i1:integer;
             begin
            s:=0.0;
             for i1:=1 to n do if(x[i1]<=b) and (x[i1]>=a) then
            begin
              m:=mp(abs(c[i1]+c[i1-1])/2);
              s:=s+m*(c[i1]+c[i1-1])/2*(x[i1]-x[i1-1])
            end;
              rs:=1/(q*s)
end;
 BEGIN
       writeln(' ');
      write(' к он цен траци яи сх. при м еси в подл ожк е в см -3? ');
       readln(ci);
      write(' тем перату ра ди ф ф у зи и в гр. Ц ел ьси я? ');
       readln(tem);
      tem:=tem+273;
      write(' врем яди ф ф у зи и в м и н у тах? ');
      readln(time);
      time:=time*60;
      write(' предел ьн аяраствори м ость ф осф ора в см -3? ');
      readln(cs);
      d:=dp(tem);
      dt:=2*sqrt(d*time);
      xmax:=4*dt; h:=xmax/200;
      i:=0;
      x[0]:=0;