We consider the bilinear Fourier integral operatorS(f, g)(x) =ZRdZRdei1(x,)ei2(x,)(x, , ) ˆ f()ˆg()d d,on modulation spaces. Our aim is to indicate this operator is well defined onS(Rd) and shall show the relationship between the bilinear operator and BFIO onmodulation spaces.