library(strucchange)

set.seed(123)
datats <- as.ts(cumsum(rnorm(20,10,40)),start=1995)
plot(datats)
fs.ts <- Fstats(datats ~ 1)
summary(fs.ts)
breakpoints(fs.ts)
lines(breakpoints(fs.ts))