#
#octave_polynomial_Hermite.m
#https://ix86.org/octave/index-main.html?octave_polynomial_Hermite
#
function H0 = octave_polynomial_Hermite( n )

	H0 = [1];
	H1 = [1,0];
	for i=1:n
		H0_1 = i * H0;
		H0_2 = conv ( [0,0,1], H0_1 );

		H2 = conv( [1,0], H1 ) - H0_2;
		
		H0 = H1;
		H1 = H2;
	end

	ZERO = roots( H0 );
	ZERO_MAX = max( ZERO );
	ZERO_MIN = min( ZERO );

	X = linspace ( ZERO_MIN, ZERO_MAX, n*10 );
	Y = polyval ( H0, X );
	plot( X, Y )

endfunction

