clc
syms s
h=s/s+1; 
u=tf('h')