function trace = Kalman_1D_driver(max_time)
%

A = 1;
B = 1;
C = 1;
R = 10;
Q = 0.01;
sigma = .1;
mu = 1;
u = 1;
x = 1;
trace = [x,x,mu,sigma];

for t = 0:max_time
    z = C*x + normal_1D(0,Q);
    [mu,sigma] = Kalman_1D(mu,sigma,u,z,A,B,R,C,Q);
    trace = [trace; x,z,mu,sigma];
    x = x + u;
end
