ref: 7e00a60abbc505720a4709e00dc7ebe7b6cb2cfc
dir: /sbsi/
#!/bin/rc
#
# Surface-based Body Shape Index (SBSI)
# H(cm)^7/4 * WC(cm)^5/6
# SBSI = ------------------------
# BSA(m²) * VTC(cm)
# where WC is waist circumference and VTC is the vertical trunk one.
#
rfork ne
vtc=()
wc=()
w=()
h=()
flagfmt=''
args='vtc wc weight height'
fn usage {
aux/usage
exit usage
}
if(! ifs=() eval `{aux/getflags $*} || ! ~ $#* 4)
usage
vtc=$1
wc=$2
w=$3
h=$4
bsa=`{bsa $w $h}
hoc <<EOF
($h ^ (7/4) * $wc ^ (5/6))/($bsa * $vtc)
EOF