Экономический и социальный советE:=q/(4*Pi*r*r*e1);{считаем напряженность,в ней первая проницаемость} D:=q/(4*Pi*r*r);{считаем диэлектрическое смещение по формуле} P:=(e1-e3)*E;{считаем поляризованность по формуле} w:=q*q/(32*Pi*e1*r*r*r*r);{считаем плотность энергии по формуле} writeln(r:4:1,' ',D,' ',E,' ',P,' ',W); r:=r+h;{по радиусу смещамся на шаг} end; for i:=300 to 600 do {другая область диэлектрика} Begin E:=q/(4*Pi*r*r*e2);{считаем напряженность,в ней вторая проницаемость} D:=q/(4*Pi*r*r);{считаем диэлектрическое смещение по формуле} P:=(e2-e3)*E;{считаем поляризованность по формуле} w:=q*q/(32*Pi*e2*r*r*r*r);{считаем плотность энергии по формуле}
r:=r+h; end; for i:=600 to 800 do {вышли за область диэлектрика в вакуум} Begin E:=q/(4*Pi*r*r*e3);{считаем напряженность,в ней проницаемость вакуума} D:=q/(4*Pi*r*r);{считаем диэлектрическое смещение по формуле} P:=(e3-e3)*E;{считаем поляризованность по формуле} w:=q*q/(32*Pi*e3*r*r*r*r);{считаем плотность энергии по формуле}
r:=r+h; end; writeln(); writeln('U');{напряженность=интеграл Е от выбранной точки до бесконечности} u1:=0;{бесконечностью приходится принебрегать приняв ее напряженность =0} U:=0;{приходится выполнять расчет отступая от самой дальней точки к центру} r:=8.0;{выбираем дальнюю точку и сначала считаем U для нее} h:=0.1; for i:=80 downto 60 do {для самой дальней области -вакуума} begin {выполняем интнгрирование по Эллеру кто желает может и по Симпсону} Ey1:=q/(4*Pi*r*r*e3);{начальное значение напряженности} Ey2:=q/(4*Pi*(r-h)*(r-h)*e3);{конечное значение напряженности} u1:=h*Ey1+(Ey2-Ey1)*h/2;{интегрируем по Эллеру} U:=U+u1;{складываем маленькие напряженности в одну большую} r:=r-h;{переходим на следущую точку} writeln(U); end; for i:=60 downto 30 do {для области второго диэлектрика} Begin Ey1:=q/(4*Pi*r*r*e2);{все как в предидущем цикле} Ey2:=q/(4*Pi*(r-h)*(r-h)*e2); u1:=h*Ey1+(Ey2-Ey1)*h/2;{можно применить и другой способ интегрирования} U:=U+u1; r:=r-h; writeln(U); end; for i:=30 downto 20 do {для области первого диэлектрика} Begin Ey1:=q/(4*Pi*r*r*e1); Ey2:=q/(4*Pi*(r-h)*(r-h)*e1); u1:=h*Ey1+(Ey2-Ey1)*h/2; U:=U+u1; r:=r-h; writeln(U); end; for i:=20 downto 0 do {для области проводника} Begin Ey1:=0; Ey2:=0; u1:=h*Ey1+(Ey2-Ey1)*h/2; U:=U+u1; r:=r-h; writeln(U); end; end.
экономический и социальный совет
|