■サマーヴィルの等面四面体(その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}
かなり面倒なので,中断.
===================================