t1 = [0; 1; 3; 4; 5.5; 6]
tau1 = t1
A1 = construct(t1, tau1)

tau2 = [0.5; 2; 2.5; 3.5; 4.5; 5.75]
y2 = [1; 1.5; 1.25; 0; 0; 1.5];
A2 = construct(t1, tau2);
x2 = mcnorm(A2, y2);
plot(t1, x2)
plot(tau2,y2,'.r')




tau3 = 0.3 * (0:20)'
y3=[0 0.6 1.4 1.7 2.1 1.9 1.6 1.4 1.4 1 0.5 0.4 -0.2 -0.8 -0.5 0 0.4 1 1.6 1.7 1.2]'
A3 = construct(t1, tau3)
x3 = mcnorm(A3, y3);
plot(t1, x3)
plot(tau3,y3,'.r')
// QR