Hi!

I'm having some trouble with the gauleg function. First the compiler complains that the arguments for arth are wrong type. I changed the 1 to 1.0 and got past that (though I am not certain if it meant something deeper problems) and the program compiled just fine. When I ran it though, I got segmentation fault. I found out that the problem is with the last code line w(n:n-m+1:-1)=w(1:m)
I call gauleg with
CALL gauleg(0.0,z,zp,dzp)
where z is of course REAL and zp and dzp are DIMENSION(np) with np=25.

Any help would be appreciated!