function c=gencheb(f)
     [n,cols]=size(f);
     fac=2.0/sym(n);
     c=sym(zeros(n,1));
     dd=digits;
     for j=1:n

sum=sym(0);
j
for k=1:n

sum=vpa(sum+f(k)*cos(sym(pi)*(sym(j)-1)*(sym(k)-sym(0.5))/sym(n)),dd);
end
c(j)=fac*sum;
end