Question: How to produce implicit plot within loop

Here is my try:

with(plots):
epsilon:=5*Pi/2:
gammaa:=0.4:
gammao:=0.1:
gammab:=0.002:
omegab:=100:
X:=60:
Omega:=6*Pi:
g:=10:
sigmag:=5*gammaa:
a:=1/(sigmag*sqrt(2*Pi)):
for deltao from 0 to 4 do
D7:=evalf(a*int(g*Omega*(gammao/2-I*Delta)*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2),Delta=-infinity..infinity)):
D8:=evalf(0.5*a*int(g^2*gammao*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
D9:=evalf(a*int(g^2*Delta*exp(-0.5*((Delta-deltao)/sigmag)^2)/(gammao^2/4+Delta^2), Delta=-infinity..infinity)):
b1:=deltao-(2*X^2*omegab*Y^2)/(omegab^2+gammab^2)-D9:
f:=epsilon-D7=Y*(gammaa+D8+I*(b1)):
P1:=implicitplot(f,Y^2=0..10,delta0=0..1,numpoints=1000,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[4, 3],linestyle=1);
end do:

display(P1);

The function (f)  is an implicit function of Y^2

so my code above is an example to plot Y^2 against deltao

note that Y is complex

I appreciate valueble comments.