Computation of complex valued integral











up vote
0
down vote

favorite












I would like to get a simple formula for the complex valued integral
$$
frac{1}{sqrt{2pisigma_Z^2}}int_{-infty}^inftyexpleft(iux+delta d(t)vert xvertright)expleft(-frac{(x-mu_Z)^2}{2sigma_Z^2}right) dx.
$$

This formula is also allowed to contain the cdf $Phi$ of a standard normal RV. The above integral appears in an ODE I am trying to solve and I would like to avoid solving an integral in each step of my ODE solver. Hence, I would like to have a simple formula for the integral only containing functions that are available in most standard programming libraries.
If I replace $iu$ with $u$, i.e. only consider the real valued case I can simplify the above integral to
$$
-expleft((u+delta d(t))mu_Z+(u+delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(frac{(u+delta d(t))sigma_Z^2+mu_Z}{sigma_Z}right)\
-expleft((u-delta d(t))mu_Z+(u-delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(-frac{(u-delta d(t))sigma_Z^2 + mu_Z}{sigma_Z}right)
$$



However, in the complex case I would need to evaluate the cdf $Phi$ at a complex number which is implemented in most computer packages like MATLAB. Is there a way to get around this?










share|cite|improve this question
























  • $sigma=sigma_Z , ?$. Also what is $delta d(t)$?
    – Diger
    Nov 23 at 23:50










  • Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
    – lbf_1994
    Nov 23 at 23:56












  • And $delta$ is another constant?
    – Diger
    Nov 23 at 23:57










  • yes exactly, it is also a constant
    – lbf_1994
    Nov 23 at 23:59










  • And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
    – Diger
    Nov 24 at 0:00















up vote
0
down vote

favorite












I would like to get a simple formula for the complex valued integral
$$
frac{1}{sqrt{2pisigma_Z^2}}int_{-infty}^inftyexpleft(iux+delta d(t)vert xvertright)expleft(-frac{(x-mu_Z)^2}{2sigma_Z^2}right) dx.
$$

This formula is also allowed to contain the cdf $Phi$ of a standard normal RV. The above integral appears in an ODE I am trying to solve and I would like to avoid solving an integral in each step of my ODE solver. Hence, I would like to have a simple formula for the integral only containing functions that are available in most standard programming libraries.
If I replace $iu$ with $u$, i.e. only consider the real valued case I can simplify the above integral to
$$
-expleft((u+delta d(t))mu_Z+(u+delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(frac{(u+delta d(t))sigma_Z^2+mu_Z}{sigma_Z}right)\
-expleft((u-delta d(t))mu_Z+(u-delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(-frac{(u-delta d(t))sigma_Z^2 + mu_Z}{sigma_Z}right)
$$



However, in the complex case I would need to evaluate the cdf $Phi$ at a complex number which is implemented in most computer packages like MATLAB. Is there a way to get around this?










share|cite|improve this question
























  • $sigma=sigma_Z , ?$. Also what is $delta d(t)$?
    – Diger
    Nov 23 at 23:50










  • Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
    – lbf_1994
    Nov 23 at 23:56












  • And $delta$ is another constant?
    – Diger
    Nov 23 at 23:57










  • yes exactly, it is also a constant
    – lbf_1994
    Nov 23 at 23:59










  • And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
    – Diger
    Nov 24 at 0:00













up vote
0
down vote

favorite









up vote
0
down vote

favorite











I would like to get a simple formula for the complex valued integral
$$
frac{1}{sqrt{2pisigma_Z^2}}int_{-infty}^inftyexpleft(iux+delta d(t)vert xvertright)expleft(-frac{(x-mu_Z)^2}{2sigma_Z^2}right) dx.
$$

This formula is also allowed to contain the cdf $Phi$ of a standard normal RV. The above integral appears in an ODE I am trying to solve and I would like to avoid solving an integral in each step of my ODE solver. Hence, I would like to have a simple formula for the integral only containing functions that are available in most standard programming libraries.
If I replace $iu$ with $u$, i.e. only consider the real valued case I can simplify the above integral to
$$
-expleft((u+delta d(t))mu_Z+(u+delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(frac{(u+delta d(t))sigma_Z^2+mu_Z}{sigma_Z}right)\
-expleft((u-delta d(t))mu_Z+(u-delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(-frac{(u-delta d(t))sigma_Z^2 + mu_Z}{sigma_Z}right)
$$



However, in the complex case I would need to evaluate the cdf $Phi$ at a complex number which is implemented in most computer packages like MATLAB. Is there a way to get around this?










share|cite|improve this question















I would like to get a simple formula for the complex valued integral
$$
frac{1}{sqrt{2pisigma_Z^2}}int_{-infty}^inftyexpleft(iux+delta d(t)vert xvertright)expleft(-frac{(x-mu_Z)^2}{2sigma_Z^2}right) dx.
$$

This formula is also allowed to contain the cdf $Phi$ of a standard normal RV. The above integral appears in an ODE I am trying to solve and I would like to avoid solving an integral in each step of my ODE solver. Hence, I would like to have a simple formula for the integral only containing functions that are available in most standard programming libraries.
If I replace $iu$ with $u$, i.e. only consider the real valued case I can simplify the above integral to
$$
-expleft((u+delta d(t))mu_Z+(u+delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(frac{(u+delta d(t))sigma_Z^2+mu_Z}{sigma_Z}right)\
-expleft((u-delta d(t))mu_Z+(u-delta d(t))^2frac{sigma_Z^2}{2}right)Phileft(-frac{(u-delta d(t))sigma_Z^2 + mu_Z}{sigma_Z}right)
$$



However, in the complex case I would need to evaluate the cdf $Phi$ at a complex number which is implemented in most computer packages like MATLAB. Is there a way to get around this?







integration complex-analysis analysis normal-distribution complex-integration






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Nov 23 at 23:55

























asked Nov 23 at 23:43









lbf_1994

10910




10910












  • $sigma=sigma_Z , ?$. Also what is $delta d(t)$?
    – Diger
    Nov 23 at 23:50










  • Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
    – lbf_1994
    Nov 23 at 23:56












  • And $delta$ is another constant?
    – Diger
    Nov 23 at 23:57










  • yes exactly, it is also a constant
    – lbf_1994
    Nov 23 at 23:59










  • And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
    – Diger
    Nov 24 at 0:00


















  • $sigma=sigma_Z , ?$. Also what is $delta d(t)$?
    – Diger
    Nov 23 at 23:50










  • Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
    – lbf_1994
    Nov 23 at 23:56












  • And $delta$ is another constant?
    – Diger
    Nov 23 at 23:57










  • yes exactly, it is also a constant
    – lbf_1994
    Nov 23 at 23:59










  • And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
    – Diger
    Nov 24 at 0:00
















$sigma=sigma_Z , ?$. Also what is $delta d(t)$?
– Diger
Nov 23 at 23:50




$sigma=sigma_Z , ?$. Also what is $delta d(t)$?
– Diger
Nov 23 at 23:50












Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
– lbf_1994
Nov 23 at 23:56






Yes, $sigma_Z=sigma$ (I've edited the question). $d(t)$ is the unknown function in the ODE I am trying to solve. You can just consider it to be a constant when computing the integral. The above integral is one summand in the right hand side of an ODE for $d(t)$.
– lbf_1994
Nov 23 at 23:56














And $delta$ is another constant?
– Diger
Nov 23 at 23:57




And $delta$ is another constant?
– Diger
Nov 23 at 23:57












yes exactly, it is also a constant
– lbf_1994
Nov 23 at 23:59




yes exactly, it is also a constant
– lbf_1994
Nov 23 at 23:59












And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
– Diger
Nov 24 at 0:00




And what do you expect from an answer? If you chuck the integral into Maple or Mathematica you will get an exponential and an error-function. Is that what you want?
– Diger
Nov 24 at 0:00















active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3010980%2fcomputation-of-complex-valued-integral%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown






























active

oldest

votes













active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Mathematics Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


Use MathJax to format equations. MathJax reference.


To learn more, see our tips on writing great answers.





Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


Please pay close attention to the following guidance:


  • Please be sure to answer the question. Provide details and share your research!

But avoid



  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.


To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3010980%2fcomputation-of-complex-valued-integral%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Wiesbaden

Marschland

Dieringhausen