function ofc_simesh_install(varargin)
  if strversion2num(version)<strversion2num('4.2.0')
    error('Using Octave version %s. Version must be >= 4.2.0. ',version)
  end
  v='0.2.0';
  v_fc_tools='0.0.18';
  v_fc_hypermesh='0.0.5';
  v_fc_oogmsh='0.0.15';
  v_fc_graphics4mesh='0.0.1';
  v_fc_siplt='0.0.1';
  p = inputParser; 
  p.addParamValue('version', v, @ischar ); 
  p.addParamValue('dir', './fc-simesh-full', @ischar );
  p.addParamValue('graphics', true, @islogical );
  p.addParamValue('fc_tools',v_fc_tools , @ischar ); % fc_tools version
  p.addParamValue('fc_oogmsh',v_fc_oogmsh , @ischar ); % fc_oogmsh version
  p.addParamValue('fc_hypermesh',v_fc_hypermesh , @ischar ); % fc_hypermesh version
  p.addParamValue('fc_graphics4mesh',v_fc_graphics4mesh , @ischar ); % fc_graphics4mesh version
  p.addParamValue('fc_siplt',v_fc_siplt , @ischar ); % fc_siplt version
  %p.addParamValue('gmsh_bin','/usr/bin/gmsh',@ischar);
  p.parse(varargin{:});
  version=p.Results.version;
  fprintf('Parts of the GNU Octave <fc-simesh> package.\nCopyright (C) 2016-2017 Francois Cuvelier <cuvelier@math.univ-paris13.fr>\n\n');

  fprintf('1- Downloading and extracting the packages\n');
  fprintf('   -> <fc-tools>[%s] ...',p.Results.fc_tools)
  fc_tools_dir=get_inst_fc_gen('tools',p.Results.fc_tools,p.Results.dir);
  fprintf(' OK\n');
  fprintf('   -> <fc-hypermesh>[%s] ...',p.Results.fc_hypermesh)
  fc_hypermesh_dir=get_inst_fc_gen('hypermesh',p.Results.fc_hypermesh,p.Results.dir);
  fprintf(' OK\n');
  fprintf('   -> <fc-oogmsh>[%s] ...',p.Results.fc_oogmsh)
  fc_oogmsh_dir=get_inst_fc_gen('oogmsh',p.Results.fc_oogmsh,p.Results.dir);
  fprintf(' OK\n');
  fprintf('   -> <fc-simesh>[%s] ...',p.Results.version)
  fc_simesh_dir=get_inst_fc_gen('simesh',p.Results.version,p.Results.dir);
  fprintf(' OK\n');
  if p.Results.graphics
    fprintf('   -> <fc-graphics4mesh>[%s] ...',p.Results.fc_graphics4mesh)
    fc_graphics4mesh_dir=get_inst_fc_gen('graphics4mesh',p.Results.fc_graphics4mesh,p.Results.dir);
    fprintf(' OK\n');
    fprintf('   -> <fc-siplt>[%s] ...',p.Results.fc_siplt)
    fc_siplt_dir=get_inst_fc_gen('siplt',p.Results.fc_siplt,p.Results.dir);
    fprintf(' OK\n');
  end
  
  fprintf('2- Setting the packages\n');
  currentdir=pwd;
  fprintf('2-a) Setting the <fc-hypermesh> package\n');
  cd(fc_hypermesh_dir)
  fc_hypermesh.configure('fc_tools_dir',fc_tools_dir);
  fprintf('2-b) Setting the <fc-oogmsh> package\n');
  cd(fc_oogmsh_dir)
  fc_oogmsh.configure('fc_tools_dir',fc_tools_dir);%,'gmsh_bin',p.Results.gmsh_bin);
  cd(fc_simesh_dir)
  fprintf('2-c) Setting the <fc-simesh> package without graphics\n');
  fc_simesh.configure('oogmsh_dir',fc_oogmsh_dir,'hypermesh_dir',fc_hypermesh_dir);
  
  if p.Results.graphics
    fprintf('2-d) Setting the <fc-graphics4mesh> toolbox\n');
    cd(fc_graphics4mesh_dir)
    fc_graphics4mesh.configure('fc_tools_dir',fc_tools_dir);
    fprintf('2-e) Setting the <fc-siplt> toolbox\n');
    cd(fc_siplt_dir)
    fc_siplt.configure('fc_graphics4mesh_dir',fc_graphics4mesh_dir,'fc_simesh_dir',fc_simesh_dir)
    fprintf('2-f) Setting the <fc-simesh> toolbox with graphics\n');
    cd(fc_simesh_dir)
    fc_simesh.configure('oogmsh_dir',fc_oogmsh_dir,'hypermesh_dir',fc_hypermesh_dir, ...
                        'siplt_dir',fc_siplt_dir);
  end

  cd(currentdir)
  fprintf('3- Using instructions \n');
  fprintf('   To use the <fc-simesh> package:\n');
  fprintf('   addpath(''%s'')\n',fc_simesh_dir);
  fprintf('   fc_simesh.init()\n');
  logfile=[currentdir,filesep,'ofc_simesh_set.m'];
  fprintf('\n   See %s\n',logfile);
  fid=fopen(logfile,'w');
  fprintf(fid,'addpath(''%s'')\n',fc_simesh_dir);
  fprintf(fid,'fc_simesh.init()\n');
  fclose(fid);
end

function fulldir=get_inst_fc_gen(name,ver,dir)
  file=sprintf('ofc-%s-%s.tar.gz',name,ver);
  urlfile=sprintf('http://www.math.univ-paris13.fr/~cuvelier/software/codes/Octave/fc-%s/%s/fc_%s-%s.tar.gz',name,ver,name,ver);
  urlwrite(urlfile,file);
  try
    untar(file,dir)
  catch
    error('Unable to get :\n  -> %s\n',urlfile);
  end
  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