% Approximate               f(z)     =  log10((1+z)/(1-z)),
% by                        phi(j,z) = z^j, j=0:2,
% on the real unit interval.
%
    clear
%   format long

    nomf      =    'log10x';


    rr=(sqrt(10)-1)/(sqrt(10)+1); a=-rr ;  b=rr ;
    n=4 ;      %degree
    np=200 ;    %   grid
    iter=20;
% coefficients de x^n , x^(n-1) , ... ,1 :
    coeff=remezp(nomf,a,b,n,iter,np);

