Question: Intersection of two curves using fsolve command

Question: You will observe that the two curves intersect at one or more places in the displayed region. Use fsolve to find the (x; y) coordinates of all the points of intersection in this region.

In theory this question has a very simple answer, but i am struggling no end finding it. 

P1 := y^2 = x^3+x^2

P2 := y = 2*x+1

This is what the graph looks like when plotted, and as we have already got the two equasions assigned as H and J, i just type in fsolve({P1,P2})? 

Apparently not, as this returns {x = 4.048917340, y = 9.097834679}, which is correct, but is only one of three points of intersection, and just so happens to be the one not shown in my plot (although it is still a valid point of intersection) and i need the ones within the region i have set for the graph (x=-1..2 and y=-1..1).

How do i get maple to show me the two points of intersection i need using the fsolve command? 

I've tried rearranging the equasions for y, rearranging for x, adding an option to solve for x [fsolve({P1,P2},x) but this returned an error] and am pretty stumped. For something so seemingly simple i dont really know where to turn.

Thank you :)