chart := (head,vals) -> 
     MATRIX([head, seq(evalf(eval(subs(head[1]=vals[i],head))),
                       i=1..nops(vals))]):
