t= 0:0.0001:1/180;
kp=10;
kd=1.3;
eqn ='0.00072*D3x+0.07206*D2x+(0.7274*kd+0.44244)*Dx+0.7274*kp*x=(0.7274*kd*Dt+0.7274*kp*t)*60*t';
inits='x(0)=0,Dx(0)=60,D2x(0)=0,t(0)=0';
x=dsolve(eqn,inits,'t');
z=eval(vectorize(x));
plot(t,z),grid on