#!/usr/bin/perl
#
# Compute values of position, velocity, acceleration, force
#   as functions of time
#   for an object in SHM
#
# MWR 4/17/2026

$debug = 1;

$mypi = 3.1415926;

$A = 2.0;
$omega = 3.5;
$phi = 0.0;

$start_t = 0.0;
$end_t = 30.0;
$d_t = 0.01;
for ($t = $start_t; $t <= $end_t; $t += $d_t) {

  $x = $A * cos($omega*$t + $phi);
  $v = - $omega * $A * sin($omega*$t + $phi);
  $a = - ($omega*$omega) * $A * cos($omega*$t + $phi);

  printf " %12.6f  %12.6f %12.6f %12.6f \n",
      $t, $x, $v, $a;
}



exit 0;