       integration - Interpolating Function

• To: mathgroup at smc.vnet.net
• Subject: [mg92015] integration - Interpolating Function
• From: marion70 at gmx.de
• Date: Wed, 17 Sep 2008 04:28:38 -0400 (EDT)

```Hi,

I have two different problems concerning Mathematica:
First: to calculate a characteristic function  \[CurlyPhi][u_] for
each u I have to evaluate a solution of a system of ordinary
differential equations depending on u:

solution[u_] := NDSolve[
{z1'[t] == -\[Kappa] \[Theta]y y1[
t] - (ly \[Mu] (y1[
t] - \[Mu]  y1[t]^2 - \[Mu] y2[t]^2))/((1 - \[Mu] y1[
t])^2 + \[Mu]^2 y2[t]^2),
z2'[t] == -\[Kappa] \[Theta]y y2[
t] - (ly \[Mu]  y2[
t])/((1 - \[Mu] y1[t])^2 + \[Mu]^2 y2[t]^2),
y1'[t] == \[Kappa] y1[t] - (\[Sigma]^2)/2 (y1[t]^2 - y2[t]^2),
y2'[t] == \[Kappa] y2[t] - \[Sigma]^2 y1[t] y2[t] - u,
z1[T] == 0,
z2[T] == 0,
y1[T] == 0,
y2[T] == 0}, {z1, z2, y1, y2}, {t, 0, T}]
a1[u_] := z1 /. solution[u][[1, 1]]
a2[u_] := z2 /. solution[u][[1, 2]]
b1[u_] := y1 /. solution[u][[1, 3]]
b2[u_] := y2 /. solution[u][[1, 4]]
\[CurlyPhi][u_] := E^(a1[u] + I a2[u] + b1[u] + I b2[u] x0)

as a result I get a set of Interpolating Functions, so
I can evaluate the characteristic functions in specific points u. Next
I have to integrate over this characteristic function:

Needs["FourierSeries`"];
f[z_]:=NInverseFourierTransform[\[CurlyPhi][u_]  ,
u,z,FourierParameters=84_{1,1}]

but Mathematica cannot solve this problem.

Does anybody know an answer or a way to fix this problem?

Second: Can anybody give me advice, how to implement efficiently a
recursion, that depends on two variables? (Does Mathematica have
another function besides RSolve which solves this kind of problems?)

I would be really grateful if someone could help me.

Best regards,
Marion

```

• Prev by Date: Re: imagesize to full width of current screen
• Next by Date: Re: Functional programming?
• Previous by thread: bug... GraphicsGrid ItemAspectRatio Spacings
• Next by thread: Variable names in Mathematica