xs=[0,1/3,2/3,1];
phi = @(x) x.*(x-2/3).*(x-1)/( (1/3)*(1/3-2/3)*(1/3-1) );
x=linspace(0,1,100);
plot(x,phi(x),xs,phi(xs),'*','LineWidth',2);
diary