//trace
t=[0;1;2]
y=[1;3;7]
T=linspace(t(1),t($),1000)';
A=constrpoly(t)
x=mcnorm(A,y)
A=constrpoly(T)
Y=A*x
plot(t,y,'.r')
plot(T,Y)


t=linspace(0,2,9)'
y =[1, 1.7, 1.95, 1.8, 3., 3.6, 4.45, 5.9, 6.6]'
T=linspace(t(1),t($),1000)'
A=constrpoly(t)
x=mcnorm(A,y)
A=constrpoly(T)
Y=A*x
plot(t,y,'.r')
plot(T,Y)