function?[deltaudeltadiffudeltaddiffudeltaPdeltadiffPdiffuddiffuu]=newmark（MCKu2diffu2timedtP2）?
clear?all
M=input（‘輸入質量M?：‘）;
C=input（‘輸入阻尼C：‘）;
K=input（‘輸入剛度K：‘）;
P=input（‘輸入荷載P：‘）;
u2=input（‘輸入初始位移u2：?‘）;
diffu2=input（‘輸入初始速度diffu2：?‘）;
time=input（‘輸入模擬時間time：‘）;
dt=input（‘輸入時間步長dt?：‘）;
n=time/dt;
ddiffu2=（P（:1）-C*diffu2-K*u2）/M;
Ke=K+（3*C）/dt+（6*M）/（dt^2）;
a=（6*M）/dt+3*C;
b=3*M+（dt*C）/2;
for?i=2:1:n+2;
????if?i==2;
????????u（i）=u2;
????