Rによる多項式の解

xRange<-c(-3,3,0.01)
yRange<-c(-50,500)
coef<-c(-45,-4,3,-9,5) #y=45-4*x+3*x^2-9*x^3+5*x^4
##############################################
buf<-round(polyroot(coef),digits=5) 
polynomialEquation<-function(x){
  tmp<-0  
  for(ccc in 1:length(coef)){
    tmp<-tmp+coef[ccc]*x^(ccc-1)
  }  
  tmp<<-tmp  
}
x<-seq(xRange[1],xRange[2],by=xRange[3])
y<-sapply(x,polynomialEquation)
plot(x,y,type="l",ylim=yRange,panel.first=grid(nx=NULL,ny=NULL,lty=2,equilogs=T))
abline(h=0,v=0,col="red")
equation<-"Equation\n"
result<-"Result\n"
for(ccc in 1:length(coef)){
    equation<-paste(equation,coef[ccc],"*x^",(ccc-1),"\n",sep="")
    if(ccc<=length(buf)){  
      result<-paste(result,buf[ccc],"\n",sep="")
    }
}
text((xRange[1]+xRange[2])/2,yRange[2],equation,adj=c(1,1))
text(xRange[2],(yRange[1]+yRange[2])/2,result,adj=c(1,0))

2016年03月29日分 シカゴIMM筋ポジション・ロング・ショート・ネットポジション数の時系列推移

2016年03月第4週 対外及び対内証券売買契約等の状況・株式投資ファンド持分・中長期債・合計のネット額時系列推移

2016年03月22日分 シカゴIMM筋ポジション・ロング・ショート・ネットポジション数の時系列推移

2016年03月第3週 対外及び対内証券売買契約等の状況・株式投資ファンド持分・中長期債・合計のネット額時系列推移

2016年03月15日分 シカゴIMM筋ポジション・ロング・ショート・ネットポジション数の時系列推移

2016年03月第2週 対外及び対内証券売買契約等の状況・株式投資ファンド持分・中長期債・合計のネット額時系列推移