Atrial fibrillation (AF) is the most frequent form of arrhythmia occurring in the industrialized world. Because of its complex nature, each identified form of AF requires specialized treatment. Thus, an in-depth understanding of the bases of these arrhythmias is essential for therapeutic development. A variety of experimental studies aimed at understanding the mechanisms of AF are performed using primary cultures of neonatal rat atrial cardiomyocytes (NRAMs). Previously, we have shown that the distinct advantage of NRAM cultures is that they allow standardized, systematic, robust re-entry induction in the presence of a constitutively-active acetylcholine-mediated K+ current (IKACh-c). Experimental studies dedicated to mechanistic explorations of AF, using these cultures, often use computer models for detailed electrophysiological investigations. However, currently, no mathematical model for NRAMs is available. Therefore, in the present study we propose the first model for the action potential (AP) of a NRAM with constitutively-active acetylcholine-mediated K+ current (IKACh-c). The descriptions of the ionic currents were based on patch-clamp data obtained from neonatal rats. Our monolayer model closely mimics the action potential duration (APD) restitution and conduction velocity (CV) restitution curves presented in our previous in vitro studies. In addition, the model reproduces the experimentally observed dynamics of spiral wave rotation, in the absence and in the presence of drug interventions, and in the presence of localized myofibroblast heterogeneities.