#!/usr/bin/perl
#
$k = 120.0;
$L = 2.0;
$theta_deg = 20.0;
$theta_rad = ($theta_deg)*(3.14159/180.0);
$m = 1.5;
$x0 = 0.20;
$muk = 0.07;
$g = 9.8;

for ($x3 = -0.30; $x3 < 0.20; $x3 += 0.001) {

  $a = 0.5*$k*$x3*$x3;
  $b = $x3*$m*$g*sin($theta_rad);
  $c = $x3*0.9669;
  $d = 1.202;

  $z = $a + $b + $c;
  $z2 = $z - $d;

  printf " %9.4f  %9.4f   %9.4f \n",
      $x3, $z, $z2;

}

exit 0;