Jacobian matrix - incorrect result calculation

Hi,

The 'sjacobi_rci_c.c' contained in the examples MKL library returns the error results:

x f1 f2

0.0000 0.0000 0.0000

0.0000 1.#INF 1.#INF

1.#INF 1.#INF 1.#INF

650751168283084220000000000000.0000 1.#INF 1.#INF

df(y0)= -431602080.0000 -431602080.0000 -431602080.0000 -431602080.0000

df(y1)= -431602080.0000 -431602080.0000 -431602080.0000 -431602080.0000

df(y2)= -431602080.0000 -431602080.0000 -431602080.0000 -431602080.0000

df(y3)= -431602080.0000 -431602080.0000 -431602080.0000 -431602080.0000

This problem does not occur for double-precision calculations:

x f1 f2

10.0000 0.0000 110.0000

10.0000 -2.2361 0.0000

10.0000 0.0000 100.0000

10.0000 0.0000 0.0000

df(y0)= 1.0000 0.0000 0.0000 0.0000

df(y1)= 10.0000 0.0000 -20.0000 0.0000

df(y2)= 0.0000 2.2361 40.0000 0.0000

df(y3)= 0.0000 -2.2361 0.0000 0.0000

I do not really know how to correct this example.

Thanks,

Stan