Recent content by Rositta

Hi, does anyone know how to calculate the total absorptin of a sample from its labdam-absorbance spectrum? I found a paper where the author find it as the average of the data, another one numerically integrate the spectrum over the whole lambda range.. but the two methods give different...

well.. the code worked with no problems at all in gfortran..

Hi TheoMcCloskey.. I've tried the modification with no change in the output. I do have a data file which I can read in Matlab as I mentioned in my first post. The flag report gives st=0. I found another Physics Forum thread talking about the same problem...

Thanks TheoMcCloskey. It is "iostat" in my original code. Apparently I made a mistake during pasting the code. I still have the problem of getting zero valued entities instead of the saved numbers.

Well.. the file is still unread.. I think I should post the codes=> =>The program I used to create the data is program xytable ! this code is to create data for the interpolation code written in console12 implicit none integer::i,j integer,parameter::m=50...

Thanks a lot for the quick answer.. I'll see if that will solve the problem..

Thanks a lot TheoMcCloskey. I have done some searching to fully understand the function of the END specifier with little success. Could you please expand a bit on the use of END?

I'm new to Microsoft Visual Studio (2010) & I'm a bit confused by the way Visual Studio handles fortran files. I saved a set of data in one console and tried to read the data from the created file in another console. I receive the next message: file opened correctly! forrt1: severe<24>...

Thaks a lot AlephZero. All the eigenvalues are calculated correctly now :)

Thaks all

Sorry.. Here is the main programme: =>program xggauss IMPLICIT NONE ! declare variables for part1 integer:: nval,i,k0 integer,parameter::kf=20 real*8,parameter::pi=3.141592653589 real*8:: x1,x2,dx,func,ss,x,r,del complex*8::sphshift ...

I'm using zheev subroutine in a do loop (up to 2000 times) that contains other subroutines and functions. The matrix to diagonalize is only 5-by-5 however zheev gives the correct eigenvalues for only few cases even after reducing the number of times in the do loop. Any hint would be highly...

Would anyone please suggest the best numerical method for integrating complicated integrands and singular functions? I find that Gauss quadratures are suggested by many references but I can't find a method to obtain the abscissas and weights for a general function.