#!/usr/bin/perl -w

use strict;

print <<EOF;
# harmonic bond

bond
N 10001

EOF

my $r=0.0;
for (my $i=1; $i < 10002; ++$i) {

  $r = 0.001*$i;
  printf "%d %.3f %.6f %.6f\n",
	$i, $r, 0.5*899.95759*($r-1.54)*($r-1.54),
    -899.95759*($r-1.54);
}