In this paper, we discuss tests of heteroscedasticity and/or autocorrelation in nonlinear models with AR(1) and symmetrical errors. The symmetrical errors distribution class includes all symmetrical continuous distributions, such as normal, Student-t, power exponential, logistic I and II, contaminated normal, so on. First, score test statistics and their adjustment forms of heteroscedasticity are derived. Then, the asymptotic properties, including asymptotic chi-square and approximate powers under local alternatives of the score tests, are studied. The properties of test statistics are investigated through Monte Carlo simulations. Finally, a real data set is used to illustrate our test methods.