function?[accveldsp]=TransResp5（kkccmmfdbcdofntdtq0dq0）
%--------------------------------------------------------------------------
%??Purpose:
%?????The?function?subroutine?TransResp5.m?calculates?transient?response?of
%?????a?structural?system?using?Newmark?integration?scheme.
%??Synopsis:
%?????[accveldsp]=TransResp5（kkccmmfdbcdofntdtq0dq0）
%??Variable?Description:
%?????Input?parameters
%???????kk?cc?mm?-?stiffness?damping?and?mass?matrices
%???????fd?-?Input?or?forcing?influence?matrix
%???????bcdof?-?Boundary?condition?dofs?vector
%???????nt?-?Number?of?time?steps
%???????dt?-?Time?step?size
%???????q0?dq0?-?Initial?condition?vectors
%?????Output?parameters
%???????acc?-?Acceleration?response
%???????vel?-?Velocity?response
%???????dsp?-?Displacement?response
%--------------------------------------------------------------------------
%??（1）?initial?condition
%--------------------------------------------------------------------------
[sdofn2]=size（kk）;
?
dsp=zeros（sdofnt）;?????????????????????????????????????????%?displacement?matrix
vel=zeros（sdofnt）;?????????????????????????????????????????????%?velocity?matrix
acc=zeros（sdofnt）;??????????????????????????????????????????%?acceleration?matrix
?
dsp（:1）=q0;???????????????????????????????????????????????%?initial?displacement
vel（:1）=dq0;??????????????????????????????????????????????????%?initial?velocity

alpha=0.5;?beta=0.5;????????????????????????????????????????%?select?the?parameters
%--------------------------------------------------------------------------
%??（2）?Newmark?integration?scheme?for?time?integration
%--------------------------------------------------------------------------
acc（:1）=inv（mm）*（fd（:1）