echo off;
global kp tp omega a;
a=2;kp=2;omega=0.5;tp=
[t,y]=ode23('ord1',[t0 tf], y0 );

for i=1:length(t);
    uc(i)=a*sin(omega*t(i));
end;
clf;
plot(t,y,'-r', t, uc, '-b');
grid;
title('rsp sist ord 1 ');