TGFB1 and cancer: 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.