s=NDSolve[{y'[x]==  -2.2067/1000000000000(y[x]^4-81*10^8),y[0]== 1200},y,{x,0,500}]

Plot the solution for different initial values:

Plot[Evaluate[y[x]/.%],{x,0,500}, Filling->Axis]