#loss function1 x=seq(-2,2,length=50) y=seq(-2,2,length=50) lossfunc1=function(x,y){ return(x*y*exp(-x^2-y^2)) } z=matrix(0,50,50) for(i in 1:50){ for(j in 1:50){ z[i,j]=lossfunc1(x[i],y[j]) } } persp(x,y,z,ticktype="detailed",theta = 30, …

func1=function(w){ return(w^2-2) } dffunc1=function(w){ return(2*w) } ddffunc1=function(w){ return(2) } newton102=function(w,times){ cnt=1 for(i in 1:times){ w_before=w w=w-dffunc1(w)/ddffunc1(w) printf("cnt:%d;w=%.8f,error=%.8f\n",cnt,w,a…