Thank you Robert and hirnyk for taking your time replying.

Iam exploring other things starting from ur codes. But since I have to submit something urgently by the end of March (and theoretical calculations by hand is taking alot of time), so can u please help me just a bit more.

I have 2 polynomial functions, say, for ex.

f(x,y) = axy + bx - c*x*(x+y)

g(x,y) = dxy + ey - f*y*(x+y)

I hav to find the roots, say they are xr and yr.

Now I need to construct a matrix J (called the Jacobian) such that the elements are ::

J11= ∂f/dx (i.e differentiating f w.r.t. x, & considering y constant)

J12= ∂f/dy 

J21= ∂g/dx

J22= ∂g/dy

and then I need to substitute the value x=xr & y=yr in the above 4 terms. This will generate the Jacobian matrix.

And finally I need to calculate the Eigenvalue of this Matrix J.

Although I have larger number of equations in my system, but from the above example, I will be able to generalize.

Immediate help will be very greatefully appreciated.

Regards

Alice