I used SymPy package to get this solution. -oo is the wrong answer. The one that matlab gives is the correct answer.
On Wednesday, May 18, 2016 at 5:49:24 AM UTC-6, j verzani wrote: > > I'm assuming this is using the `SymPy` package. I can confim your -oo > answer, but have no sense if it is correct or not. I would guess that there > are some differences in parsing floating point numbers between MATLAB and > Julia that could give a difference here. I don't know if it is possible, > but if you used rational coefficients you might see the answer you want. > > On Wednesday, May 18, 2016 at 1:33:57 AM UTC-4, Sreenath Chalil Madathil > wrote: >> >> We have the below symbolic parameters with symbolic values as below >> ~~~~~~~~~~~ >> s >> Pr1, Pr2,...., Pr14 >> a = (-25.7727272727273 - 51.5454545454545*(0.0758377425044092 + >> 0.37037037037037/s)/((-0.0758377425044092 - 0.37037037037037/s)*(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s))))*(-0.0388007054673721*Pr4/s + >> 0.0388007054673721*Pr5/s + 3.0/s + 30.0/s^2) + (-25.7727272727273 - >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)/((- >> 0.0758377425044092 - 0.37037037037037/s)*(9.0 - 0.0388007054673721*( >> 4.90909090909091 + 19.0909090909091/s)/(-0.0758377425044092 - >> 0.37037037037037/s))))*(-0.0388007054673721*Pr2/s + 0.0388007054673721* >> Pr4/s) + 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 >> - 0.0388007054673721*(25.7727272727273 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr6/s + >> 0.0388007054673721*Pr7/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(51.5454545454545 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr7/s + >> 0.0388007054673721*Pr8/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(77.3181818181818 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr8/s + >> 0.0388007054673721*Pr9/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(103.090909090909 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(0.0388007054673721*Pr10/s - >> 0.0388007054673721*Pr9/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(128.863636363636 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr10/s + >> 0.0388007054673721*Pr11/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(154.636363636364 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr11/s + >> 0.0388007054673721*Pr12/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(180.409090909091 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr12/s + >> 0.0388007054673721*Pr13/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1.0 - >> 0.0388007054673721*(206.181818181818 - 1.0*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))/(9.0 - >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/(- >> 0.0758377425044092 - 0.37037037037037/s)))*(-0.0388007054673721*Pr13/s + >> 0.0388007054673721*Pr14/s)/(-0.0758377425044092 - 0.37037037037037/s) + >> 25.7727272727273*(0.0758377425044092 + 0.37037037037037/s)*(1 + >> 0.0388007054673721*(4.90909090909091 + 19.0909090909091/s)/((- >> 0.0758377425044092 - 0.37037037037037/s)*(9.0 - 0.0388007054673721*( >> 4.90909090909091 + 19.0909090909091/s)/(-0.0758377425044092 - >> 0.37037037037037/s))))*(-0.0388007054673721*Pr5/s + 0.0388007054673721* >> Pr6/s - 3.0/s - 30.0/s^2)/(-0.0758377425044092 - 0.37037037037037/s) - >> 1133.0*(0.0758377425044092 + 0.37037037037037/s)/(s*(-0.0758377425044092 >> - 0.37037037037037/s)*(9.0 - 0.0388007054673721*(4.90909090909091 + >> 19.0909090909091/s)/(-0.0758377425044092 - 0.37037037037037/s))) >> >> >> >> When we apply the code >> ~~~~~~~~~~~~~~ >> limit(s*a,s,0) we get the value of -∞ >> >> >> ~~~~~~~~~~~~~~~~ >> But when we apply the same code in matlab, we get the below result >> ~~~~~~~~~~~~~~ >> (10*Pr2)/11 - Pr4/11 - Pr6/11 - Pr7/11 - Pr8/11 - Pr9/11 - Pr10/11 - Pr11 >> /11 - Pr12/11 - Pr13/11 - Pr14/11 + 1052/11 >> >> >> >>
