function PF_DPF_PRSS;
w=tf([0 0 1],[50 15 1]);
step(w);
z=c2d(w,4);
step(z);
sys=ss(z);
step(sys)