This study uses the optimal auxiliary function method to approximate solutions for fractional-order non-linear partial differential equations, utilizing Riemann–Liouville’s fractional integral and the Caputo derivative. This approach eliminates the need for assumptions about parameter magnitudes, offering a significant advantage. We validate our approach using the time-fractional Cahn–Hilliard, fractional Burgers–Poisson, and Benjamin–Bona–Mahony–Burger equations. Comparative testing shows that our method outperforms new iterative, homotopy perturbation, homotopy analysis, and residual power series methods. These examples highlight our method’s effectiveness in obtaining precise solutions for non-linear fractional differential equations, showcasing its superiority in accuracy and consistency. We underscore its potential for revealing elusive exact solutions by demonstrating success across various examples. Our methodology advances fractional differential equation research and equips practitioners with a tool for solving non-linear equations. A key feature is its ability to avoid parameter assumptions, enhancing its applicability to a broader range of problems and expanding the scope of problems addressable using fractional calculus techniques.

In an increasingly data‐driven world, statistical literacy is a necessity yet statistical learning is often inhibited by statistics anxiety. Using the Auzmendi Scale to Measure Attitude toward Statistics (ASMAS), this study examines how statistics anxiety in university students is related to other dimensions of their attitudes toward statistics and how statistics anxiety and other dimensions change following introductory statistics instruction. Based on data collected from Spain, Canada, and Australia, this study finds that anxiety is negatively related to security–confidence, pleasantness, and motivation. The structure of these relationships is consistent across countries and disciplines and remains in place after statistics instruction. Further, by the end of an introductory statistics course, students report higher security–confidence and pleasantness but lower anxiety. Results thus suggest where efforts to improve students' experience with statistics might need to be directed, and the paper concludes with a discussion of the implications of these results for statistics instruction.

In this article, we investigate the utilization of Riemann–Liouville’s fractional integral and the Caputo derivative in the application of the Optimal Auxiliary Function Method (OAFM). The extended OAFM is employed to analyze fractional non-linear coupled ITO systems and non-linear KDV systems, which feature equations of a fractional order in time. We compare the results obtained for the ITO system with those derived from the Homotopy Perturbation Method (HPM) and the New Iterative Method (NIM), and for the KDV system with the Laplace Adomian Decomposition Method (LADM). OAFM demonstrates remarkable convergence with a single iteration, rendering it highly effective. In contrast to other existing analytical approaches, OAFM emerges as a dependable and efficient methodology, delivering high-precision solutions for intricate problems while saving both computational resources and time. Our results indicate superior accuracy with OAFM in comparison to HPM, NIM, and LADM. Additionally, we enhance the accuracy of OAFM through the introduction of supplementary auxiliary functions.