evaluating double integral iteratively
1 次查看(过去 30 天)
显示 更早的评论
I am trying to discretize a kernel K(x,t) for solving an integral equation. I need to implement
for i=1...N
for j=1..N
a(i,j) = Int1 Int2 K(x,t)dx dt
where Int1 has limits i/N to (i+1)/N and Int2 has limits j/N to (j+1)/N
I tried the symbolic computation but it takes a lot of time. I want to use the numerical computation and speed it up. thanks
0 个评论
回答(1 个)
Star Strider
2014-7-14
If K(x,t) is a function, you can use integral2.
2 个评论
Star Strider
2014-7-14
My pleasure!
That’s my understanding of how integral2 works. See the More About in the documentation, and its friends (linked to at the end of the page). Note that only integral can take array-valued functions, so if K(x,t) in one such, you’ll have to nest calls to integral. If it’s scalar-valued, the integral functions will be happy and will return a scalar.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!