[t,y]=ode45('ejem1',[0,10],5);
plot(t,y)
xlabel('t')
ylabel('y')