We have proposed a new mathematical method, SEIHCRD-Model that is an extension of the SEIR-Model adding hospitalized and critical twocompartments. SEIHCRD model has seven compartments: susceptible (S), exposed (E), infected (I), hospitalized (H), critical (C), recovered (R), and deceased or death MESHD (D), collectively termed SEIHCRD. We have studied COVID- 19 cases of six countries, where the impact of this disease MESHD in the highest are Brazil, India, Italy, Spain, the United Kingdom, and the United States. SEIHCRD model is estimating COVID-19 spread and forecasting under uncertainties, constrained by various observed data in the present manuscript. We have first collected the data for a specific period, then fit the model for death MESHD cases, got the values of some parameters from it, and then estimate the basic reproduction number TRANS over time, which is nearly equal to real data, infection MESHD rate, and recovery rate of COVID-19. We also compute the case fatality rate over time of COVID-19 most affected countries. SEIHCRD model computes two types of Case fatality rate one is CFR daily and the second one is total CFR. We analyze the spread and endpoint of COVID-19 based on these estimates. SEIHCRD model is time-dependent hence we estimate the date and magnitude of peaks of corresponding to the number of exposed cases, infected cases, hospitalized cases, critical cases, and the number of deceased cases of COVID-19 over time. SEIHCRD model has incorporated the social distancing parameter, different age groups TRANS analysis, number of ICU beds, number of hospital beds, and estimation of how much hospital beds and ICU beds are required in near future.