The damping work fine by following damping ratio rule
ζ=C/Cc, Cc = 2√km
program dashpot
implicit none
integer,parameter :: N=10001,IN=300
integer :: i
real :: pi,dt,k,m,c,U,d,wb,am
real,dimension (N)::T,V,X,f
!initial condition= position,velocity,time, pi
pi = 3.14159265359
dt =...