function [C,S]=trig_val(indx,eval_pts)
kk=size(indx);
if kk(1)~=1,
   error('First argument is not a row vector.');
end
t=eval_pts;
C=[ones(1,length(t));cos(indx'*t)];
S=sin(indx'*t);