
Re: Double Integral error
Posted:
Nov 28, 2012 12:44 PM


Hi Steve,
collecting all the data I have, I came up with the idea that dblquad does not work cause k1 is a huge array and not just a scalar value.
Hence, I should call dblquad for each k1(i) with i = 1,2,...,numel(k1).
Unfortunately I have no clue how to perform this.
Furthermore, I was thinking of a nested quad integration:
phi_33 = @(k2,k3) E_k0./((4*pi.*k.^4)).*(k1.^2 + k2.^2);
for i = 2:numel(k1)
PHI_33_k3(i) = quad(phi_33,k3(i1),k3(i));
PHI_33 = quad(PHI_k3,k2(i1),k2(i)); end
but,again, I don't know how to perform all this (taking into account k1(i) too).
I really look forward to hearing from you again.
BR, Francesco

