(* dotprod.wl -- Dot Product for Polynomials *) dp=Function[{p,q},Integrate[p*q,{x,-1,1}]] nm=Function[p,Sqrt[dp[p,p]]] dp[x^2+x-1,x^4] nm[x^2+x-1]