Displaying 1 result from an estimated 1 matches for "sersic".
2010 Feb 09
2
Double Integral Minimization Problem
...ius_eff){
thita=x[1]
radiusb=x[2]
counts2=(1-ellip)*(vita-1)/(pi*alpha^2)*counts0*radiusb*exp(-(2*index-0.324)*(radiusb/radius_eff)^(1/index))*(1+(((majaxis_pix^2)+(radiusb^2)-2*majaxis_pix*radiusb*cos(-thita)+(ellip^2-2*ellip)*((radiusb*sin(thita))^2))/(alpha^2)))^(-vita)
return(counts2)
}
sersic=function(p)
{
counts0=p[1]
index=p[2]
radius_eff=p[3]
value1=adapt(2,c(0,0),c(2*pi,200),functn=intensity,minpts=1000,maxpts=NULL,eps=0.01,counts0=counts0,index=index,radius_eff=radius_eff)
test=value1$value
f=sum(((counts-test)/countserr)^2)
return(f)
}
out<-nlminb(c(16000.0, 3.0,10.0),...