#!/bin/csh


#	dm "x2" "x1" "x1 * y2" "x1 * y3" < trapping.tmp | \
#	regress trap SA SA^2 SA^3 | \
#	sed -e "1,15d" -e "22,24d" -e "20d" | awk '{print $3}' | \
#	transpose  >! reg.parms.tmp

#	awk '{print "dm \""}' '{print $1 * x1+$2 * x1 * x1+$3 * x1 * x1 * x1+$4}' < reg.parms.tmp
#	awk '{print "dm \""}' '{print $1}' < reg.parms.tmp

	set c1 = `awk '{print $1}' < reg.parms.tmp`
	set c2 = `awk '{print $2}' < reg.parms.tmp`
	set c3 = `awk '{print $3}' < reg.parms.tmp`
	set c4 = `awk '{print $4}' < reg.parms.tmp`


gentable 1 90  ${c4} ${c1} ${c2} ${c3}

#	@ i=1
#	while (${i} < 90)
##		echo ${i} ${c1} ${c2} ${c3} ${c4}
##		@ val = `expr ${c1} \* ${i} + ${c2} \* ${i} \* ${i} + ${c3} \* ${i} \* ${i} \* ${i} + ${c4} `
#		echo expr ${c1} \* ${i} + ${c2} \* ${i} \* ${i} + ${c3} \* ${i} \* ${i} \* ${i} + ${c4} 
#		echo `expr ${c1} \* ${i} + ${c2} \* ${i} \* ${i} + ${c3} \* ${i} \* ${i} \* ${i} + ${c4} `
#		echo ${i} ${val}
#		@ i++
#	end
#
#	
##	awk '{ print $1 $2 $3 $4}'<reg.parms.tmp
##	awk '{ for (i=1;i<90;i++) print $i $1 * $i+$2 * $i * $i+$3 * $i * $i * $i+$4}'<reg.parms.tmp
##	awk '{ print "$1 * i+$2 * i * i+$3 * i * i * i+$4"}' <reg.parms.tmp | echo
##	echo "for (i=1;i<10i++) ${c1} * i+${c2} * i * i+${c3} * i * i * i+${c4}" | bc
##	echo "for (i=1;i<10i++) `awk '{ print $1 * i+$2 * i * i+$3 * i * i * i+$4}' <reg.parms.tmp`"
#
#
##set noglob
##bc <<DONE
##for(i=1;i<85;i++) ${c1} * i+${c2} * i * i+${c3} * i * i * i+${c4}
##DONE: >! tmp