Компьютерные технологии в физике. Часть 1. Компьютерное моделирование физических процессов. Красов В.И - 86 стр.

UptoLike

89
Bfield(x1-delta1,y0-y1);
fl1:=fl1+2*Pi*y1*Bx*delta/2;
y1:=y1+delta;
Bfield(x1-delta1,y0-y1);
fl1:=fl1+2*Pi*y1*Bx*delta/2;
end;
{вычисление потока в области I}
if (x<x1-delta1)or(x>x2+delta1) then
begin
y1:=0;
delta:=(y0-y)/N;
for i:=1 to N do
begin
Bfield(x,y0-y1);
fl:=fl+2*Pi*y1*Bx*delta/2;
y1:=y1+delta;
Bfield(x,y0-y1);
fl:=fl+2*Pi*y1*Bx*delta/2;
end
end;
else
{вычисление потока в области II }
begin
x3:=x1-delta1;
delta:=(x-x3)/N;
fl2:=0;
{поток по боковой поверхности цилиндра}
for i:=1 to N do
begin
Bfield(x3,y0-R-delta1);
fl2:=fl2+2*Pi*(R+delta1)*By*delta/2;
x3:=x3+delta;
Bfield(x3,y0-R-delta1);
fl2:=fl2+2*Pi*(R+delta1)*By*delta/2;
end;
{поток по кольцу}
y1:=R+delta1;
delta:=(y0-y)/N;
fl3:=0;
while y1<=y0-y do
begin
Bfield(x,y0-y1);
PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com
                    Bfield(x1-delta1,y0-y1);
                    fl1:=fl1+2*Pi*y1*Bx*delta/2;
                    y1:=y1+delta;
                    Bfield(x1-delta1,y0-y1);
                    fl1:=fl1+2*Pi*y1*Bx*delta/2;
                  end;
                    {вычисление потока в области I}
                if (xx2+delta1) then
                  begin
                    y1:=0;
                    delta:=(y0-y)/N;
                    for i:=1 to N do
                      begin
                       Bfield(x,y0-y1);
                       fl:=fl+2*Pi*y1*Bx*delta/2;
                       y1:=y1+delta;
                       Bfield(x,y0-y1);
                       fl:=fl+2*Pi*y1*Bx*delta/2;
                      end
                  end;
                else
                    {вычисление потока в области II }
                  begin
                    x3:=x1-delta1;
                    delta:=(x-x3)/N;
                    fl2:=0;
                   {поток по боковой поверхности цилиндра}
                    for i:=1 to N do
                      begin
                       Bfield(x3,y0-R-delta1);
                       fl2:=fl2+2*Pi*(R+delta1)*By*delta/2;
                       x3:=x3+delta;
                       Bfield(x3,y0-R-delta1);
                       fl2:=fl2+2*Pi*(R+delta1)*By*delta/2;
                      end;
                   {поток по кольцу}
                    y1:=R+delta1;
                    delta:=(y0-y)/N;
                    fl3:=0;
                    while y1<=y0-y do
                      begin
                       Bfield(x,y0-y1);

                                                     89


PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com