سلام
کاری نداره
clear all;
close all;
clc;
num=1;
den=poly([-3 -3 -3 -3 -3 -3]);
G=tf(num,den);
step(G)
t=-10:0.01:10;
u=zeros(size(t));
u(t>0)=1;
[y, t]=lsim(G,u,t);
hold on
plot(t,y,'r')
legend('ba dastore step','bedone dastore step')
%%% --- mesal
% -- y ra be eza x=5 hesab konid
[z,v]=find(t==5);
y(v)