数値計算

• 連立方程式の計算
• 一般化、大規模化のためにこちらの改良
  • Rの逆行列を用いている

x<-c(0,0)
T<-100
d<-0.01

dim<-length(x)

f<-function(x){
z<-rep(0,dim)
z[1]<-(x[1]-3)^2+x[2]^2-3
z[2]<-sin(x[1])+exp(x[2]-1)-1
return(z)
}


g<-function(i){
	v<-rep(0,dim)
	v[i]<-1
	return(v)
	}

A<-matrix(0,nr=dim,nc=dim)


for(t in 1:T){

for(i in 1:dim){
for(j in 1:dim){
	A[i,j]<-(f(x+d*g(j))[i]-f(x-d*g(j))[i])/(2*d)
}}


x<-x-solve(A)%*%f(x)

}

print(x)
print(f(x))