An　ejection　dynamics　mathematical　model　of　human　left　ventricle　(LV)　based　on　physiological　data　of　human　heart　is　proposed　in　this　study.　The　mathematical　equations　were　expressed　in　terms　of　vorticity-stream　function　equations　in　a　prolate　spheroidal　coordinate　system.　These　equations　combined　with　specified　boundary　conditions　were　numerically　solved　by　using　an　alternating-direction-implicit　(ADI)　algorithm　with　second　order　accuracy.　The　unsteady　aspects　of　the　ejection　process　were　subsequently　introduced　into　the　numerical　simulation.　The　numerical　results　have　shown　that　the　present　ellipsoidal　model　could　be　available　to　simulate　the　ejection　process　of　the　human　LV.　Such　a　model　combined　with　cardiac　muscle　mechanics　could　be　studied　further　to　determine　altered　left　ventricular　function　in　cardiac　diseases.