Abstract
This paper, presents a new perspective on the numerical solution of a fuzzy linear system of differential equations, where initial values
and constant coefficients are fuzzy numbers. To do this, the matrix of
coefficients was first decomposed into two matrices, then the variational
iteration method was presented and applied to them. Afterwards, the convergence of the above mentioned method was proved. Finally, some examples are presented which confirm the applicability, accuracy, and efficiency of the method
