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
>
>
>
>

Reply via email to