function ofc_graphics4mesh_install(varargin)
  if strversion2num(version)<strversion2num('4.2.0')
    error('Using Octave version %s. Version must be >= 4.2.0. ',version)
  end
  warning('off','Octave:shadowed-function');more off
  v='0.0.3';
  v_fc_mesh='0.0.1';
  p = inputParser; 
  p.addParamValue('version', v, @ischar );
  p.addParamValue('dir', './fc-graphics4mesh-full', @ischar );
  p.addParamValue('verbose', 1, @(x) ismember(x,0:2) ); % level of verbosity 
  p.addParamValue('fc_mesh',v_fc_mesh , @ischar ); % fc_mesh version
  p.parse(varargin{:});
  version=p.Results.version;verbose=p.Results.verbose;
  vprintf(verbose,1,'Parts of the GNU Octave <fc-graphics4mesh> package.\nCopyright (C) 2017 Francois Cuvelier <cuvelier@math.univ-paris13.fr>\n\n');
  currentdir=pwd;
  vprintf(verbose,1,'**************************************\n');
  vprintf(verbose,1,'Downloading and installing the package\n');
  vprintf(verbose,1,'    <fc-mesh>[%s]\n',p.Results.fc_mesh)
  vprintf(verbose,1,'*************** start ****************\n');
  fc_mesh_dir=get_install('mesh',p.Results.fc_mesh,p.Results.dir,'graphics',true,'verbose',verbose);
  vprintf(verbose,1,'**************** end *****************\n');
  vprintf(verbose,1,'*** Using instructions \n');
  vprintf(verbose,1,'   To use the <fc-graphics4mesh> toolbox:\n');
  vprintf(verbose,1,'   addpath(''%s'')\n',fc_mesh_dir);
  vprintf(verbose,1,'   fc_mesh.init()\n');
  logfile=[currentdir,filesep,'ofc_graphics4mesh_set.m'];
  vprintf(verbose,1,'\n   See %s\n',logfile);
  fid=fopen(logfile,'w');
  fprintf(fid,'warning(''off'',''Octave:shadowed-function'');more off\n');
  fprintf(fid,'addpath(''%s'')\n',fc_mesh_dir);
  fprintf(fid,'fc_mesh.init()\n');
  fclose(fid);
end

function fulldir=get_install(name,ver,dir,varargin)
  mfile=sprintf('ofc_%s_install',name);
  urlfile=sprintf('http://www.math.univ-paris13.fr/~cuvelier/software/codes/Octave/fc-%s/%s/%s.m',name,ver,mfile);
  try
    urlwrite(urlfile,[mfile,'.m']);
  catch
    error('Unable to get :\n  -> %s\n',urlfile);
  end
  eval([mfile,'(''dir'',dir,varargin{:})']) % run installation 
  fulldir=fullfile(pwd,sprintf('%s%sfc_%s-%s',dir,filesep,name,ver));
end

function ver=strversion2num(strv)
% convert version string 'x.y.z' in x*10^6+y*10^3+z
S=strsplit(strv,'.');
assert(length(S)==3);
z=str2num(S{3});assert(z<10^3);
ver=z;
y=str2num(S{2});assert(y<10^3);
ver=ver+10^3*y;
x=str2num(S{1});assert(x<10^3);
ver=ver+10^6*x;
end

function vprintf(verbose,level,varargin)
  if verbose>=level, fprintf(varargin{:});end
end
