@Kitonum

Hi,

Please find below a set of equations I am trying to obtain the solution for. It has three lines:

Line 1: focdeltapioptS2Tbeta_eg :=

Line 2: focbetapioptS2Tbeta_eg :=

Line 3: optS2Tset := {solve({focbetapioptS2Tbeta_eg=0, focdeltapioptS2Tbeta_eg=0}, {beta, delta})};

Additional information that can be used in solving the two equations is is 0<= beta, delta <= 1. If the constraints cannot be incorporated explicitly for any reasons, then at least obtaining all possible sets of solutions for beta and delta would help.

Thank you,

Omkar

*************************************

focdeltapioptS2Tbeta_eg := 2.*delta*(-8.016120437*10^13*beta-5.599041156*10^11*beta^4-9.950892840*10^12*beta^3+4.172202042*10^9*beta^2*delta^6-3.593992717*10^11*beta^5*delta^2-5.153172141*10^12*beta^4*delta^2+1.100201852*10^9*beta^4*delta^8+3.667339507*10^8*beta^6*delta^8-3.667339507*10^8*beta^3*delta^8-1.100201852*10^9*beta^5*delta^8+2.670209307*10^9*beta^4*delta^4+1.032944471*10^10*beta^3*delta^4+2.114848520*10^13*beta^2*delta^2-1.562440115*10^13*beta^3*delta^2-6.630411895*10^9*beta^3*delta^6-1.582188745*10^10*beta*delta^4+1.713992189*10^9*beta^5*delta^6+7.442176619*10^8*beta^4*delta^6+2.822233434*10^9*beta^2*delta^4+2.464874192*10^10*beta*delta^2+2.000000000*10^10*delta^2-5.383607038*10^13*beta^2)/(-2636.981242*beta^2*delta^2+2636.981242*beta*delta^2-4108.123654*beta-10000.)^3:

focbetapioptS2Tbeta_eg := (.8333333333*(1.784795533*10^13*beta^4*delta^6-1.780394727*10^13*beta^5*delta^6-5.910514827*10^13*beta^2*delta^4-5.501009263*10^10*beta^8*delta^6+1.701903476*10^14*beta*delta^2+1.100201852*10^10*beta^9*delta^6+1.100201852*10^11*beta^7*delta^6+5.842965601*10^12*beta^6*delta^6-3.878841795*10^13*beta+2.205824535*10^11*beta^5+1.735006720*10^11*beta^4+4.159889289*10^10*beta^6+7.379516540*10^12*beta^3-6.234247772*10^11*beta^6*delta^4+5.141976568*10^10*beta^8*delta^4-1.126878340*10^12*beta^5*delta^2-3.757524927*10^13*beta^4*delta^4+1.556869544*10^13*beta^5*delta^4-8.051300147*10^10*beta^7*delta^4+1.496707827*10^11*beta^6*delta^2+8.010627926*10^10*beta^7*delta^2+1.459035087*10^13-9.619344526*10^13*delta^2+2.590921095*10^11*beta^2-5.952985789*10^12*beta^3*delta^6+2.536603101*10^13*beta*delta^4+1.142130779*10^13*beta^4*delta^2-3.282217618*10^13*beta^2*delta^2-5.169893249*10^13*beta^3*delta^2+5.639818880*10^13*beta^3*delta^4))/((-1.+beta)^2*(-2636.981242*beta^2*delta^2+2636.981242*beta*delta^2-4108.123654*beta-10000.)^3):

optS2Tset := {solve({focbetapioptS2Tbeta_eg=0, focdeltapioptS2Tbeta_eg=0}, {beta, delta})};