based on an analysis of samples from natural populations of drosophila pseudoobscura collected from three localities on mount san jacinto ( california ) , dobzhansky et al . ( 1942 ) showed that homozygotes for genes located on the second and third chromosomes in this species had a lower developmental rate than expected . these authors suggested that , under controlled culture conditions , the viability and developmental rates were dependent on the demographic density of the larvae . in contrast to the genetic finding , the latter observation was subsequently confirmed by several authors ( bakker , 1959 , 1969 ; de witt , 1960 ; robertson , 1963 , 1964 ; barker and podger , 1970 ; huang et al . , 1971 ; tosic and ayala , 1981 ; mather and caligari , 1981 ; mnsua and moya , 1983 ; adell et al . , 1988 ) , not only for drosophilid populations but for other insects as well . as a result , all subsequent research on this subject has centered on the analysis of ecological models based on competition among pre - adult individuals . oliveira and cordeiro ( 1981 ) detected a major gene effect on delayed development in drosophila melanogaster that was later associated with a recessive gene located on the second chromosome of this fly ( oliveira et al . , 1991 ) . since the development of molecular genetics in the early 1990s many authors have examined the role of genes involved in the regulation of developmental timing and body size determination in drosophila . most of these studies focused on genes related to apoptosis and to the production or release of ecdysone and other ecdysteroids ( molting hormones ) . the results of some relevant articles ( all dealing with this hormone in d. melanogaster ) are summarized below . other authors described the importance of additional gene - controlled systems , such as those involved in the regulation of adenosine deaminase - related growth factor adgf - a , in developmental timing ( dolezal et al . , 2005 ) and in the regulation of the activities of oxidative phosphorylation complexes through nuclear - encoded mitochondrial tyrosyl - trna synthetases ( meiklejohn et al . , 2013 ) . oldham et al . ( 2000 ) and quinn et al . ( 2012 ) provide appropriate reviews of this subject . sliter and gilbert ( 1992 ) showed that loss of function mutations at the dre4 gene caused ecdysteroid deficiency and developmental delay or arrest . 2002 ) showed that isoform - specific null mutants for the orphan member e75a gene had a reduced ecdysteroid titer that resulted in developmental delay . 2007 ) showed that prothoracicotropic ( ptth ) hormone production was not essential for molting or metamorphosis but that its lack resulted in significantly delayed larval development . ( 2010 ) showed that the gap gene giant ( gt ) regulated ecdysone production through specification of ptth and that mutants for this gene exhibited delayed larval development ( in addition to other effects ) . the purpose of the present theoretical investigation was to analyze the population genetic effects of a delay in the larval developmental rate produced by an autosomal recessive mode of inheritance . using elements from the theory of finite differences , difference and differential equations , and dynamic systems , we examined the problem based on a model with discrete generations and with the total number of individuals ( adults plus larvae ) kept constant . we also summarize the results obtained with alternative models involving an increase or decrease in the total population number . the case of continuous generations ( to be presented in detail elsewhere ) is also briefly discussed . the central principle of this model is that the total number of individuals ( adults + larvae ) is kept constant in a population with discrete generations . for this , let a and a be a pair of segregating alleles at an autosomal locus . in the homozygous state the recessive gene a induces a delay in larval development ( increase in the diapause phase ) , measured in an integer number of generations : crosses therefore take place randomly ( panmixia ) only among individuals belonging to the same generation . the population size , kept rigorously constant throughout generations , is assumed to be so large that the effects of random sampling genetic drift are negligible . the effects of differential migrations , selection pressure and mutation are also considered negligible . in the equations describing the model we will use the following parameters : n0 = number of adults in the initial population at generation 0 ; this number equals the total number of individuals ( adults plus larvae ) at any time or generation k  0 ; q0 = initial frequency of the recessive allele a [ q0 = f0(a ) ] ; p0 = initial frequency of the normal allele a [ p0= f0(a ) = 1q0 ] ; t = delay in larval development determined by the genotype aa , measured in an integer number of generations ; k = index corresponding to an integer number of generations ; nk = number of adult individuals in the population at the k - th generation ; nk = number of larvae in the population at the k - th generation ; qk = frequency of the recessive gene a among adults at the k - th generation ; pk = frequency of the normal allele , also among adults belonging to the k - th generation ; qk = frequency of the recessive gene in the total population ( adults plus larvae ) at the k - th generation ; pk = frequency of the normal allele , also in the total population , at the k - th generation ( pk= 1  qk ) . based on the foregoing definitions : ( a ) the total number of adults at generation k is given by 
 nk = nk1(1qk12 )  if  ktor 
 nk = nk1(1qk12)+nk1tqk1t2  if  k > t ; ( b ) the genotypic proportions among aa , aa and aa adult individuals at generation k are , respectively : 
 ( nk1/nk)pk12,(nk1/nk)2pk1qk1 ,  and  0  if  ktor 
 ( nk1/nk)pk12,(nk1/nk)2pk1qk1 ,  and(nk1t / nk)qk1t2  if  k > t ; ( c ) the number of larvae aa present at generation k after the maturation of normal individuals is given by 
 nk = n0q02+n1q12+  +nk1qk12  if  kt  ornk = nktqkt2+nkt+1qkt+12+ +nk1qk12  ifk > t . using the foregoing equations , it is not difficult to show that nk+1 + nk+1= nk+ nk= ... = n0 + 0 = n0 and that , if pk= f(a ) in the total population , then pk= ( nk / n0)pk= ( nk / n0)(n0p0/nk ) = p0 , i.e. , allele frequencies and the number of individuals remain constant when the whole population ( adults plus larvae ) is considered , as assumed by the model . for k t , we have , among adults of the population , a mechanism similar to the classic model of total selection against homozygous recessive individuals aa , despite the absence of selection against any individual of the population ( when adults and larvae are considered together ) , as shown by the equations derived in the previous paragraph . this results in qk+1= qk/(1+qk ) , a recurrence equation that has the general solution in simple analytical form qk= q0/(1+kq0 ) . for k > t , we have 
 nk+1=nk(1qk2)+nktqkt2=  nk(1qk)(1+qk)+nktqkt2=  nkpk(2pk)+nkt(1pkt)2 . since pk+1= ( nkpk+nkpkqk)/nk+1= nkpk / nk+1 and , hence , nk= n0p0/pk , it follows that n0p0/pk+1 = n0p0(2pk ) + n0p0(1pkt)/pkt and 
 1/pk+1=1/pkt+pktpk , pk+1/pk+1=pkt+1/pkt , or , since pk= n0p0/nk , we also obtain 
 nk+1=nkt+(n0p0)2(1/nkt1/nk ) . below , we analyze the behavior of allele frequencies pk and qk . for k  t , pk= ( p0 + kq0)/(1 + kq0 ) , i.e. , pk is a monotonically increasing function of k , whereas for k > t , 1/pk+1= pkt+ 1/pkt pk . x < 1 shows that f(x ) = 1  1/x < 0 ; f(x ) is therefore a monotonically decreasing function of x in this interval . if pk > pkt , then we have pk+ 1/pk < pkt+ 1/pkt= pk+ 1/pk+1 and , therefore , pk+1 < pk ; similarly , pk < pkt implies that pk+1 > pk , and pk= pkt implies that pk+1= pk . in the interval 0 < k  t , pk increases such that when k surpasses t , we have pk > pkt ; pk then decreases to a new point where pk < pkt is reached and we will then have another increase in pk , and so on . thus , there will be an oscillation in gene frequencies among adults of the population . with the passing of generations , i.e. , when k increases , the variations become smaller . in the limit case , when k tends to infinity , the gene frequencies reach an equilibrium ( pinf= pe , qinf= qe ) . figures 13 show numerical examples of oscillatory convergence to equilibrium points ( pe ) over 10 generations of random crosses , for populations starting with initial gene frequencies p0 = 0.1 , ... , 0.9 and t values of t = 1 , 2 , and 3 . these graphs and others shown in this paper were prepared using routines and packages from mathematicav.8.0 ( wolfram research ) . the preliminary global analysis presented above already shows a point of real theoretical interest , namely , that despite intense selective pressure among adult individuals ( but not among the total population of adults plus larvae ) , the recessive allele is never eliminated from the population , but rather is kept in equilibrium at a point qe > 0 . this is true except when the number of generations involved in the developmental delay is very large or tends to infinity such that allele a is practically eliminated from the adult population and all larvae fail to develop . in a generic instant we have the following population composition ( figure 4 ) . if we consider an equilibrium situation where 
 nk+1=nk = nk1=nk2=  =nkt = neqk+1=qk = qk1=qk2=  =qkt = qe=1pe , from nk+ nk= n0 it follows that ne+ tneqe= n0 . since ne= n0p0/pe , it follows that n0p0(1+tqe)/(1qe ) = n0 , and , hence , p0 + tp0qe= 1qe and , finally , 
 qe=[(1 + 4tp0q0)1/21]/(2tp0)=  2q0/[(1 + 4tp0q0)1/2 + 1 ] . as expected , qe= q0 if t = 0 ( simple case of panmixia without selection ) and qe= 0 if t =  ( special case of complete selection against recessive individuals aa ) . considering that : 
 pk = nk1pk1/nk = n0p0/nk  for  any  k , nk = nk1(1qk12 )  for  any  kt , and 
 nk = nk1(1qk12)+nkt1qkt12  for  any  k > t , it is not difficult to show that the recurrence equation 1/pk+1= pkt+1/pktpk can be rewritten in the more convenient form 
 1/pk=1/p0[j=1,t]{(1pkj)2/pkj}. this recurrence equation allows straightforward derivation of the equilibrium point pe . indeed , when k tends to infinity , the limit for the right - hand side of the above equation is 1/p0 t(1pe)/pe and it follows that pe= p0(1+tqe ) and pe= 1  qe= 1  [ ( 1 + 4tp0q0)1]/(2tp0 ) . table 1 shows the sets of equilibrium points qe for several values of t , as a function of the initial values q0 = 0.0 , 0.1 , ... , 0.9 , 1.0 . figure 5 shows the sets of equilibrium values qe as a function of the initial frequency values q0 and t varying from 0 to 15 . for the case t = 1 , the recurrence equation above can be rewritten as 
 1/pk+1=1/p0(1pk)2/pk  or  pk+1=  p0pk/[pkp0(1pk)2 ] . the eigenvalue r responsible for the convergence of the series { p0 , ... , pe } is obtained by directly differentiating pk+1= f(pk ) at equilibrium point pe : 
 r = dpk+1/dpk|pk = pe=  p02(qe22qe)/(pep0qe2)2=pe21 . since |r| = |pe1| < 1 for any value of pe in the interval 0 < pe < 1 , the equilibrium point pe is asymptotically stable and convergence to this point occurs at a geometric rate . since r is always smaller than zero in this interval , convergence to equilibrium is oscillatory . figure 6 shows the values of r and qe as functions of the initial gene frequency q0 for the case t = 1 . for a generic value of t the method that follows is used for determining convergence rates as well as equilibrium properties . as shown above , for any k > t , 
 1/pk=1/p0[j=1,t]{(1pkj)2/pkj }  orpk=1/[1/p0[j=1,t]{(1pkj)2/pkj } ] . we then define a set of t variables as follows : 
 x1(k)pkx2(k)pk1 .. xk= (  .  ) ..... xt(k)pkt+1 since 
 pk+1=[1/p0[j=1,t]{(1pkj+1)2/pkj+1}]1we have : 
  considering the equilibrium point , at which 
 x1(k)pex2(k)pe .. xe= (  . ) =(  .  ) , .... xt(k)pewe have , in an infinitesimal neighborhood of this point , 
 where 
 a = df / dx1(k)|x = xe = df / dx2(k)|x = xe= the eigenvalues r1 , r2 , ... , rt are the solutions of the characteristic equation 
 det(ari)=(1)t(1+r+r2+ +rt1rt / a)=0or 
 rt / a=1+r+r2+r3+  +rt2+rt1,which can also be rewritten as r / a = ( 1r)/(1r ) , r  1 . based on the function f(r ) = r ( a+1)r+ a , r  1 , 1 < a < 0 , it can be shown that f(0 ) = a < 0 , f(1 ) = 0 and df(r)/dr = f(r ) = ( t+1)r t(a+1)r . also , as r tends to +  , f(r ) tends to + inf and as r tends to  , f(r ) tends to +  if t is odd , and to  if t is even . when t is even , f(r ) is an increasing function of r in the open interval (  , 0 ) , decreasing in the open interval ( 0 , ( a+1)t/(t+1 ) ) and again increasing for r > 0 . therefore , r = 0 is a relative maximum and r = ( a+1)t/(t+1 ) a relative minimum . when t is odd , f(r ) is a decreasing function of r in the open interval ( inf , ( a+1)t/(t+1 ) ) and an increasing function of r for r > ( a+1)t/(t+1 ) . therefore , r = 0 in the case when t is odd is an inflection point , whereas r = ( a+1)t/(t+1 ) is a relative minimum ( as in the case when t is even ) . hence , we conclude that when t is even , the equation rt+1 ( a+1)rt+ a = 0 , r  1 , 1 < a < 0 does not have any real root since f(r ) tends to  as r tends to  ; when t is odd , the above equation has just one real negative root since f(r ) tends to +  as r tends to . for the cases t = 1 to 4 it is possible to obtain explicit solutions in simple analytical form for the values of the eigenvalues r ; for t > analytical procedures applied in the case t < 5 and extensive numerical analysis of the equation for any values of t > 4 have shown that ( a ) when t is odd , there is always one real eigenvalue , with |ri| < 1 and 1 < ri < 0 , and ( t1)/2 pairs of complex conjugates rj , k = a  ib , with ( a+b ) < 1 for any initial gene frequency p0 ; ( b ) when t is even , there are always only t/2 pairs of complex conjugates with the same properties as stated before . we therefore conclude that the set of equilibrium points pe= f(p0 , t ) is a stable one , with convergence always occurring to it for any initial value p0 in the open interval ( 0 , 1 ) and for any integer value of t. based on the foregoing definitions : ( a ) the total number of adults at generation k is given by 
 nk = nk1(1qk12 )  if  ktor 
 nk = nk1(1qk12)+nk1tqk1t2  if  k > t ; ( b ) the genotypic proportions among aa , aa and aa adult individuals at generation k are , respectively : 
 ( nk1/nk)pk12,(nk1/nk)2pk1qk1 ,  and  0  if  ktor 
 ( nk1/nk)pk12,(nk1/nk)2pk1qk1 ,  and(nk1t / nk)qk1t2  if  k > t ; ( c ) the number of larvae aa present at generation k after the maturation of normal individuals is given by 
 nk = n0q02+n1q12+  +nk1qk12  if  kt  ornk = nktqkt2+nkt+1qkt+12+ it is not difficult to show that nk+1 + nk+1= nk+ nk= ... = n0 + 0 = n0 and that , if pk= f(a ) in the total population , then pk= ( nk / n0)pk= ( nk / n0)(n0p0/nk ) = p0 , i.e. , allele frequencies and the number of individuals remain constant when the whole population ( adults plus larvae ) is considered , as assumed by the model . for k t , we have , among adults of the population , a mechanism similar to the classic model of total selection against homozygous recessive individuals aa , despite the absence of selection against any individual of the population ( when adults and larvae are considered together ) , as shown by the equations derived in the previous paragraph . this results in qk+1= qk/(1+qk ) , a recurrence equation that has the general solution in simple analytical form qk= q0/(1+kq0 ) . for k > t , we have 
 nk+1=nk(1qk2)+nktqkt2=  nk(1qk)(1+qk)+nktqkt2=  nkpk(2pk)+nkt(1pkt)2 . since pk+1= ( nkpk+nkpkqk)/nk+1= nkpk / nk+1 and , hence , nk= n0p0/pk , it follows that n0p0/pk+1 = n0p0(2pk ) + n0p0(1pkt)/pkt and 
 1/pk+1=1/pkt+pktpk , pk+1/pk+1=pkt+1/pkt , or , since pk= n0p0/nk , we also obtain 
 nk+1=nkt+(n0p0)2(1/nkt1/nk ) . below , we analyze the behavior of allele frequencies pk and qk . for k  t , pk= ( p0 + kq0)/(1 + kq0 ) , i.e. , pk is a monotonically increasing function of k , whereas for k > t , 1/pk+1= pkt+ 1/pkt pk . analysis of the function f(x ) = x + 1/x in the interval 0 < x < 1 shows that f(x ) = 1  1/x < 0 ; f(x ) is therefore a monotonically decreasing function of x in this interval . if pk > pkt , then we have pk+ 1/pk < pkt+ 1/pkt= pk+ 1/pk+1 and , therefore , pk+1 < pk ; similarly , pk < pkt implies that pk+1 > pk , and pk= pkt implies that pk+1= pk . in the interval 0 < k  t , pk increases such that when k surpasses t , we have pk > pkt ; pk then decreases to a new point where pk < pkt is reached and we will then have another increase in pk , and so on . thus , there will be an oscillation in gene frequencies among adults of the population . with the passing of generations , i.e. , when k increases , the variations become smaller . in the limit case , when k tends to infinity , the gene frequencies reach an equilibrium ( pinf= pe , qinf= qe ) . figures 13 show numerical examples of oscillatory convergence to equilibrium points ( pe ) over 10 generations of random crosses , for populations starting with initial gene frequencies p0 = 0.1 , ... , 0.9 and t values of t = 1 , 2 , and 3 . these graphs and others shown in this paper the preliminary global analysis presented above already shows a point of real theoretical interest , namely , that despite intense selective pressure among adult individuals ( but not among the total population of adults plus larvae ) , the recessive allele is never eliminated from the population , but rather is kept in equilibrium at a point qe > 0 . this is true except when the number of generations involved in the developmental delay is very large or tends to infinity such that allele a is practically eliminated from the adult population and all larvae fail to develop . in a generic instant we have the following population composition ( figure 4 ) . if we consider an equilibrium situation where 
 nk+1=nk = nk1=nk2=  =nkt = neqk+1=qk = qk1=qk2=  =qkt = qe=1pe , from nk+ nk= n0 it follows that ne+ tneqe= n0 . since ne= n0p0/pe , it follows that n0p0(1+tqe)/(1qe ) = n0 , and , hence , p0 + tp0qe= 1qe and , finally , 
 qe=[(1 + 4tp0q0)1/21]/(2tp0)=  2q0/[(1 + 4tp0q0)1/2 + 1 ] . as expected , qe= q0 if t = 0 ( simple case of panmixia without selection ) and qe= 0 if t =  ( special case of complete selection against recessive individuals aa ) . considering that : 
 pk = nk1pk1/nk = n0p0/nk  for  any  k , nk = nk1(1qk12 )  for  any  kt , and 
 nk = nk1(1qk12)+nkt1qkt12  for  any  k > t , it is not difficult to show that the recurrence equation 1/pk+1= pkt+1/pktpk can be rewritten in the more convenient form 
 1/pk=1/p0[j=1,t]{(1pkj)2/pkj}. this recurrence equation allows straightforward derivation of the equilibrium point pe . indeed , when k tends to infinity , the limit for the right - hand side of the above equation is 1/p0 t(1pe)/pe and it follows that pe= p0(1+tqe ) and pe= 1  qe= 1  [ ( 1 + 4tp0q0)1]/(2tp0 ) . table 1 shows the sets of equilibrium points qe for several values of t , as a function of the initial values q0 = 0.0 , 0.1 , ... , 0.9 , 1.0 . figure 5 shows the sets of equilibrium values qe as a function of the initial frequency values q0 and t varying from 0 to 15 . for the case t = 1 , the recurrence equation above can be rewritten as 
 1/pk+1=1/p0(1pk)2/pk  or  pk+1=  p0pk/[pkp0(1pk)2 ] . the eigenvalue r responsible for the convergence of the series { p0 , ... , pe } is obtained by directly differentiating pk+1= f(pk ) at equilibrium point pe : 
 r = dpk+1/dpk|pk = pe=  p02(qe22qe)/(pep0qe2)2=pe21 . since |r| = |pe1| < 1 for any value of pe in the interval 0 < pe < 1 , the equilibrium point pe is asymptotically stable and convergence to this point occurs at a geometric rate . since r is always smaller than zero in this interval , convergence to equilibrium is oscillatory . figure 6 shows the values of r and qe as functions of the initial gene frequency q0 for the case t = 1 . for a generic value of t the method that follows is used for determining convergence rates as well as equilibrium properties . as shown above , for any k > t , 
 1/pk=1/p0[j=1,t]{(1pkj)2/pkj }  orpk=1/[1/p0[j=1,t]{(1pkj)2/pkj } ] . we then define a set of t variables as follows : 
 x1(k)pkx2(k)pk1 .. xk= (  .  ) ..... xt(k)pkt+1 since 
 pk+1=[1/p0[j=1,t]{(1pkj+1)2/pkj+1}]1we have : 
  considering the equilibrium point , at which 
 x1(k)pex2(k)pe .. xe= (  .  ) =(  .  ) , .... xt(k)pewe have , in an infinitesimal neighborhood of this point , 
 where 
 a = df / dx1(k)|x = xe = df / dx2(k)|x = xe=  =  = df / dxj(k)|x = xe=  = pe21 . the eigenvalues r1 , r2 , ... , rt are the solutions of the characteristic equation 
 det(ari)=(1)t(1+r+r2+ +rt1rt / a)=0or 
 rt / a=1+r+r2+r3+  +rt2+rt1,which can also be rewritten as r / a = ( 1r)/(1r ) , r  1 . based on the function f(r ) = r ( a+1)r+ a , r  1 , 1 < a < 0 , it can be shown that f(0 ) = a < 0 , f(1 ) = 0 and df(r)/dr = f(r ) = ( t+1)r t(a+1)r . also , as r tends to +  , f(r ) tends to + inf and as r tends to  , f(r ) tends to +  if t is odd , and to  if t is even . when t is even , f(r ) is an increasing function of r in the open interval (  , 0 ) , decreasing in the open interval ( 0 , ( a+1)t/(t+1 ) ) and again increasing for r > 0 . therefore , r = 0 is a relative maximum and r = ( a+1)t/(t+1 ) a relative minimum . when t is odd , f(r ) is a decreasing function of r in the open interval ( inf , ( a+1)t/(t+1 ) ) and an increasing function of r for r > ( a+1)t/(t+1 ) . therefore , r = 0 in the case when t is odd is an inflection point , whereas r = ( a+1)t/(t+1 ) is a relative minimum ( as in the case when t is even ) . hence , we conclude that when t is even , the equation rt+1 ( a+1)rt+ a = 0 , r  1 , 1 < a < 0 does not have any real root since f(r ) tends to  as r tends to  ; when t is odd , the above equation has just one real negative root since f(r ) tends to +  as r tends to . for the cases t = 1 to 4 it is possible to obtain explicit solutions in simple analytical form for the values of the eigenvalues r ; for t > analytical procedures applied in the case t < 5 and extensive numerical analysis of the equation for any values of t > 4 have shown that ( a ) when t is odd , there is always one real eigenvalue , with |ri| < 1 and 1 < ri < 0 , and ( t1)/2 pairs of complex conjugates rj , k = a  ib , with ( a+b ) < 1 for any initial gene frequency p0 ; ( b ) when t is even , there are always only t/2 pairs of complex conjugates with the same properties as stated before . we therefore conclude that the set of equilibrium points pe= f(p0 , t ) is a stable one , with convergence always occurring to it for any initial value p0 in the open interval ( 0 , 1 ) and for any integer value of t. in order to study the effect of an increase or decrease in population size , it is enough to introduce the additional parameter c : n1p1= cn0p0 into equation n1p1= n0p0 of the previous section . consideration of this parameter leads straightforwardly the equation pk+1= 1/[2pk+(1pkt)/(ctpkt ) ] , k > t. this is the same equation obtained in the analysis of the constant population size model if c = 1 . at equilibrium , pk+1= pk= pkt= pe and the above equation reduces to pe(1pe)(1ct ) = 0 , c  1 . when c  1 , the only possible solutions for the foregoing equation are pe = 0 and pe= 1 , i.e. , the normal allele is completely lost or fixed . intuitively , it is easy to show that pe= 0 when c < 1 ( decreasing population size ) . thus , for any k > t , the proportion of recessive genes introduced by larvae ending their diapause period will become increasingly greater as generations go by , in relation to the normal allele present in the adult population , the size of which is decreasing . the inverse , i.e. , when c > 1 ( population increasing in size ) , is also obviously true since proportionally fewer recessive genes will be introduced into the adult population by emerging larvae . let p(t ) , r(t ) , q(t ) and l(t ) be the normalized proportions { p(t ) + r(t ) + q(t ) + l(t ) = 1 } of adults aa , aa and aa , and larvae aa in generation t. among adults , the normalized genotypic proportions are given by : 
 f(aa)=p=p/(1l)f(aa)=r=r/(1l)f(aa)=q=q/(1l)and the corresponding allelic frequencies f(a ) and f(a ) by 
 f(a)=p(t)=p+r/2f(a)=q(t)=q+r/2 . in the total population ( adults + larvae ) , after a time interval dt , the frequency of allele a among adults of the population is given by 
 p(t+dt)=[n(t).p2(t)+n(t).2p(t)q(t)/2]/n(t+dt)=={n(t).p(t)[p(t)+q(t)]}/n(t+dt)==n(t).p(t)/n(t+dt),with the restriction n(t+dt).p(t+dt ) = n(t).p(t ) = n0p0 = constant , as in the discrete model . the expression for p(t+dt ) can be used for quantitative analysis ( to be presented elsewhere ) in the case of continuous generations . one expects that equilibrium values for any fraction 0 <  < 1 in t < t+ < t+1 would interpolate between the equilibrium curves that represent the set of equilibrium points for the integer cases of t and t+1 generations in delayed larval development .
observations by dobzhansky s group in the 1940s suggesting that the presence of recessive genotypes could account for lower larval developmental rates in drosophila melanogaster were not confirmed at the time and all subsequent investigations on this subject focused on the analysis of ecological models based on competition among pre - adult individuals . 
 however , a paper published in this journal in 1991 eventually confirmed the finding made by dobzhansky and his co - workers . in this report 
 , we provide a theoretical analysis of the population genetic effects of a delay in the rate of larval development produced by such a genetic mechanism .