Global Power Prediction Systems prototype version 2021 is presented with its system decomposition, scope, geographical/administrative/power grid decompositions, and similar. “Welcome”, “sign-up”, “log-in”, and “non-registered user main” web-interfaces are designed as draft on Quant UX. Map canvas is given as world political map with/without world power grid layers on QGIS 3.16.7-Hannover. Data input file is prepared based on several sources (1971-2018). It includes minimum and maximum values due to source value differences. 70/30 principle is applied for train/test splitting (training/testing sets: 1971-2003/2004-2018). 10 models are prepared on R version 4.1.1 with RStudio 2021.09.0+351. These are R::base(lm), R::base(glm), R::tidymodels::parsnip(engine("lm")), R::tidymodels::parsnip(engine("glmnet")) with lasso regularization, R::tidymodels::parsnip(engine("glmnet")) with ridge regularization, R::forecast(auto.arima) auto autoregressive integrated moving average (ARIMA), R::forecast(arima) ARIMA(1,1,2), and ARIMA(1,1,8). Electricity demand in kilowatt-hours at the World level zone for up to 500-years (2019-2519) prediction period with only 1-year interval is forecasted. The best model is the auto ARIMA (mean absolute percentage error MAPE and symmetric mean absolute percentage error SMAPE for minimum and maximum electricity consumption respectively 1,1652; 6,6471; 1,1622; 6,9043). Ex-post and ex-ante plots with 80%-95% confidence intervals are prepared in R::tidyverse::ggplot2. There are 3 alternative scripts (long, short, RStudio Cloud). Their respective runtimes are 41,45; 25,44; and 43,33 seconds. Ex-ante 500-year period (2019-2519) is indicative and informative.