R help - Sep 2008 - FW: odesolve dynload example

HI R Gurus,
> This is my first foray into using c-code with R, so ...
> I had a look at the archives and did not find anything on this, so
> hopefully I am not doubling up.
> 
I have previously used the following approach where I needed some very
small numbers/large (using Brobdingnag):


surfacewithdiff <- function(t, y, p)
  {
    const=p["const"]
    kay =p["kay"]
    psii=p["psii"]
    m1=p["m1"]
    m2=p["m2"]
    M=p["M"]
    N=p["N"]
    T=p["T"]

    bexp<-as.brob((const*(psii-m1*y[1]-m2*M/N))/T)
    yd1<-as.numeric((kay*exp(bexp))/(kay*t*exp(bexp)+1)-
    	(kay^2*t*exp(2*bexp))/(kay*t*exp(bexp)+1)^2)
    	
    list(c(yd1))
  }


Now I guess this could be done with gmp.

Now I use this stuff in the function that odesolve's lsoda refers to,
and 
I have found that my approach using pure R code is not fast enough, so I
need to convert 
my function to C code. This is all on a windows system.

Is it possible to access gmp or Brobdingnag dll's from the C code?
Or do I need to install and compile the gmp library itself?
> Advice appreciated.
> 
> Kind regards, 
> 
> Matt Redding
********************************DISCLAIMER**************...{{dropped:15}}