The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Double Integral error
Replies: 5   Last Post: Nov 28, 2012 1:24 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Francesco Perrone

Posts: 39
Registered: 5/2/12
Re: Double Integral error
Posted: Nov 28, 2012 12:44 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

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(i-1),k3(i));

PHI_33 = quad(PHI_k3,k2(i-1),k2(i));

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.


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.