Displaying 1 result from an estimated 1 matches for "grad2".
Did you mean:
grade
2007 Sep 22
0
error messages
...ubdivisions=1000000)$value)
))
}
gr<-function(u){
k1<-u[1]
k2<-u[2]
k12<-u[3]
g<-u[4]
grad1<-sum(integrate(function(y)integrate(function(a)(cos(x)-cos(a))*exp(k2*cos(y)+(k1+k12*cos(y))*cos(a)),subdivisions=10000000,-pi,pi)$value,subdivisions=10000000,g,pi)$value)
grad2<-sum(integral(function(y)cos(y)*exp((k2+k12*cos(x))*cos(y)),subdivisions=10000000,g,pi)$value*integrate(function(y)besselI(k1+k12*cos(y),0)*exp(k2*cos(y)),subdivisions=10000000,g,pi)$value-integrate(function(y)exp((k2+k12*cos(x))*cos(y)),subdivisions=10000000,g,pi)$value*integrate(function(y)bes...