       NIntegrate within NDSolve-optimization?

• To: mathgroup at smc.vnet.net
• Subject: [mg128169] NIntegrate within NDSolve-optimization?
• From: maves <vesel.mateja at gmail.com>
• Date: Thu, 20 Sep 2012 00:25:26 -0400 (EDT)
• Delivered-to: l-mathgroup@mail-archive0.wolfram.com
• Delivered-to: l-mathgroup@wolfram.com
• Delivered-to: mathgroup-newout@smc.vnet.net
• Delivered-to: mathgroup-newsend@smc.vnet.net

```Hello,
I am trying to numerically solve a set of differential equations, and the problem is that the expression also contains numerical integrals, I believe this is why I am getting large numerical errors.

Firstly, I define three integrals (for instance for b=6, c=0.5):

int1[x_?NumericQ, z_?NumericQ] :=
NIntegrate[-(R (x - R Cos[\[Phi]]))/
((x - R Cos[\[Phi]])^2 + (R Sin[\[Phi]])^2 + (z -
Z)^2)^(3/2), {\[Phi], 0, 2 Pi}, {R, b, b + c}, {Z,
R - (b + c), -R + (b + c)}];

int2[x_?NumericQ, z_?NumericQ] :=
NIntegrate[-(R (z-Z))/
((x - R Cos[\[Phi]])^2 + (R Sin[\[Phi]])^2 + (z -
Z)^2)^(3/2), {\[Phi], 0, 2 Pi}, {R, b, b + c}, {Z,
R - (b + c), -R + (b + c)}];

int3[x_?NumericQ, z_?NumericQ] :=
NIntegrate[-R/((x - R Cos[\[Phi]])^2 + (R Sin[\[Phi]])^2 + (z -
Z)^2)^(1/2), {\[Phi], 0, 2 Pi}, {R, b, b + c}, {Z,
R - (b + c), -R + (b + c)}];

And then I use them in NDSolve:

s = NDSolve[{
x''[t] == int1[x[t], z[t]],
z''[t] == int2[x[t], z[t]],
x == 7.1, x' == 0, z' == 0, z' == 1}, {x, x=
', z,
z'}, {t, 0, 5}][]

And plot the solution:
ParametricPlot[Evaluate[{x[t], z[t]} /. s], {t, 0, 5},
PlotRange -> All]

Plot[Evaluate[0.5*((x'[t])^2+(z'[t])^2+int3[x[t],z[t]]) /. s], {t, 0, 5}]

I change the initial conditions and try to find something "useful", but this second plot shows me, that something just isn't right, as it should be constant, but it oscillates with a huge amplitude. Could this be due to numerical errors and how you think I should optimize my code to get a good result? Is this combination of NIntegrate and NDSolve even acceptable?

Thank you.

But the solution just doesn't seem right

```

• Prev by Date: Re: Changing typeface using the find and replace dialog
• Next by Date: About CDF
• Previous by thread: Solved: Weighting in simultaneous NonlinearModelFits