.MCD 25000 0
.CMD SURFACEFORMAT rot=10 tilt=35 vScale=20 size=15,30
.CMD SKETCHFORMAT mag=1.000000,1.000000 center=0.500000,0.500000 size=15,30 box=y
.CMD PLOTFORMAT logs=0,0 subdivs=1,1 size=5,15 type=l
.CMD FORMAT rd=d ct=10 im=i et=5 zt=15 pr=3 масса длина вpемя заряд
.CMD SET ORIGIN 1
.CMD SET TOL 0.001000
.CMD MARGIN 0
.CMD LINELENGTH 78
.CMD SET PRNCOLWIDTH 8
.CMD SET PRNPRECISION 4
.EQN 0 55 2 13
р:4*10^-6
.EQN 0 12 2 14
E:1.7*10^5
.EQN 1 -67 9 12
X:({9,1}Ў2*lЎlЎ0Ў2*lЎlЎ0Ў2*lЎlЎ0)
.EQN 0 12 9 12
Y:({9,1}Ў2*hЎ2*hЎ2*hЎhЎhЎhЎ0Ў0Ў0)
.EQN 0 12 8 16
Z:({8,3}Ў8Ў7Ў9Ў8Ў5Ў4Ў6Ў5Ў9Ў8Ў6Ў5Ў6Ў5Ў3Ў2Ў5Ў4Ў5Ў4Ў2Ў1Ў2Ў1)
.EQN 0 19 1 8
l~20
.EQN 1 0 1 10
h~12.5
.EQN 0 12 3 10
A:(l*h)/2
.EQN 0 12 1 11
ц:0.27
.EQN 1 -24 1 9
t~0.7
.EQN 0 24 1 11
t1:100
.EQN 1 0 1 10
t2:20
.EQN 2 -20 5 27
D:E/(1-ц^2)*({3,3}Ў(1-ц)/2Ў0Ў0Ў0Ў1ЎцЎ0ЎцЎ1)
.EQN 5 -47 4 45
b(z,i):Y[(Z[(z,mod(i,3)+1){49} )-Y[(Z[(z,mod(i+1,3)+1){49} )
.EQN 1 66 8 11
T:({8,1}Ўt2Ўt2Ўt2Ўt2Ўt2Ўt1Ўt2Ўt1)
.EQN 4 -66 4 45
c(z,i):X[(Z[(z,mod(i+1,3)+1){49} )-X[(Z[(z,mod(i,3)+1){49} )
.EQN 5 0 3 62
B(z):1/(2*A)*({3,6}Ўb(z,3)Ўc(z,3)Ў0Ўc(z,3)Ў0Ўb(z,3)Ўb(z,2)Ўc(z,2)Ў0Ўc(z,2)Ў0Ўb(z,2)Ўb(z,1)Ўc(z,1)Ў0Ўc(z,1)Ў0Ўb(z,1))
.EQN 4 0 1 26
k(z):B(z){51} *D*B(z)*t*A
.EQN 2 0 6 28
f(z):(р*E*T[z*t)/(2*(1-ц))*({6,1}Ўc(z,3)Ўb(z,3)Ўc(z,2)Ўb(z,2)Ўc(z,1)Ўb(z,1))
.EQN 0 34 1 12
r:1;8
.EQN 1 16 2 15
K[(9*2,9*2):0
.EQN 0 16 2 11
Sum[9:0
.EQN 1 -32 1 13
nx:1;3
.EQN 1 0 1 14
ix:-1;0
.EQN 1 0 1 13
ny:1;3
.EQN 0 16 2 11
F[(9*2):0
.EQN 0 13 2 11
ы[(9*2):0
.EQN 1 -29 1 14
iy:-1;0
.EQN 2 -34 3 72
K[(Z[(r,ny)*2+iy,Z[(r,nx)*2+ix):K[(Z[(r,ny)*2+iy,Z[(r,nx)*2+ix)+k(r)[(ny*2+iy,nx*2+ix)
.EQN 3 53 3 18
▀(z):р*T[z*({3,1}Ў0Ў1Ў1)
.EQN 1 -53 3 42
F[(Z[(r,ny)*2+iy):F[(Z[(r,ny)*2+iy)+f(r)[(ny*2+iy)
.EQN 4 0 1 65
({1,18}Ў18Ў17Ў16Ў15Ў14ЎЎ12Ў11Ў10Ў9Ў8ЎЎЎ5ЎЎ3ЎЎ)
.EQN 2 0 1 53
UV:({1,12}Ў18Ў17Ў16Ў15Ў14Ў12Ў11Ў10Ў9Ў8Ў5Ў3){51}
.EQN 1 0 1 20
i:1;(9*2-6)
.EQN 0 24 1 20
j:1;(9*2-6)
.EQN 2 -24 3 19
KK[(i,j):K[(UV[i,UV[j)
.EQN 0 24 4 13
FF[i:F[(UV[i{49} )
.EQN 0 32 2 15
ыы:KK^-1*FF
.EQN 3 0 4 13
ы[(UV[i{49} ):ыы[i
.EQN 1 -56 18 38
хх(z):D*B(z)*({6,1}Ўы[(Z[(z,3)*2)Ўы[(Z[(z,3)*2-1)Ўы[(Z[(z,2)*2)Ўы[(Z[(z,2)*2-1)Ўы[(Z[(z,1)*2)Ўы[(Z[(z,1)*2-1))-D*▀(z)
.EQN 4 46 4 26
Sum[(Z[(r,nx){49} ):Sum[(Z[(r,nx){49} )+1
.EQN 6 -1 1 12
n:1;9
.EQN 2 0 2 12
х[(nx,n):0
.EQN 7 -45 4 34
х[(nx,Z[(r,ny){49} ):х[(nx,Z[(r,ny){49} )+хх(r)[nx
.EQN 0 48 5 16
х[(nx,n):х[(nx,n)/Sum[n