We present a flexible, open-source Python package for the accurate simulation of the z-propagation dynamics of ultrashort optical pulses in nonlinear waveguides, especially valid for few-cycle pulses and their interaction. The simulation approach is based on unidirectional propagation equations for the analytic signal. The provided software allows to account for dispersion, attenuation, four-wave mixing processes including, e.g., third-harmonic generation, and features various models for the Raman response. The propagation equations are solved on a periodic temporal domain. For z-propagation, a selection of pseudospectral methods is available. Propagation scenarios for a custom propagation constant and initial field pulses can either be specified in terms of a HDF5 based input file format or by direct implementation using a python script. We demonstrate the functionality for a test-case for which an exact solution is available, by reproducing exemplary results documented in the scientific literature, and a complex propagation scenario involving multiple pulses. Program summary Program Title: py-fmas CPC Library link to program files: https://doi.org/10.17632/7s2cv9kjfs.1 Developer's repository link: https://github.com/omelchert/py-fmas Code Ocean capsule: https://codeocean.com/capsule/8221780 Licensing provisions: MIT Programming language: Python3 Supplementary material: Reference manual, extended user guide, and usage examples are hosted on gitHub pages under https://omelchert.github.io/py-fmas. Nature of problem: Solves for the z-propagation dynamics of spectrally broad ultrashort optical pulses in single mode nonlinear waveguides in terms of propagation models for the analytic signal of the optical field [1-3]. The implemented models include, e.g., third-harmonic generation and the Raman effect. Solution method: The initial real-valued optical field is defined on a periodic one-dimensional temporal grid and converted to the complex-valued analytic signal. z-stepping is performed via spectral methods. The software implements a selection of algorithms with fixed or adaptive stepsize, commonly used in nonlinear optics for solving nonlinear Schrodinger type equations. Additional comments including restrictions and unusual features: The range of applicability of the provided software is equivalent to that of the forward Maxwell equation [4]. For reasonably chosen initial conditions, it can be used beyond the unidirectional approximation as a bidirectional model for a complex field, allowing to describe forward and backward waves coupled through nonlinear interactions [1]. The software implements various models for the Raman response and allows to calculate spectrograms, detailing the time-frequency composition of the analytic signal. Additionally, a convenience class for analyzing propagation constants is provided. (C) 2021 Elsevier B.V. All rights reserved.