%solves equations described in hw4_1_fun.m, for time from 1 to 20, with
%initial condition Y(0)=3
[t,Y]=ode45(@hw4_1_fun,[0 20],[3]);
figure(1)
plot(t,Y,'-')