Hi Petr, In case you weren't aware, there's a Julia package for computation of arbitrary-order Gauss quadrature rules: https://github.com/billmclean/GaussQuadrature.jl
-Alex On Tuesday, December 16, 2014 11:31:01 PM UTC-6, Petr Krysl wrote: > > Hello everybody, > > In case you're interested, here is the implementation so far. Switch to > the top folder, and "include" one of the examples. > If you have paraview, you could check the graphics output (enable graphics > export in the source). > > As always, I'm keen on getting feedback, so don't hesitate please. > > https://github.com/PetrKryslUCSD/jfineale_for_trying_out > > Petr > >