DSolve[u''=2u(3u+1),u]