#!/usr/bin/perl -w

$pi = 3.1415926;
$window = 500;

for ($i=-$window; $i<$window; $i++) {
	$x = $i / 10.0;
	if ($i != 0) {
		$sinc[$i+$window] = sin($pi * $x) / ($pi * $x);
	} else {
		$sinc[$i+$window] = 1.0;
	}
}

for ($i=0; $i<400; $i++) {
	$out[$i] = 0.0;
}

while (<>) {
	chomp;
	if (/^(.*)\t(.*)$/) {
		$t = $1 + 0;
		$a = $2 + 0;
		for ($i = -$window; $i < $window; $i++) {
			if ($t*10 + $i > 0) {
				$out[$t*10 + $i] += $a * $sinc[$i+$window];
			}
		}
	} else {
		die "WTF?\n";
	}
}

for ($i=0; $i<1000; $i++) {
	print(($i/10)."\t".$out[$i]."\n");
}
