fid = fopen('output','w')
      a=linspace(1,100,100);
      fprintf(fid,'%g degC = %g degF\n',a,32+1.8*a);
fid=fopen('output','r')
     x=fscanf(fid,'%g degC = %g degF');
     size(x)
     x=reshape(x,(length(x)/2),2)