Some seismic design codes require to determine earthquake motions at ground surface, although we need sometimes a motion at engineering basement of ground. For example, it is necessary for optimal design of structures to introduce effects of soil-structure interaction and to input ground motions to engineering basement, which should satisfy the motions at surface required by a design code. This is simplified as a problem to find a possible input signal of a known non-linear-single-input and single-output (SISO) system with a given output signal. This study proposes a very simple algorithm to find a possible input signal for a known SISO system satisfying a given output signal from the system. The proposed algorithm searches for a possible input signal, whose system response is similar in shape to the given output signal, without any system identifications. Thus, it is not an inversion technique; it is an algorithm used to search for one of the possible input signals. Herein, the proposed algorithm is described and its performance is demonstrated using numerical examples.