%
% Entrada-salida a y desde archivos
%

f1 = fopen('datos_ejemplo3.txt', 'w');

for i = 1:100
  fprintf(f1, '%6d\n', i);
end

status = fclose(f1);

f2 = fopen('datos_ejemplo3.txt', 'r');
f1 = fopen('datos_ejemplo3b.txt', 'w');

clear v;

npoints = 0;
while (~ feof(f2))
  [x, ndat] = fscanf(f2, '%d', 1);
  if (ndat > 0)
    npoints = npoints + 1;
    x2 = x * x;
    xs = sqrt(x);
    fprintf(f1,'%16.8g %16.8g %16.8g\n', x, x2, xs);
  end
end

fclose(f2); fclose(f1);

fprintf('El número de datos procesados es %d\n', npoints);

