20/07/08 07:34:45.40 I3BoIViR.net
>>967
思考停止のプログラムでの数値解
> PQ <- function(xy){
+ x=xy[1]
+ y=xy[2]
+ P=c(x,x^2+1)
+ Q=c(2*y^2+2,y)
+ sqrt(sum((P-Q)^2))
+ }
>
> opt=optim(par=c(0,0),fn=PQ,method='Nelder')
> x=opt$par[1]
> (P=c(x,x^2+1))
[1] 0.6189828 1.3831398
> y=opt$par[2]
> (Q=c(2*y^2+2,y))
[1] 2.0814249 0.2017733
> PQ(opt$par)
[1] 1.87999