Functions returns variable values at each iteration

So I have this function where I'm calculating the values of the U matrix at each t value (0,1,2,etc.). The problem I'm having is that the function only returns the final value for U, and does not return t value. To illustrate what I'm trying to get, here's an example:

t   x y
0   2    1
1   2.6  0.9
2   3.458 0.891

etc.

function [U,t] = IVP(a,b,c)
syms x y h t;
h=c;
x=a;
y=b;    
for t=0:1:30
    f=[1.2*x-0.6*x*y;0.3*x*y-0.8*y];
    u=[x;y];
    U=u+h*f;
    x=U(1);
    y=U(2);
end
end

Answer 1

[Azzi Abdelmalek]

U=[]; 
for t=0:1:30
  f=[1.2*x-0.6*x*y;0.3*x*y-0.8*y];
  u=[x;y];
  w=u+h*f;
  U=[U w];
  x=w(1);
  y=w(2);
end
t=(0:1:38)'