Применение магнитного поля типа «магнитный кокон»

В биологии и медицине.

УДК 577.23, УДК 577.24

Медицина. {14.04.2017 03.00 МСК} Теория, требующая совершенно безопасной и относительно дешёвой экспериментальной проверки, предсказывает, что при помощи особой конфигурации такого магнитного поля, становится возможным лечить болезни суставов, кровеносную, лимфатическую и нервную системы организма за счёт глубинного эндофореза. Я, также, предполагаю, что в определённых условиях такое поле способно содействовать омоложению организма и значительному увеличению продолжительности жизни за счёт значительного продления молодости и зрелости. Известно, что живые организмы это коллоидные растворы с упорядоченным и целенаправленным движением материи (атомов и молекул, взаимодействующих при помощи электрического и магнитного полей) и информации с помощью энергетики. Об этом знают биофизики и специалисты в области информации. Так вот, при помощи такого поля, можно направленно управлять энергетикой и процессами, происходящими в организме на молекулярном и атомарном уровнях, а именно, сшивать повреждённые участки белков, РНК и ДНК и механизма, компенсирующего укорочение теломер (теломеразы).

Приложение 1

Январь 1997 года

Впервые в науке - расчет магнитного поля внутри системы проводников,

эквивалентной по электромагнитным свойствам полому замкнутому проводнику.

Рассчитывается напряженность магнитного поля вдоль координатной оси X.

Поверхностный ток Is, текущий по полому замкнутому проводнику между его полюсами А и В, можно представить в виде "нитей тока" i [И.Е. Тамм, Основы теории электричества, стр. 140], токов, текущих по математическим меридианам сферы от B к А.

Это значительно упрощает программу расчета магнитного поля внутри системы проводников, эквивалентной по электромагнитным свойствам (в необходимом приближении) сферическому полому замкнутому проводнику (ПЗП ).

Закон Био-Савара-Лапласа

в векторной форме

Видно, что в пределах точности машинного счёта, МП внутри ПЗП равно нулю.

Приложение 2

Программа расчёта интенсивности магнитного поля внутри сферического полого замкнутого проводника вдоль координаты Х на языке программирования Quick BASIC

1995 год

'SFERA

SCREEN 12

TIMER ON

LOCATE 12, 1

COLOR 15

PRINT "Calculation of a magnetostatic field inside a Spherical Conductors System, which";

LOCATE 14, 2

PRINT "is equivalent by electromagnetic properties to a Closed and Hollow Conductor."

ON TIMER(7) GOSUB 1

DO WHILE INKEY$ = "": LOOP

CLS

Nl = 150 'number of lines with current in Spherical Conductor System

Nds = 50 'number of elements of current in line of SCS

DIM SG(Nds + 1), CG(Nds + 1)

C$ = " +###.##"

Pi = 3.14159265359#

I = 10 'SCS line current

IL = I * Nl 'current of Linear Conductor

Ndl = Nds 'number of elements of current in Linear Conductor

Rs = 20

L = Rs 'length of Linear Conductor strength

dl = L / Ndl

dg = 2 * Pi / Nds

ds = Rs * dg

FOR j = 0 TO Nds - 1

SG(j) = SIN(j * dg + dg / 2)

CG(j) = COS(j * dg + dg / 2)

NEXT j

LINE (500, 10)-(500, 60), 9

LOCATE 3, 64

COLOR 9

PRINT "LC"

CIRCLE (500, 110), 50, 15

LINE (500, 110)-(550, 110), 8

LOCATE 7, 66

COLOR 8

PRINT "Rs"

COLOR 5

LINE (20, 40)-(320, 40) 'X-axis

LINE (20, 105)-(320, 105) 'Y-axis

LINE (20, 165)-(320, 165) 'Z-axis

LINE (20, 230)-(320, 230) 'SCS module-axis

LINE (20, 300)-(320, 300) 'Linear current module-axis

LINE (320, 10)-(320, 350), 8, , &HFF00

LOCATE 23, 40

COLOR 8

PRINT "Rs (radius of SCS,"

LOCATE 24, 30

PRINT "along which calculates a magnetic field)"

FOR X = 0 TO Rs - 1 STEP 1

HVx = 0

HVy = 0

HVz = 0

Hs = 0

HCx = 0

HCy = 0

HCz = 0

HL = 0

FOR v = 0 TO Nds - 1

HUx = 0

HUy = 0

HUz = 0

FOR u = 0 TO Nl - 1

Rx = X - Rs * SG(v) * COS(u)

Ry = -Rs * SG(v) * SIN(u)

Rz = -Rs * CG(v)

R3 = (SQR(Rx * Rx + Ry * Ry + Rz * Rz)) ^ 3

dsx = -ds * CG(v) * COS(u)

dsy = -ds * CG(v) * SIN(u)

dsz = ds * SG(v)

Hsx = (dsy * Rz - Ry * dsz) / R3

Hsy = (dsz * Rx - Rz * dsx) / R3

Hsz = (dsx * Ry - Rx * dsy) / R3

HUx = HUx + Hsx

HUy = HUy + Hsy

HUz = HUz + Hsz

NEXT u

HVx = HVx + HUx

HVy = HVy + HUy

HVz = HVz + HUz

Hs = SQR(HVx ^ 2 + HVy ^ 2 + HVz ^ 2)

RLx = X

'RLy = 0

RLz = Rs + dl * v + dl / 2

RL3 = (SQR(RLx * RLx + RLz * RLz)) ^ 3

'dlx = 0

'dly = 0

dlz = dl

'Hlx = (dly * RLz - RLy * dlz) / RL3

Hly = (dlz * RLx - RLz * dlx) / RL3

'Hlz = (dlx * RLy - RLx * dly) / RL3

'HCx = HCx + Hlx

HCy = HCy + Hly

'HCz = HCz + Hlz

'HL = SQR(HCx ^ 2 + HCy ^ 2 + HCz ^ 2)

HL = HCy

NEXT v

LOCATE 2, 1

COLOR 14

PRINT "Hsx ="; USING C$; I * HVx

PSET (X * 300 / Rs + 20, HVx + 40)

COLOR 10

LOCATE 6, 1

PRINT "Hsy ="; USING C$; I * HVy

PSET (X * 300 / Rs + 20, HVy + 105)

COLOR 11

LOCATE 10, 1

PRINT "Hsz ="; USING C$; I * HVz

PSET (X * 300 / Rs + 20, HVz + 165)

COLOR 15

LOCATE 14, 1

PRINT "Hs ="; USING C$; I * Hs

PSET (X * 300 / Rs + 20, Hs + 230)

LOCATE 14, 25

PRINT "Magnetic intensity in Spherical Conductors System (SCS)"

'LOCATE 14, 25

'PRINT "in Spherical Conductors System (SCS)"

COLOR 9

LOCATE 18, 1

PRINT "HL ="; USING C$; HL * IL * 10

PSET (X * 300 / Rs + 20, HL * IL * 10 + 300)

LOCATE 18, 25

PRINT "Magnetic intensity from Linear Conductor current (LC)"

NEXT X

FOR X = Rs TO Rs * 2 STEP 1

HCx = 0

HCy = 0

HCz = 0

HL = 0

FOR v = 0 TO Nds - 1

RLx = X

'RLy = 0

RLz = Rs + dl * v + dl / 2

RL3 = (SQR(RLx * RLx + RLz * RLz)) ^ 3

'dlx = 0

'dly = 0

dlz = dl

'Hlx = (dly * RLz - RLy * dlz) / RL3

Hly = (dlz * RLx - RLz * dlx) / RL3

'Hlz = (dlx * RLy - RLx * dly) / RL3

'HCx = HCx + Hlx

HCy = HCy + Hly

'HCz = HCz + Hlz

'HL = SQR(HCx ^ 2 + HCy ^ 2 + HCz ^ 2)

HL = HCy

NEXT v

PSET (X * 300 / Rs + 20, HL * IL * 10 + 300)

NEXT X

STOP

COLOR 15

END

Приложение 3

Движение одной заряженной частицы в магнитном поле, создаваемом двумя другими движущимися частицами и между ними.

SCREEN 12

TIMER ON

COLOR 3

LOCATE 2, 12

PRINT " Motion of one charged corpuscle (OMC)"

LOCATE 4, 12

PRINT " in a magnetic field (MF) such as a Magnetic; Cocoon,"

LOCATE 6, 12

PRINT "which is created by two others mobile charged corpuscles (TMC)"

LOCATE 8, 25

PRINT "Etmc = 5 eV; Eomc = 50 eV"

ON TIMER(5) GOSUB 1

DO WHILE INKEY$ = "": LOOP

CLS

C$ = " +#.#####^^^^"

Pi = 3.141593

Mo = 4 * Pi * .0000001

K = Mo / (4 * Pi) 'SI

M = 1.67E-27 * 2 '[kg]

Q = 1.6E-19 '[C]

LOCATE 6, 2

COLOR 9

INPUT "distance between TMC [sm](10 - 50)="; A0

a = A0 * .01 '[m] distance between TMC

Xk1 = 0

Yk1 = 0

Xk2 = Xk1 + a

Yk2 = 0

CIRCLE (Xk1 * 1000 + 50, Yk1 + 240), 10, 9

CIRCLE (Xk2 * 1000 + 50, Yk2 + 240), 10, 9

Ek = 5 '[eV] energy of TMC

TVk = 0 'angle of moving of TMC

Vk = SQR(2 * Ek * ABS(Q) / M) 'velocity of TMC

Vkx = Vk * COS(TVk)

Vky = 0

Vkz = 0

LOCATE 1, 60

PRINT "V tmc ="; USING C$; Vk

COLOR 13

Eq = 50 '[eV] energy of entry of OMC

LOCATE 8, 2

INPUT "angle of entry of OMC (0 - 360)="; TVq

Vq = SQR(2 * Eq * ABS(Q) / M) 'velocity of OMC

Vqx = Vq * COS(TVq)

Vqy = Vq * SIN(TVq)

Vqz = 0

LOCATE 3, 60

PRINT "V omc ="; USING C$; Vq

dlq = .0001# 'step of moving of OMC

Xq = Xk1 + .5 * a 'coordinate began movements of OMC

LOCATE 10, 2

INPUT "Coordinate Y began movements of OMC [sm](-25/+25)="; Yq0

COLOR 8

Yq = Yq0 * .01

Zq = 0

Fx = 0

Fy = 0

Fz = 0

FOR t = 0 TO 100000

dt = dlq / Vq

Rx1 = Xq - Xk1

Rx2 = Xq - Xk2

Ry1 = Yq - Yk1

Ry2 = Yq - Yk2 'Rz1 = 0 Rz2 = 0

R1 = SQR(Rx1 * Rx1 + Ry1 * Ry1)

R2 = SQR(Rx2 * Rx2 + Ry2 * Ry2)

'B = K * Q * [V x R] / R^3 MF, created by TMC [Tl]

'Љ—1

'Bx1 = K * Q * (Vky * Rz1 - Ry1 * Vkz) / R1 ^ 3 == 0

'By1 = K * Q * (Vkz * Rx1 - Rz1 * Vkx) / R1 ^ 3 == 0

Bz1 = K * Q * (Vkx * Ry1 - Rx1 * Vky) / R1 ^ 3

'Љ—2

'Bx2 = K * Q * (Vky * Rz2 - Ry2 * Vkz) / R2 ^ 3 == 0

'By2 = K * Q * (Vkz * Rx2 - Rz2 * Vkx) / R2 ^ 3 == 0

Bz2 = K * Q * (Vkx * Ry2 - Rx2 * Vky) / R2 ^ 3

'Lorentz force which influence on OMC FL = Q*[V x B]

'Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)

'Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)

'Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)

Fx = Vqy * (Bz1 + Bz2)

Fy = -(Bz1 + Bz2) * Vqx

'Fz = (Vqx * (By1 + By2) - (Bx1 + Bx2) * Vqy) ==0

'F = M*dV/dt = FL = Q*[V x Bk] ==> dV=PX*dt/M

dVqx = Fx * dt / M

dVqy = Fy * dt / M

dVqz = 0 'dVz = Fz * dt / M

Vqxt = Vqx + dVqx

Vqyt = Vqy + dVqy

Vqzt = 0 'Vqzt = Vqz + dVqz

Vqt = SQR(Vqxt ^ 2 + Vqyt ^ 2) '+ Vqzt ^ 2)

kV = Vq / Vqt

'kV - because velocity of OMC must be constant

Vqx = Vqxt * kV

Vqy = Vqyt * kV

Vqz = 0 'Vqz = Vqzt * kV

V = SQR(Vqx ^ 2 + Vqy ^ 2) '+ Vqz ^ 2)

dXq = Vqx * dt

dYq = Vqy * dt

dZq = 0 'dZq = Vz * dt

Xq = Xq + dXq

Yq = Yq + dYq

Zq = 0 'Zq = Zq + dZq

PSET (Xq * 1000 + 50, Yq * 1000 + 240), 13

NEXT t

END

Приложение 4

Движение заряженной частицы в магнитном поле типа «Магнитный кокон».

SCREEN 12

TIMER ON

COLOR 3

LOCATE 3, 17

PRINT " Motion of one deuteron in a magnetic field "

LOCATE 5, 17

PRINT " such as a Magnetic Cocoon "

LOCATE 7, 17

PRINT " Ed = 500 eV"

COLOR 15

ON TIMER(3) Gosub 1

DO While INKEY$ = "": LOOP

CLS

Dim X1(80), X2(80)

C$ = " +#.#####^^^^"

Pi = 3.141593

Mo = 4 * Pi * .0000001

M = 1.67E-27 * 2

Q = 1.6E-19

E = 500 '[eV]

2 CLS

LOCATE 2, 10

INPUT "ANGLE (50 - 175)="; AN

If 175 < AN Goto 2

If AN < 50 Goto 2

TV = AN * Pi / 180

V = Sqr(2 * E * Abs(Q) / M)

Vx = V * Cos(TV)

Vy = V * Sin(TV)

Vz = 0

LOCATE 2, 60

PRINT "V ="; USING C$; V

dlq = .0002

'dt = dlq / V

A = .04

YL1 = 0

YL2 = 0

Xv = 0

Yv = .001

Zv = 0

L = .08

DL = .001

N = L / DL

IL = 3000000

K = Mo / (4 * Pi)

For I = 0 To N - 1

X1(I) = -I * DL - DL / 2 - A

X2(I) = I * DL + DL / 2 + A

Next I

For t = 0 To 3000

dt = dlq / V

PXx = 0

PXy = 0

PXz = 0

For I = 0 To N - 1

Rx1 = Xv - X1(I)

Rx2 = Xv - X2(I)

Ux = DL

Uy = 0

Uz = 0

Ry1 = Yv + YL1

Ry2 = Yv + YL2

Rz = Zv

R1 = Sqr(Rx1 * Rx1 + Ry1 * Ry1 + Rz * Rz)

R2 = Sqr(Rx2 * Rx2 + Ry2 * Ry2 + Rz * Rz)

Hy1 = (-Rz * Ux) / R1 ^ 3

Hz1 = (Ux * Ry1) / R1 ^ 3

Hy2 = (-Rz * Ux) / R2 ^ 3

Hz2 = (Ux * Ry2) / R2 ^ 3

Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)

Py = -(Hz1 + Hz2) * Vx

Pz = Vx * (Hy1 + Hy2)

PXx = PXx + Px

PXy = PXy + Py

PXz = PXz + Pz

Next I

dVx = Q * K * IL * PXx * dt / M

dVy = Q * K * IL * PXy * dt / M

dVz = Q * K * IL * PXz * dt / M

Vxt = Vx + dVx

Vyt = Vy + dVy

Vzt = Vz + dVz

Vt = Sqr(Vxt ^ 2 + Vyt ^ 2 + Vzt ^ 2)

kV = V / Vt

'kV - the particle must have a constant value modulo of velocity

Vx = Vxt * kV

Vy = Vyt * kV

Vz = Vzt * kV

V = Sqr(Vx ^ 2 + Vy ^ 2 + Vz ^ 2)

dXv = Vx * dt

dYv = Vy * dt

dZv = Vz * dt

Xv = Xv + dXv

Yv = Yv + dYv

Zv = Zv + dZv

PSET (Xv * 10000 + 320, Yv * 10000 + 240), 12

Next t

End

Приложение 5

Движение двух заряженных частиц в магнитном поле типа «Магнитный кокон» с учётом Кулоновского взаимодействия.

SCREEN 12

TIMER ON

COLOR 3

LOCATE 3, 17

PRINT " Motion of two deuterones in a magnetic field "

LOCATE 5, 17

PRINT " such as a Magnetic Cocoon "

COLOR 15

ON TIMER(3) GOSUB 1

DO WHILE INKEY$ = "": LOOP

CLS

DIM XP1(150), XP2(150)

C$ = " #.###^^^^"

Pi = 3.141593

Eo = 8.85E-12

Mo = 4 * Pi * .0000001

M = 1.67E-27 * 2

Q = 1.6E-19

E = 700! '[eV]

J = 8E+08

LOCATE 3, 10

INPUT "ANGLE (30 - 170)="; AN

IF 170 < AN GOTO 1

IF AN < 30 GOTO 1

TV = AN * Pi / 180

V = SQR(2 * E * ABS(Q) / M)

Vx1 = V * COS(TV)

Vy1 = V * SIN(TV)

Vz1 = 0

Vx2 = Vx1

Vy2 = -Vy1

dlq = .0001

A = .007

Xv1 = 0

Yv1 = .0003

Xv2 = 0

Yv2 = -.0003

L = .25

DL = .01

N = L / DL

Du = .06

S = Pi * Du ^ 2 / 4

IL = S * J

K = Mo / (4 * Pi)

FOR I = 0 TO N - 1

XP1(I) = -I * DL - DL / 2 - A

XP2(I) = I * DL + DL / 2 + A

NEXT I

FOR t = 0 TO 8000

IL = S * J

K1 = 5000000!

dt = dlq / V

LOCATE 2, 65

PRINT "t ="; USING C$; t * dt

PXx1 = 0

PXy1 = 0

PXx2 = 0

PXy2 = 0

PCx1 = 0

PCy1 = 0

PCx2 = 0

PCy2 = 0

FOR I = 0 TO N - 1

Rx11 = Xv1 - XP1(I)

Rx12 = Xv1 - XP2(I)

Rx21 = Xv2 - XP1(I)

Rx22 = Xv2 - XP2(I)

Ux = DL

REM Uy = 0 Uz = 0

Ry1 = Yv1

Ry2 = Yv2

R11 = SQR(Rx11 * Rx11 + Ry1 * Ry1)

R12 = SQR(Rx12 * Rx12 + Ry1 * Ry1)

R21 = SQR(Rx21 * Rx21 + Ry2 * Ry2)

R22 = SQR(Rx22 * Rx22 + Ry2 * Ry2)

RGx1 = Xv1 - Xv2

RGy1 = Yv1 - Yv2

RG1 = SQR(RGx1 ^ 2 + RGy1 ^ 2)

RGx2 = -RGx1

RGy2 = -RGy1

RG2 = RG1

REM Hx1 = (Uy * Rz - Ry * Uz) / R1 ^ 3

REM Hy1 = (Uz * Rx1 - Rz * Ux) / R1 ^ 3

REM Hz1 = (Ux * Ry - Rx1 * Uy) / R1 ^ 3

Hz11 = (Ux * Ry1) / R11 ^ 3

Hz12 = (Ux * Ry1) / R12 ^ 3

Hz21 = (Ux * Ry2) / R21 ^ 3

Hz22 = (Ux * Ry2) / R22 ^ 3

REM Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)

REM Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)

REM Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)

PGx1 = RGx1 / (ABS(RG1 ^ 3)) * K1

PGy1 = RGy1 / (ABS(RG1 ^ 3)) * K1

PGx2 = RGx2 / (ABS(RG2 ^ 3)) * K1

PGy2 = RGy2 / (ABS(RG2 ^ 3)) * K1

'Coefficient "K1" imitates group of particles

Px1 = Vy1 * (Hz11 + Hz12)

Py1 = -(Hz11 + Hz12) * Vx1

Px2 = Vy2 * (Hz21 + Hz22)

Py2 = -(Hz21 + Hz22) * Vx2

PCx1 = PCx1 + PGx1

PCy1 = PCy1 + PGy1

PCx2 = PCx2 + PGx2

PCy2 = PCy2 + PGy2

PXx1 = PXx1 + Px1

PXy1 = PXy1 + Py1

PXx2 = PXx2 + Px2

PXy2 = PXy2 + Py2

NEXT I

dVx1 = Q * K * IL * PXx1 * dt / M

dVy1 = Q * K * IL * PXy1 * dt / M

dVx2 = Q * K * IL * PXx2 * dt / M

dVy2 = Q * K * IL * PXy2 * dt / M

Vxt1 = Vx1 + dVx1

Vyt1 = Vy1 + dVy1

Vt1 = SQR(Vxt1 ^ 2 + Vyt1 ^ 2)

kV1 = V / Vt1

Vxt2 = Vx2 + dVx2

Vyt2 = Vy2 + dVy2

Vt2 = SQR(Vxt2 ^ 2 + Vyt2 ^ 2)

kV2 = V / Vt2

'kV - the particle must have a constant value modulo of velocity

dVCx1 = Q ^ 2 * PCx1 * dt / M / (4 * Pi * Eo)

dVCy1 = Q ^ 2 * PCy1 * dt / M / (4 * Pi * Eo)

dVCx2 = Q ^ 2 * PCx2 * dt / M / (4 * Pi * Eo)

dVCy2 = Q ^ 2 * PCy2 * dt / M / (4 * Pi * Eo)

Vx1 = Vxt1 * kV1 + dVCx1

Vy1 = Vyt1 * kV1 + dVCy1

Vx2 = Vxt2 * kV2 + dVCx2

Vy2 = Vyt2 * kV2 + dVCy2

Vt = SQR(Vx1 ^ 2 + Vy1 ^ 2 + Vz1 ^ 2)

REM LOCATE 23, 20

REM PRINT "Vt ="; USING C$; Vt

dXv1 = Vx1 * dt

dYv1 = Vy1 * dt

dXv2 = Vx2 * dt

dYv2 = Vy2 * dt

Xv1 = Xv1 + dXv1

Yv1 = Yv1 + dYv1

Xv2 = Xv2 + dXv2

Yv2 = Yv2 + dYv2

PSET (Xv1 * 70000 + 320, Yv1 * 60000 + 240), 9

PSET (Xv2 * 70000 + 320, Yv2 * 60000 + 240), 12

IF t = 2000 THEN CLS 1

IF t = 5000 THEN CLS 1

IF t = 6500 THEN CLS 1

NEXT t

LOCATE 29, 5

INPUT "repeat under other angle - 1 "; U1

IF U1 = 1 THEN GOTO 1

IF U1 <> 1 THEN END

Приложение 6

Движение двух заряженных частиц в нарастающем магнитном поле типа «Магнитный кокон» - «Сферический пинч».

SCREEN 12

TIMER ON

COLOR 3

LOCATE 3, 12

PRINT " Motion of two deuterones"

LOCATE 5, 12

PRINT " in a increasing magnetic field such as a Magnetic Cocoon "

LOCATE 7, 12

PRINT " (Spherical pinch)"

COLOR 15

ON TIMER(4) GOSUB 1

DO WHILE INKEY$ = "": LOOP

CLS

DIM XP1(80), XP2(80)

C$ = " ##.###^^^^"

Pi = 3.141593

Eo = 8.85E-12

Mo = 4 * Pi * .0000001

M = 1.67E-27 * 2

Q = 1.6E-19

E = 10000 '[eV]

J = 7E+10

LOCATE 2, 10

INPUT "ANGLE (60 - 300)="; AN

IF 300 < AN GOTO 1

IF AN < 60 GOTO 1

TV = AN * Pi / 180

V = SQR(2 * E * ABS(Q) / M)

Vx1 = V * COS(TV)

Vy1 = V * SIN(TV)

Vz1 = 0

Vx2 = Vx1

Vy2 = -Vy1

dlq = .0003

A = .015

Xv1 = 0

Yv1 = .001

Xv2 = 0

Yv2 = -.001

L = .08

DL = .001

N = L / DL

Du = .01

s = Pi * Du ^ 2 / 4

IL = s * J

K = Mo / (4 * Pi)

FOR i = 0 TO N - 1

XP1(i) = -i * DL - DL / 2 - A

XP2(i) = i * DL + DL / 2 + A

'PSET ((XP1(I) * 2000 + 320 - A * 5000), 240)

'PSET ((XP2(I) * 2000 + 320 + A * 5000), 240)

NEXT i

tt = 0

FOR t = 0 TO 8000

dt = dlq / V

tt = tt + dt

LOCATE 2, 65

PRINT "t ="; USING C$; tt

LOCATE 4, 65

PRINT "I ="; USING C$; IL

PXx1 = 0

PXy1 = 0

PXx2 = 0

PXy2 = 0

PCx1 = 0

PCy1 = 0

PCx2 = 0

PCy2 = 0

FOR i = 0 TO N - 1

Rx11 = Xv1 - XP1(i)

Rx12 = Xv1 - XP2(i)

Rx21 = Xv2 - XP1(i)

Rx22 = Xv2 - XP2(i)

Ux = DL

REM Uy = 0 Uz = 0

Ry1 = Yv1

Ry2 = Yv2

R11 = SQR(Rx11 * Rx11 + Ry1 * Ry1)

R12 = SQR(Rx12 * Rx12 + Ry1 * Ry1)

R21 = SQR(Rx21 * Rx21 + Ry2 * Ry2)

R22 = SQR(Rx22 * Rx22 + Ry2 * Ry2)

RGx1 = Xv1 - Xv2

RGy1 = Yv1 - Yv2

RG1 = SQR(RGx1 ^ 2 + RGy1 ^ 2)

RGx2 = -RGx1

RGy2 = -RGy1

RG2 = RG1

REM Hx1 = (Uy * Rz - Ry * Uz) / R1 ^ 3

REM Hy1 = (Uz * Rx1 - Rz * Ux) / R1 ^ 3

REM Hz1 = (Ux * Ry - Rx1 * Uy) / R1 ^ 3

Hz11 = (Ux * Ry1) / R11 ^ 3

Hz12 = (Ux * Ry1) / R12 ^ 3

Hz21 = (Ux * Ry2) / R21 ^ 3

Hz22 = (Ux * Ry2) / R22 ^ 3

REM Px = (Vy * (Hz1 + Hz2) - (Hy1 + Hy2) * Vz)

REM Py = (Vz * (Hx1 + Hx2) - (Hz1 + Hz2) * Vx)

REM Pz = (Vx * (Hy1 + Hy2) - (Hx1 + Hx2) * Vy)

PGx1 = RGx1 / (ABS(RG1 ^ 3))

PGy1 = RGy1 / (ABS(RG1 ^ 3))

PGx2 = RGx2 / (ABS(RG2 ^ 3))

PGy2 = RGy2 / (ABS(RG2 ^ 3))

Px1 = Vy1 * (Hz11 + Hz12)

Py1 = -(Hz11 + Hz12) * Vx1

Px2 = Vy2 * (Hz21 + Hz22)

Py2 = -(Hz21 + Hz22) * Vx2

PCx1 = PCx1 + PGx1

PCy1 = PCy1 + PGy1

PCx2 = PCx2 + PGx2

PCy2 = PCy2 + PGy2

PXx1 = PXx1 + Px1

PXy1 = PXy1 + Py1

PXx2 = PXx2 + Px2

PXy2 = PXy2 + Py2

NEXT i

dVx1 = Q * K * IL * PXx1 * dt / M

dVy1 = Q * K * IL * PXy1 * dt / M

dVx2 = Q * K * IL * PXx2 * dt / M

dVy2 = Q * K * IL * PXy2 * dt / M

Vxt1 = Vx1 + dVx1

Vyt1 = Vy1 + dVy1

Vt1 = SQR(Vxt1 ^ 2 + Vyt1 ^ 2)

kV1 = V / Vt1

Vxt2 = Vx2 + dVx2

Vyt2 = Vy2 + dVy2

Vt2 = SQR(Vxt2 ^ 2 + Vyt2 ^ 2)

kV2 = V / Vt2

'kV - the particle must have a constant value modulo of velocity

kq = 3000 'Magnification of distance of a Coulomb interaction

'imitates group of particles

dVCx1 = (kq * Q) ^ 2 * PCx1 * dt / M / (4 * Pi * Eo)

dVCy1 = (kq * Q) ^ 2 * PCy1 * dt / M / (4 * Pi * Eo)

dVCx2 = (kq * Q) ^ 2 * PCx2 * dt / M / (4 * Pi * Eo)

dVCy2 = (kq * Q) ^ 2 * PCy2 * dt / M / (4 * Pi * Eo)

Vx1 = Vxt1 * kV1 + dVCx1

Vy1 = Vyt1 * kV1 + dVCy1

Vx2 = Vxt2 * kV2 + dVCx2

Vy2 = Vyt2 * kV2 + dVCy2

Vt = SQR(Vx1 ^ 2 + Vy1 ^ 2 + Vz1 ^ 2)

REM LOCATE 23, 20

REM PRINT "Vt ="; USING C$; Vt

dXv1 = Vx1 * dt

dYv1 = Vy1 * dt

dXv2 = Vx2 * dt

dYv2 = Vy2 * dt

Xv1 = Xv1 + dXv1

Yv1 = Yv1 + dYv1

Xv2 = Xv2 + dXv2

Yv2 = Yv2 + dYv2

PSET (Xv1 * 25000 + 320, Yv1 * 25000 + 240), 9

PSET (Xv2 * 25000 + 320, Yv2 * 25000 + 240), 4

Y12 = Yv1 - Yv2

IF Y12 < .0000000000002# THEN dlq = .00000000000001#

IF Y12 > .0000000000002# THEN dlq = .00000000000002#

IF Y12 > .000000000002# THEN dlq = .0000000000002#

IF Y12 > .00000000002# THEN dlq = .000000000002#

IF Y12 > .0000000002# THEN dlq = .00000000002#

IF Y12 > .000000002# THEN dlq = .0000000002#

IF Y12 > .00000002# THEN dlq = .000000002#

IF Y12 > .0000002# THEN dlq = .00000002#

IF Y12 > .000002 THEN dlq = .0000002#

IF Y12 > .00002 THEN dlq = .000004#

IF Y12 > .0002 THEN dlq = .00004#

IF Y12 > .002 THEN dlq = .0001#

IF t = 7000 THEN CLS 1

J = J + 8E+07

IL = s * J

NEXT t

LOCATE 29, 5

INPUT "Repeat with other angle - 1 "; U1

IF U1 = 1 THEN GOTO 1

IF U1 <> 1 THEN END

Приложение 7

Движение электрона в магнитном поле тора с полоидальным током – «Ондулятор».

SCREEN 12

COLOR 3

LOCATE 5, 5

PRINT " Motion of an electron in a magnetic field "

LOCATE 7, 5

PRINT " such as a Magnetic Cocoon (undulator) "

COLOR 15

C$ = " +#.#####^^^^"

Pi = 3.14159265359#

M = 9.1E-31

Q = 1.6E-19

jc = 6E+08

dt = 3E-10

Uq = 400!

Du = .001

TV = 30 * Pi / 180

K = 4 * Pi * .0000001 * (-Q) * (jc * Du ^ 2 / 4) / (4 * Pi)

V = SQR(2 * Q * Uq / M)

Vx = V * COS(TV)

Vy = V * SIN(TV)

Vz = 0

LOCATE 3, 60

PRINT "V= "; USING C$; V

A = .03

X1 = -.05

X2 = .05

dX = (X2 - X1) / 100

dG = Pi / 50

Xv = -.32

Yv = 0

Zv = 0

PSET (Xv * 1200 + 420, Yv * 4000 + 300)

DIM X(100)

FOR i = 0 TO 99

X(i) = X1 + i * dX + dX / 2

PSET (X(i) * 1200 + 420, A * 4000 + 300)

PSET (X(i) * 1200 + 420, -A * 4000 + 300)

NEXT i

PGx = 0

PGy = 0

PGz = 0

FOR t = 1 TO 100

FOR j = 0 TO 99

CG = COS(j * dG)

SG = SIN(j * dG)

PXx = 0

PXy = 0

PXz = 0

Ry = -A * CG + Yv

Rz = -A * SG + Zv

Ry2 = Ry * Ry

Rz2 = Rz * Rz

Ux = dX 'Uy = 0 Uz = 0

FOR i = 0 TO 99

Rx = -X(i) + Xv 'R = SQR(Rx * Rx + Ry2 + Rz2)

R3 = (SQR(Rx * Rx + Ry2 + Rz2)) ^ 3

Hx = (Vy * Rz - Vz * Ry) / R3

Hy = (Vz * Rx - Vx * Rz) / R3

Hz = (Vx * Ry - Vy * Rx) / R3

'Px = Uy * Hz - Uz * Hy

Py = -Ux * Hz 'Py = Uz * Hx - Ux * Hz

Pz = Ux * Hy 'Pz = Ux * Hy - Uy * Hx

'PXx = PXx + Px

PXy = PXy - Py

PXz = PXz - Pz

NEXT i

'PGx = PGx + PXx

PGy = PGy + PXy

PGz = PGz + PXz

NEXT j

'dVx = K * PGx * dt / M

dVy = K * PGy * dt / M

dVz = K * PGz * dt / M

Vxt = Vx + dVx

Vyt = Vy + dVy

Vzt = Vz + dVz

Vt = SQR(Vxt ^ 2 + Vyt ^ 2 + Vzt ^ 2)

kV = V / Vt

'kV - the particle must have a constant value modulo of velocity

Vx = Vxt * kV

Vy = Vyt * kV

Vz = Vzt * kV

V = SQR(Vx ^ 2 + Vy ^ 2 + Vz ^ 2)

dXv = Vx * dt

dYv = Vy * dt

dZv = Vz * dt

Xv = Xv + dXv

Yv = Yv + dYv

Zv = Zv + dZv

PSET (Xv * 1200 + 420, Yv * 4000 + 300), 12

NEXT t

Приложение 8

Дата добавления: 2018-09-20; просмотров: 365; Мы поможем в написании вашей работы!

Поделиться с друзьями:

⇐ Предыдущая 1 2 3 4 5 678 Следующая ⇒

Мы поможем в написании ваших работ!