u := [2,1]: v:= [-3,4]: mix := c -> evalm(c*u+(1-c)*v); m0:=mix(0); m1:= mix(1/4); m2:= mix(3/4); m3 := mix(1);