Variational iterative method for system fractional DE with Katugampola derivative?

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: