dy/dx = by^a, y(0)=d