real sx(10000), sy(10000)
        open(1,file='system',status='old')

	i=1
1	read(1,*,end=99) sx(i),sy(i)
        i=i+1
	goto 1

99	ns=i-1

	open(1,file='filters/f300w',status='old')
101	read(1,*,end=199) fx, fy
	call linint(sx,sy,ns,fx,r)
	print *, fx, fy*r   !, fy, r
	goto 101

199	continue

	stop
	end