Tuesday, January 24, 2012

Thursday, January 05, 2012

"MAYOR" COTA INFERIOR del COMBUSTIBLE MÍNIMO usado

Para cada iteración de cálculo:

aa = RootOf(16*GG*MT*DTL^2-32*AA*DTL^2*vvi*DT*rri+96*AA*DTL*vvi^2*DT^2*rri+96*AA*DTL*vvi*DT*rri^2-32*GG*MT*DTL*vvi*DT+32*GG*MT*vvi*DT*rri-32*GG*ML*vvi*DT*rri-16*GG*ML*rri^2-16*AA*DTL^2*rri^2+32*AA*DTL*rri^3-16*AA*vvi^4*DT^4+16*GG*MT*rri^2-16*AA*rri^4-16*AA*DTL^2*vvi^2*DT^2+32*AA*DTL*vvi^3*DT^3-64*AA*vvi^3*DT^3*rri-96*AA*vvi^2*DT^2*rri^2-64*AA*vvi*DT*rri^3-32*GG*MT*DTL*rri+16*GG*MT*vvi^2*DT^2-16*GG*ML*vvi^2*DT^2+(8*DT^6*rri-AA*DT^8-4*DTL*DT^6+8*DT^7*vvi)*_Z^4+(-24*DTL*DT^4*rri-24*DTL*DT^5*vvi+4*DTL^2*DT^4-8*AA*DT^6*rri+24*DT^4*rri^2-8*AA*DT^7*vvi+4*AA*DTL*DT^6+48*DT^5*vvi*rri+24*DT^6*vvi^2)*_Z^3+_Z^5*DT^8+(24*AA*DTL*DT^4*rri+16*DTL^2*DT^3*vvi+16*DTL^2*DT^2*rri+32*DT^5*vvi^3-4*GG*ML*DT^4+32*DT^2*rri^3+96*DT^3*vvi*rri^2-48*DTL*DT^2*rri^2-48*AA*DT^5*vvi*rri+96*DT^4*vvi^2*rri-96*DTL*DT^3*vvi*rri-4*AA*DTL^2*DT^4+24*AA*DTL*DT^5*vvi+4*GG*MT*DT^4-48*DTL*DT^4*vvi^2-24*AA*DT^4*rri^2-24*AA*DT^6*vvi^2)*_Z^2+(64*vvi^3*DT^3*rri+96*vvi^2*DT^2*rri^2+64*vvi*DT*rri^3+16*rri^4+16*vvi^4*DT^4+16*DTL^2*rri^2-32*DTL*rri^3-32*AA*DT^5*vvi^3-32*AA*DT^2*rri^3+16*DTL^2*vvi^2*DT^2-32*DTL*vvi^3*DT^3+32*DTL^2*vvi*DT*rri-96*DTL*vvi^2*DT^2*rri-96*DTL*vvi*DT*rri^2-16*AA*DTL^2*DT^2*rri+48*AA*DTL*DT^4*vvi^2+48*AA*DTL*DT^2*rri^2-96*AA*DT^4*vvi^2*rri-96*AA*DT^3*vvi*rri^2-16*GG*MT*DTL*DT^2+16*GG*MT*DT^2*rri-16*GG*ML*DT^3*vvi-16*GG*ML*DT^2*rri-16*AA*DTL^2*DT^3*vvi+16*GG*MT*DT^3*vvi+96*AA*DTL*DT^3*vvi*rri)*_Z)

mg = DT*mm*AA/(EV+DT*AA)

vvf = aa*DT+vvi

rrf = (1/2)*aa*DT^2+vvi*DT+rri

NOTA: Podemos mejorar la cota haciendo DT tender a cero, lo cual incrementa el número de iteraciones.