Two things: In your DSolve[] command, change b to b[a]. Second, if it's not plotting properly, it's for two reasons. You're using a lin-log plot to plot the answer, and the answer has imaginary components. Also, you have a second order equation with only one boundary condition. I assume that...