parameter( maxnt = 100000 )
	real tt(maxnt),u(maxnt)
	character *70 fn,outfn
 
    	write(*,'(a)') 'please input input  file name'
	read(*,'(a)') fn
	write(*,'(a)') 'please input output file name'
        read(*,'(a)') outfn
	open(10,file = outfn)


             	
         call rsac1(fn,u,npts,beg,delt,maxnt,nerr)
         do i = 1, npts
	 write(10,*) beg+(i-1)*delt,u(i)	   
	 enddo
	
	end