Stuff

function [T,X]=disparo(a,b,h,Y0,tol,factor)
T=a:h:b;
N=length(T);
options=odeset('RelTol',tol);
[T,X1]=ode45(@sist,T,[Y0(1) 0],options,factor);
Z1=X1(:,1);
[T,X2]=ode45(@sist,T,[0 1],options,factor);
Z2=X2(:,1);
X=Z1+((Y0(2)-Z1(N))/Z2(N))*Z2;
plot(T,X)
shg
function y=sist(t,x,factor)
y=[x(2),2*factor*(0.5)*x(2) + [(0.5).^2]*x(1) + 10*(factor*cos(0.5+0.5*(t+1))-sin(0.5+0.5*(t+1)))]';