jw_cad 外部変形 - (1026) tclshで多角形の面積を計算する -

外部変形は データのやり取りをテキストファイルで行うので プログラム言語は 自由に選ぶことができます。図形は機能的かつシンプルなため、数多くのユーザーに受け入れられています。

 

tclshで多角形の面積を計算する

:tclshで多角形の面積を計算する
@echo off
for /f "delims=:" %%n in ('findstr /n "^#!" %0') do (
  copy jwc_temp.txt myfiles > nul
  more +%%n %0 | tclsh - > jwc_temp.txt
)
goto:eof

REM #jww
REM #1-%d 点を指示してください
REM #99#%d
REM #e

#!この次の行からプログラムを書いてください
set tcl_precision 16; #文字列を数値演算するとき有効桁数 16 (デフォルト 12)

proc area {{a 0} {sx 0} {sy 0} {ix 0} {iy 0}} {
  global hp
  set n [array size hp]
  if {$hp(1) != $hp($n)} {
    set hp([expr $n + 1]) $hp(1)
    incr n
  }
  for {set i 1} {$i < $n} {incr i} {
    lassign "$hp($i) $hp([expr $i + 1])" x1 y1 x2 y2
    set a [expr $a + 0.5 * ($x1 - $x2) * ($y1 + $y2)]
    set sx [expr $sx + ($x1 - $x2) * ($y1 * (2 * $x1 + $x2) + $y2 * (2 * $x2 + $x1)) / 6.0]
    set sy [expr $sy + ($y2 - $y1) * ($x1 * (2 * $y1 + $y2) + $x2 * (2 * $y2 + $y1)) / 6.0]
    set ix [expr $ix + ($x1 - $x2) * ($y1 * pow(2 * $x1 + $x2, 2) + $y2 * pow(2 * $x2 + $x1, 2) + pow($x1 - $x2, 2) * ($y1 + $y2) / 2) / 18.0]
    set iy [expr $iy + ($y2 - $y1) * ($x1 * pow(2 * $y1 + $y2, 2) + $x2 * pow(2 * $y2 + $y1, 2) + pow($y1 - $y2, 2) * ($x1 + $x2) / 2) / 18.0]
  }
  return "$a $sx $sy $ix $iy"
}

foreach 0 [split [read [open myfiles]] \n] {
  switch -regexp -matchvar n -- $0 {
    {^hp(\d+)} {
      lassign $0 1 2 3
      set hpn [lindex $n 1]
      set hp($hpn) "$2 $3"
    }
  }
}

set 0 [area]
lassign $0 a sx sy ix iy
set gx [expr $sx / $a]
set gy [expr $sy / $a]
set ix [expr abs($ix) - abs($a) * $gx * $gx]
set iy [expr abs($iy) - abs($a) * $gy * $gy]
puts [format "h#Area = %.3f Iy = %.3f Ix = %.3f" $a $ix $iy]