連立方程式の計算 一般化、大規模化のためにこちらの改良 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,…