- 01
- 02
- 03
- 04
- 05
- 06
- 07
- 08
- 09
- 10
- 11
- 12
- 13
- 14
- 15
- 16
- 17
- 18
- 19
- 20
- 21
- 22
- 23
- 24
- 25
- 26
- 27
- 28
- 29
- 30
- 31
- 32
- 33
- 34
- 35
- 36
- 37
- 38
- 39
- 40
- 41
procedure PGSK_VSK(xia,xtia:dimer3; kw:tnmatrix; var kw_v:tnmatrix);
var
VV,NN,BB:dimer3;
modu:extended;
k, i, m, L: Cardinal;
Mvv:mat33;
begin
modu:=md(xtia);
VV[1]:=xtia[1]/modu; VV[2]:=xtia[2]/modu; VV[3]:=xtia[3]/modu;
NN[1]:=xtia[2]*xia[3]-xia[2]*xtia[3];
NN[2]:=xtia[3]*xia[1]-xia[3]*xtia[1];
NN[3]:=xtia[1]*xia[2]-xia[1]*xtia[2];
modu:=md(NN);
NN[1]:=NN[1]/modu; NN[2]:=NN[2]/modu; NN[3]:=NN[3]/modu;
BB[1]:=NN[2]*VV[3]-NN[3]*VV[2];
BB[2]:=NN[3]*VV[1]-NN[1]*VV[3];
BB[3]:=NN[1]*VV[2]-NN[2]*VV[1];
for k:=1 to 3 do Mvv[1][k]:=VV[k];
for k:=1 to 3 do Mvv[2][k]:=NN[k];
for k:=1 to 3 do Mvv[3][k]:=BB[k];
rm661:=MM0; kw_v:=MM0;
for k:=1 to 3 do for l:=1 to 3 do for m:=1 to 3 do
rm661[k][l]:=rm661[k][l]+Mvv[k][m]*kw[m][l];
for k:=1 to 3 do
for l:=k to 3 do begin
for m:=1 to 3 do
kw_v[k][l]:=kw_v[k][l]+Mvv[k][m]*rm661[l][m];
kw_v[l][k]:=kw_v[k][l];
end;
for k:=1 to 3 do for l:=1 to 3 do for m:=1 to 3 do
rm661[k+3][l+3]:=rm661[k+3][l+3]+Mvv[k][m]*kw[m+3][l+3];
for k:=1 to 3 do
for l:=k to 3 do begin
for m:=1 to 3 do
kw_v[k+3][l+3]:=kw_v[k+3][l+3]+Mvv[k][m]*rm661[l+3][m+3];
kw_v[l+3][k+3]:=kw_v[k+3][l+3];
end;
end;
Мой первый пост - это код преподавателя... я встал и вышел...