We present a Maple implementation of the well known global approach to time series analysis and some further developments designed to improve the computational efficiency of the forecasting capabilities of the approach. This global approach can be summarized as being a reconstruction of the phase space, based on a time ordered series of data obtained from the system. After that, using the reconstructed vectors, a portion of this space is used to produce a mapping, a polynomial fitting, through a minimization procedure, that represents the system and can be employed to forecast further entries for the series. In the present implementation, we introduce a set of commands, tools, in order to perform all these tasks. For example, the command VecTS deals mainly with the reconstruction of the vector in the phase space. The command GfiTS deals with producing the minimization and the fitting. ForecasTS uses all these and produces the prediction of the next entries. For the non-standard algorithms, we here present two commands: IforecasTS and NiforecasTS that, respectively deal with the one-step and the N-step forecasting. Finally, we introduce two further tools to aid the forecasting. The commands GfiTS and AnalysTS, basically, perform an analysis of the behavior of each portion of a series regarding the settings used on the commands just mentioned above. Program summaryProgram title: TimeSCatalogue identifier: AERW_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AERW_v1_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 3001No. of bytes in distributed program, including test data, etc.: 95018Distribution format: tar.gzProgramming language: Maple 14.Computer: Any capable of running MapleOperating system: Any capable of running Maple. Tested on Windows ME, Windows XP, Windows 7.RAM: 128 MBClassification: 4.3, 4.9, 5Nature of problem:Time series analysis and improving forecast capability.Solution method:The method of solution is partially based on a result published in [1].Restrictions:If the time series that is being analyzed presents a great amount of noise or if the dynamical system behind the time series is of high dimensionality (Dim≫3), then the method may not work well.Unusual features:Our implementation can, in the cases where the dynamics behind the time series is given by a system of low dimensionality, greatly improve the forecast.Running time:This depends strongly on the command that is being used.