Mathematica: Convolution Integral

  • #1
Hi all! I'm new to Mathematica.
I have written a code for performing a convolution integral (as follows) but it seems to be giving out error messages:

My code is:

a[x_?NumericQ] := PDF[NormalDistribution[40, 2], x]
b[k_?NumericQ, x_?NumericQ] := 0.0026*Sin[1.27*k/x]^2
c[k_?NumericQ, x_?NumericQ] := {a[x]*b[k, x]}
d[k_?NumericQ] := NIntegrate[c[k, x], {x, 0, Infinity}]
Plot [d[k], {k, 0, 350}]

It gives the following error message multiple times:

ntegrand c[0.00715,x] is not numerical at {x} = {124.67}


And the plot dosen't show up..!
Does anyone have any suggestions on how to fix it ?
 
Last edited:

Answers and Replies

  • #2
281
138
The argument of the sin function is wrong. It must be dimentionless, i.e. k multiplied by x not divided by x.
 

Related Threads on Mathematica: Convolution Integral

Replies
9
Views
2K
Replies
3
Views
3K
Replies
3
Views
3K
Replies
1
Views
2K
Replies
2
Views
1K
Replies
7
Views
3K
Replies
1
Views
3K
Replies
2
Views
9K
Top