The malware spreading in Wireless Sensor Network (WSN) has lately attracted the attention of many researchers as a hot problem in nonlinear systems. WSN is a collection of sensor nodes that communicate with each other wirelessly. These nodes are linked in a decentralised and distributed structure, allowing for efficient data collection and communication. Due to their decentralised architecture and limited resources, WSN is vulnerable to security risks, including malware attacks. Malware can attack sensor nodes, causing them to malfunction and consume more energy. These attacks can spread from one infected node to others in the network, making it essential to protect WSN against malware attacks. In this paper, we focus on the analysis of a novel fractional epidemiology model, specifically the fractional order SEIVR epidemic model in the sense of Caputo's fractional derivative of order 0< α≤1 with the goal of examining the efficacy of vaccination strategies and the heterogeneity of a scale-free network on epidemic spreading. First, using the next-generation technique and obtain the basic reproduction number of the proposed epidemic model, which is essential for determining both the locally asymptotically stable equilibrium point of the wormfree system and the unique existence of the endemic equilibrium point. To numerically solve the model, the Adam-Bashforth-Moulton predictor-corrector (ABM) method is applied. The fractional calculus enables us to deal directly with the ''memory effect'' of numerous phenomena, taking into account the system's dependence on previous stages. This method provides the results of a complex system. Additionally, research demonstrates that vaccine treatments are quite effective at preventing the spread of malware. The outcome of the study reveals that the applied ABM predictor-corrector method is computationally strong and effective to analyse fractional order dynamical systems in the SEIVR epidemic model for malware propagation in WSN. The results show that the order of the fractional derivative has a significant effect on the dynamic process. Also, from the result, it is obvious that the memory effect is zero for α = 1. When the fractional order α is decreased from 1, the memory effect appears, and its dynamics vary according to the time. This memory effect points out the difference between derivatives of fractional and integer orders. The theorems and their proofs are presented to validate the validity of the proposed model. To validate the proposed model, extensive theoretical study and computational analysis have also been applied.The associate editor coordinating the review of this manuscript and approving it for publication was Stefano Scanzio .