v 10.1

P=ε0(χE+χ2E2+χ3E3+)