function zarray = myEulerSolver(functionHandle, tarray, z0)
%Solves a system numerically using Euler's Method
%Inputs:
%   functionHandle      function handle of myrhs
%   tarray              array of time values
%   z0                  Initial conditions
%Output:
%   zarray              state of system over time

%Create zarray
n = length(tarray);
zarray = zeros(n,length(z0));

zarray(1,:) = z0(:);

for i = 1:n-1
    dt = tarray(i+1) - tarray(i);
    t = tarray(i);
    znow = zarray(i,:);
    
    zdot = functionHandle(t,znow);
    znext = znow + zdot.'*dt;
    zarray(i+1,:) = znext;
end
             
end