def IndustWater(time):
pop = (Gaussian(Fit.PopulationFit(time),pPara)/Gaussian(Fit.PopulationFit(years[-1]),pPara))
pcg = (Gaussian(Fit.PCGDPFit(time),gPara)/Gaussian(Fit.PCGDPFit(years[-1]),gPara))
ele = (Gaussian(Fit.ElectricityFit(time),ePara)/Gaussian(Fit.ElectricityFit(years[-1]),ePara))
spr = (Gaussian(Fit.SteelProductFit(time),sPara)/Gaussian(Fit.SteelProductFit(years[-1]),sPara))
return Data.WaterUseIndustry[years[-1]]*(pop/4+pcg/4+ele/4+spr/4)
plt.show()
plt.close()
plt.figure()
plt.xlabel("Population /10k people")
plt.ylabel("Industry Water Usage/ 100m m^3")
plt.scatter([Data.Population[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.PopulationFit(x),Gaussian(Fit.PopulationFit(x),pPara),'g-',label="Fit")
plt.legend()
plt.show()
plt.close()
plt.figure()
plt.xlabel("PCGDP /CNY")
plt.ylabel("Industry Water Usage/ 100m m^3")
plt.scatter([Data.PCGDP[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.PCGDPFit(x),Gaussian(Fit.PCGDPFit(x),gPara),'g-',label="Fit")
plt.legend()
plt.show()
plt.close()
plt.figure()
plt.xlabel("Electricity /100M kWh")
plt.ylabel("Industry Water Usage/ 100m m^3")
plt.scatter([Data.Electricity[t] for t in years],[Data.WaterUseIndustry[t] for t in years],marker="^",s=50,label='Raw Data')
plt.plot(Fit.ElectricityFit(x),Gaussian(Fit.ElectricityFit(x),ePara),'g-',label="Fit")
plt.legend()
plt.show()
plt.close()
