			function psi = recpsi(Alist, V)
			N = length(Alist) + 1;

			% manage the trivial case, Alist is empty
			if N < 2
			psi = V;
			return;
			end

			indices = cell(1,N);

			% A1 index
			indices{1} = [-1, -2, 1];

			% Ak index
			for k = 2:(N-1)
			indices{k} = [k-1, -(k+1), k];
			end

			% V index
			indices{N} = [N-1, -(N+1), -(N+2)];

			% reconstruct of psi
			psi = ncon({Alist{:}, V}, indices);
			end

compphys/exo4/recpsi.m
09de8b211167master

recpsi.m
View Options