Model joint variation equations