The modeling procedure starts with presenting nonlinear dynamics of cancer cells and CAFs using ordinary differential equations based on TGFβ, CXCL12, and LIF signaling pathways. The gene discussed is CXCL12; the disease is cancer.