newton_quiz.m
View Options

			function newton(N)

			A = (sparse(2:N,1:N-1,-1,N,N)+ sparse(1:N,1:N,2,N,N) + sparse(1:N-1,2:N,-1,N,N));
			b = ???;
			x_0 = ones(N,1);
			u = x_0;
			dd = sparse(N,N);
			discrep=1;
			k=0;
			maxIt=10000;

			while(discrep>0.000001 && k < maxIt)
			for i=1:N
			dd(i,i)=???;
			end
			mat=A+dd;
			rhs=???;
			y=mat\rhs;
			discrep=norm(y);
			u=u+y;
			k=k+1;
			end

			discrep
			k
			u


			end

master-raffo/advnum/newton_quiz.m8daa7c7f64acmaster