Reciban un cordial saludo.

De manera muy comedida les pido ayuda con la función indicada abajo. La usan
en un artículo de inventarios. No veo dónde estoy errado.

Gracias de antemano.

César Escalante C.

> lamb<-50
> L<-1
> h.<-10
> p<-25
> K<-5
> integ<-function(y){integrate(function(x)ppois(x,lamb*L),0,y)$value}
> G<-function(y){(h.+p)*integ(y) + p*(lamb*L-y)}
> G(0:43)
 [1] 1250 1225 1200 1175 1150 1125 1100 1075 1050 1025 1000  975  950  925
900  875  850  825  800  775  750  725  700  675  650  625  600  575  550
525
[31]  500  475  450  425  400  375  350  325  300  275  250  225  200  175
> G(44)
Error en integrate(function(x) ppois(x, lamb * L), 0, y) :
  maximum number of subdivisions reached
> G(500)
[1] 4499.785

	[[alternative HTML version deleted]]

