エルミート多項式 http://en.wikipedia.org/wiki/Hermite_polynomials

HermitePolynomials.go

package main

import (
	"fmt"
	"math"
)

func HermitePolynomialsProb(x float64, n int) float64 {
	norm := float64(factorial(n))
	y := 0.0
	bound := math.Floor(float64(n) / 2.0)
	for m := 0; m <= int(bound); m++ {
		y += math.Pow(-1, float64(m)) * math.Pow(x, float64(n-2*m)) /
			(float64(factorial(m)) * float64(factorial(n-2*m)) *
				math.Pow(2, float64(m)))
	}
	return y * norm
}

func factorial(x int) int {
	if x == 0 {
		return 1
	}
	return x * factorial(x-1)
}

func main() {
	fmt.Println("Hello Hermite Polynomials.")
	for i := 0; i < 10; i++ {
		fmt.Println(HermitePolynomialsProb(0, i))
	}
}

scipy.special.hermitenorm 的なの

Derivative of Gaussianの計算に使う

http://play.golang.org/p/0ukIb_vASq
最初からこうすればよかった