[t,Y]=ode45(@hw4_3_fun,[0 2],[0.5 0.5]);
figure(1)
plot(t,Y(:,1),'-',t,Y(:,2),'--')
figure(2)
plot(Y(:,1),Y(:,2))