//mode(-1);

getf('Unstab.sci');
getf('funcs.sci');

tmax=1;
ok=1;
x_message(['Résolution de l''équation différentielle', 'y'' = y, sur (0,1)', ...  
    'par une méthode instable'])
while ok do
  h=evstr(x_dialog('Choisir une valeur de h', '0.1'));  
  xbasc();
  ok = (h~=[]);
  if ok then
    [tvals, yvals]=Unstab(f2,h,tmax);
  end
  end