A new method for computer analysis of nystagmic eye movements in vestibulo-ocular (VOR) and optokinetic (OKN) reflexes is developed. A fuzzy set theoretical approach is used to construct the slow cumulative eye position (SCEP) curve by eliminating fast components (saccades) from eye movement signal. These procedures are able to perform automatically some pattern recognition tasks traditionally used--in classical interactive programs--when human operators distinguish between fast- and slow-phases of eye movements. The structure of the algorithm is as follows. 1) A fuzzy clustering of slow- and fast-phases is made. An iterative method is used to refine the membership function of slow-phases, step by step, until a sufficiently discriminating membership function is obtained. 2) Saccades are detected and removed from the eye position signal. 3) SCEP is then built by interpolating between slow phases. 4) A weighted least-squares curve fitting is made. Weighting coefficients are obtained from the last membership function resulting from iterations in step 1). This curve fitting is referenced to the SCEP and the parameters of VOR and OKN are calculated using this last curve. This approach permits an analysis of nystagmic eye movements with high reliability even when the data are of modest quality. The definitive innovative feature of the program is that it allows entirely automatic analysis without participation of a human operator. The main algorithm being independent of the shape of stimulus, the program can be generalized to fit any type of simulation.