QuantumKnoll
- 1
- 0
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 doesn't show up..!
Does anyone have any suggestions on how to fix it ?
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 doesn't show up..!
Does anyone have any suggestions on how to fix it ?
Last edited: