[t,x]=ode45(@sedacka,[0 5],[0 0]); plot(t,x(:,1))