Displaying 2 results from an estimated 2 matches for "psipol".
Did you mean:
phipol
2011 Jan 26
0
Bivariate polynomials in R
...pol')
psi <-function (order){
psi<-matrix(data=polynomial.coefficients(legendre.polynomials(order)[[order+1]]),ncol=1)
}
phi <-function (order){
phi<-matrix(data=polynomial.coefficients(legendre.polynomials(order)[[order+1]]),nrow=1)
}
# Then I convert them to multinomials.
psipol<-as.multipol(psi(2))
phipol<-as.multipol(phi(2))
print(psipol)
print(phipol)
# Here error occurs
psipol*phipol
The multiplication above fails.
According the multipol pdf
http://cran.r-project.org/web/packages/multipol/vignettes/multipol.pdf
this works: (which does)
require(multipol)
a &...
2011 Jan 27
0
adaptIntegral takes too much time
...,.9999))
-but this runs to infinity:
(I ll leave computer open all night just in case it terminates)
adaptIntegrate(f,lowerLimit=c(-1,-1),upperLimit=c(1,1))
I do not really want to bother you with the details of my code so do you might have any idea what might cause this?
f is
f <- as.function(psipol*phipol) # Multivariate polynomial to function
and
f(c(1,1))
returns immediately results.
Any clue what should I try to check?
I would like to thank you for your advice
Best Regards
Alex