Following is the problem I am trying to solve. This function will
generate
the weights for Gauss Hermite integration.
hermite_w[n](r):='block(
map(lambda([x],(2^(n-1))*(n!)*bfloat(sqrt(%pi))/((n*hermite_h[n-1](rhs(x)))^2)),r))$
So I execute
roots(hermite_h(4,x))
and get
(d7) [x = - 1.65068012388578d0, x = - 0.52464762327528d0, x =
0.52464762327528d0, x = 1.65068012388578d0]
Then I execute hermite_w[4](d7);
and get:
[((1.7724538509055160273b0 * 2^(n - 1) * n!)/((hermite_h[n - 1]( -
1.65068012388578d0))^2 * n^2)),
((1.7724538509055160273b0 * 2^(n - 1) * n!)/((hermite_h[n - 1]( -
0.52464762327528d0))^2 * n^2)), ((1.7724538509055160273b0 * 2^(n - 1) *
n!)/((hermite_h[n - 1](0.52464762327528d0))^2 * n^2)),
((1.7724538509055160273b0 * 2^(n - 1) * n!)/((hermite_h[n -
1](1.65068012388578d0))^2 * n^2))]
So for some reason n is not known in the function. What I want is that
this expression should evaluate to a list of floating point numbers
only.
Dan