predict model(the global pollution of ocean)_ocean_ICM2020年E题代码_

//this project is a predict model for the global sea's pollution #include<iostream> #include<math.h> using namespace std; double G(double x){ double y=0; y=0.548*x+0.031*x*x+0.001*x*x*x-0.4+2;//G(0)=2 (Wt) return y; }; double g(int x, double b){ double y=0; y=G(x)*b; return y; } double F(int t, int k, double a, double b){ double x, count=0; double c=1;//Discard rate 汕 x=t-1950; for(int i=0;i<=k-1;i++) { count+=g(i,b)*c; cout<<"year:"<<i+1950<<'\n'<<g(i,b)*c*a<<endl; } if(x>=k) { for(int j=k;j<=x;j++) { if(j>=38)//start recycling { if(j==38) { c-=0.006; cout<<"c="<<c<<endl; count=count+(g(j,b)+g(j,1-b))*c; cout<<"j="<<j<<endl; cout<<"year:"<<j+1950<<'\n'<<(g(j,b)*c+(1-b)*c*G(j-k))*a<<endl; } else { c-=0.007; cout<<"c="<<c<<endl; count=count+(g(j,b)*c+(1-b)*c*G(j-k)); cout<<"j="<<j<<endl; cout<<"year:"<<j+1950<<'\n'<<(g(j,b)*c+(1-b)*c*G(j-k))*a<<endl; } } else if(j>=31 && j<38)//Start incineration { if(j==31) { c-=0.017; cout<<"c="<<c<<endl; count+=g(j,b)*c; cout<<"j="<<j<<endl; cout<<"year:"<<j+1950<<'\n'<<g(j,b)*c*a<<endl; } else { c-=0.007; cout<<"c="<<c<<endl; count+=g(j,b)*c; cout<<"j="<<j<<endl; cout<<"year:"<<j+1950<<'\n'<<g(j,b)*c*a<<endl; } } else if(j<31)//all discard { cout<<"c="<<c<<endl; count+=g(j,b)*c; cout<<"j="<<j<<endl; cout<<"year:"<<j+1950<<'\n'<<g(j,b)*c*a<<endl; } } } return count; }; void main() { int t; int k; double a,b; cout<<"please input the final year:"<<endl; cin>>t; cout<<"please input the Probability of going into the sea:"<<endl;//(污/汕) cin>>a; cout<<"please input the Scrap rate:"<<endl;//汐 cin>>b; cout<<"please input the Average use time of plastics:"<<endl; cin>>k; cout<<"The final weight of plastic waste:"<<F(t,k,a,b)<<endl; }