Пособие по практике "Теория электромагнитного поля"

USES GRAPH; TYPE VECT0R=ARRAY[1..7] OF REAL; INTVEC=ARRAY[1..7] OF INTEGER; VAR AL,C,G,H,II,JJ,LL, MS,MX,OM,R,T,TO,U,WW: REAL; I,J,M,N,XO,XG,XS, GY,GD,GM: INTEGER; D,DD,Y,YO,MA: VECTOR; YG,YS,YO:INTVEC; {$I D:\BP\BIN\KOORSET.PAS} PROCEDURE DER; BEGIN AL:=Y[1]; 0M:=Y[2]; II:=Y[3]; WW:=Y[4]; D[1]:=0M; D[2]:=(C*II-MS)/JJ; D[3]:=(U-C*0M-R*II)/LL; D[4]:=U*II END; PROCEDURE RUNCUT; BEGIN DER; FOR I:=l TO M DO BEGIN DD[I]:=D[I]; YO[I]:=Y[I]; Y[I]:=YO[I]+G*D[I] END; T:=T+G; DER; FOR I:=l TO M DO BEGIN Y[I]:=YO[I]+G*D[I]; DD[I]:=DD[I]+2*D[I] END; DER; FOR I:=l TO M DO BEGIN Y[I]:=YO[I]+H*D[I]; DD[I]:=DD[I]+2*D[I] END; T:=T+G; DER; F0RI:=1 TOM DO Y[I]:=Y0[I]+H/6*(DD[I]+D[I]) END; BEGIN T0:=0; M:=4; N:=1000; H:=0.002; R:=l; C:=2; LL:=0.1; JJ:=1; MS:=2; U:=27; X0:=10; MX:=600/(N*H); 32