To effectively represent photovoltaic (PV) modules while considering their dependency on changing environmental conditions, three novel mathematical and empirical formulations are proposed in this study to model PV curves with minimum effort and short timing. The three approaches rely on distinct mathematical techniques and definitions to formulate PV curves using function representations. We develop our models through fractional derivatives and stochastic white noise. The first empirical model is proposed using a fractional regression tool driven by the Liouville-Caputo fractional derivative and then implemented by the Mittag-Leffler function representation. Further, the fractional-order stochastic ordinary differential equation (ODE) tool is employed to generate two effective generic models. In this work, multiple commercial PV modules are modeled using the proposed fractional and stochastic formulations. Using the experimental data of the studied PV panels at different climatic conditions, we evaluate the proposed models’ accuracy using two effective statistical indices: the root mean squares error (RMSE) and the determination coefficient (R2). Finally, the proposed approaches are compared to several integer-order models in the literature where the proposed models’ precisely follow the real PV curves with a higher R2 and lower RMSE values at different irradiance levels lower than 800 w/m2, and module temperature levels higher than 50 °C.