■サマーヴィルの等面四面体(その30)

E55=−D05+(D15−D05)(D51−D01)/D01

E44=−D04+(D14−D04)(D41−D01)/D01

F55=E55+(D25−D05)(D52−D02)/D02

G55=F55+(D35−D05)(D53−D03)/D03

F44=E44+(D24−D04)(D42−D02)/D02

E33=−D03+(D13−D03)(D31−D01)/D01

E22=−D02+(D12−D02)(D21−D01)/D01

F33=E33−(D23−D03)(D32−D02)/E22

G44=F44−(D34−D04)(D43−D03)/F33

H55=G55−(D45−D05)(D54−D04)/G44

===================================

D01=D12=D23=D34=D45

D02=D13=D24=D35

D03=D14=D25

D04=D15=D02

D05=D01

E55=−D01+(D02−D01)(D02−D01)/D01*

E44=−D02+(D03−D02)(D03−D01)/D01*

F55=E55+(D03−D01)(D03−D02)/D02

=−D01+(D02−D01)(D02−D01)/D01+(D03−D01)(D03−D02)/D02*

G55=F55+(D02−D01)(D02−D03)/D03

=−D01+(D02−D01)(D02−D01)/D01+(D03−D01)(D03−D02)/D02+(D02−D01)(D02−D03)/D03*

F44=E44+(D02−D02)(D02−D02)/D02=E44*

E33=−D03+(D02−D03)(D02−D01)/D01*

E22=−D02+(D01−D02)(D01−D01)/D01=−D02*

F33=E33−(D01−D03)(D01−D02)/E22

=−D03+(D02−D03)(D02−D01)/D01+(D01−D03)(D01−D02)/D02*

G44=F44−(D01−D02)(D01−D03)/F33

H55=G55−(D01−D01)(D01−D01)/G44=G55*

D01E22=−D01D02

F33G44=F33F44−(D01−D02)(D01−D03)

F33G44H55=F33F44G55−(D01−D02)(D01−D03)G55

D01E22F33G44H55=−D01D02F33F44G55+(D01−D02)(D01−D03)G55D01D02

G55D01D02=

=−D01^2D02+(D02−D01)^2D02+(D03−D01)(D03−D02)D01+(D02−D01)(D02−D03)D01D02/D03

G55D01D02(D01−D02)(D01−D03)=

=−D01^2D02(D01−D02)(D01−D03)−(D01−D03)(D02−D01)^3D02−(D03−D01)^2(D01−D02)(D03−D02)D01−(D02−D01)^2(D01−D03)(D02−D03)D01D02/D03

F33F44=

={−D02D03+(D02−D03)(D02−D01)D02/D01+(D01−D03)(D01−D02)}{(D03−D02)(D03−D01)/D01}

=−D02D03(D03−D02)(D03−D01)/D01

−(D02−D03)^2(D02−D01)(D03−D01)D02/D01^2−(D01−D03)^2(D01−D02)(D03−D02)/D01}

 かなり面倒なので,中断.

===================================