%%%%central?difference?methrod
clc
clear?all
close?all
%%%%intial?values?of?the?system
m=0.2533;
c=0.1592;
k=10;
%%%%the?expression?of?excitation
syms?t
p=10*sin（pi*t/0.6）;
%%%%intial?caculation
t0=0.6;?????%the?end?of?the?excitation
total=1.0;???%totla?caculation?time
dt=0.1;?????%time?step
n=total/dt;
n1=t0/dt;
u1=0;
v1=0;
a1=（subs（p0）-c*v1-k*u1）/m;
uu=u1-dt*v1+0.5*dt^2*a1;
keq=m/dt^2+c/2/dt;