Introduction:　Hemodynamic　diagnosis　indexes　(HDIs)　can　comprehensively　evaluate　the　health　status　of　the　cardiovascular　system　(CVS),　particularly　for　people　older　than　50　years　and　prone　to　cardiovascular　disease　(CVDs).　However,　the　accuracy　of　non-invasive　detection　remains　unsatisfactory.　We　propose　a　non-invasive　HDIs　model　based　on　the　non-linear　pulse　wave　theory　(NonPWT)　applied　to　four　limbs.Methods:　This　algorithm　establishes　mathematical　models,　including　pulse　wave　velocity　and　pressure　information　of　the　brachial　and　ankle　arteries,　pressure　gradient,　and　blood　flow.　Blood　flow　is　key　to　calculating　HDIs.　Herein,　we　derive　blood　flow　equation　for　different　times　of　the　cardiac　cycle　considering　the　four　different　distributions　of　blood　pressure　and　pulse　wave　of　four　limbs,　then　obtain　the　average　blood　flow　in　a　cardiac　cycle,　and　finally　calculate　the　HDIs.Results:　The　results　of　the　blood　flow　calculations　reveal　that　the　average　blood　flow　in　the　upper　extremity　arteries　is　10.78　ml/s　(clinically:　2.5-12.67　ml/s),　and　the　blood　flow　in　the　lower　extremity　arteries　is　higher　than　that　in　the　upper　extremity.　To　verify　model　accuracy,　the　consistency　between　the　clinical　and　calculated　values　is　verified　with　no　statistically　significant　differences　(p　＜　0.05).　Model　IV　or　higher-order　fitting　is　the　closest.　To　verify　the　model　generalizability,　considering　the　risk　factors　of　cardiovascular　diseases,　the　HDIs　are　recalculated　using　model　IV,　and　thus,　consistency　is　verified　(p　＜　0.05　and　Bland-Altman　plot).Conclusion:　We　conclude　our　proposed　algorithmic　model　based　on　NonPWT　can　facilitate　the　non-invasive　hemodynamic　diagnosis　with　simpler　operational　procedures　and　reduced　medical　costs.