I want to know how to write a Mathematica code for variational iterative method for system of fractional differential equations with caputo katugampola derivative.
Start with these documentation pages and you can put them together: https://reference.wolframcloud.com/language/tutorial/FractionalCalculus.html https://reference.wolfram.com/language/guide/FractionalCalculus.html