A method for automated modeling of plasma shape control response in tokamaks