If you have to do it the "hard way"; it's not even hard. In fact, you don't have to explicitly find any antiderivatives:
Let $c = \frac{1}{\sigma\sqrt{2\pi}}$ .
Then, using the linearity of integration, and noting all the integrals converge:
$\eqalign{ \Bbb E(aX+b)&=\int_{-\infty}^\infty (ax+b) f_X(x)\,dx \cr &=c\int_{-\infty}^\infty (ax+b) \exp \textstyle({-(x-\mu)^2\over 2\sigma^2})\,dx \cr &= a\cdot c\int_{-\infty}^\infty x \exp {\textstyle({-(x-\mu)^2\over 2\sigma^2}})\,dx +b\cdot c\int_{-\infty}^\infty x \exp (\textstyle{-(x-\mu)^2\over 2\sigma^2})\,dx \cr &=a\,\Bbb E(X)+b\cdot 1\cr &=a\mu+b. } $
The last integral being $1$, since we are integrating a density function.
It should be remarked this is one way the general formula $\Bbb E(aX+b)=a\Bbb E(X)+b$ can be obtained for continuous $X$. Nothing special about the normal distribution was used.