computeNumericalGradient<-function(J,theta){ interval<-0.1 perturb<-0.1 loss1<-J(theta-perturb) loss2<-J(theta+perturb) numgrad<-(loss2-loss1)/(2*interval) return(numgrad) }