Plotting an infinite integral
4 次查看(过去 30 天)
显示 更早的评论
How do I plot an infinite integral? I am trying to plot a solution to Laplaces equation where -inf < x < inf.
What I have tried is:
f = 100;
b = 50;
x = linspace(-100, 100, .1);
y = linspace(0, b, .01);
k = (0:.1:100);
kk = (-100:.2:100);
[K, KK] = meshgrid(k, kk);
T = zeros(length(x), length(y));
T = T + trapz(K, trapz(KK, sinh(K.*(b - y))./sinh(K.*b)*f.*...
cos(K.*(KK - x))/pi));
imagesc(x, y, flipud(T))
0 个评论
回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!