I do agree with you. But to answer where is the probability of finding maximum, simply evaluating probability density wouldn't do the job. For 1s we may evaluate probability of finding between r and r+dr using volume of shell to get the answer. But say we have to do it for 2p. Then in this case, we cannot directly use density to compare, since even in 1s, using simply density gives points at nucleus as points of maximum probability, but there isn't any volume there (volume of shell becomes zero at nucleus), meaning there aren't any points there too. How should I do this for 2p orbital?