Skip to main content
added 163 characters in body
Source Link
Ulrich Neumann
  • 56.8k
  • 2
  • 26
  • 60

Change the integration variable x to xdh=x/h , in this way parameter h dissappears

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

pic=ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

List of solutions found:

pic[[1, 1]][[1]]
(*{{0.8, 0.0416416},{0.804065,0.0422764}, 
{0.807732,0.0428571}, {0.8375,0.0477974}, 
{0.838872, 0.0480182},...}*)

Alternatively try direct solution NSolve

sol[Vi_] := {Vi, b} /. NSolve[{0 ==eq[[1]], 0 < b < 1}, b][[1]]   
Table[sol[Vi], {Vi, .8, 5, .1}]  

Change the integration variable x to xdh=x/h , in this way parameter h dissappears

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

pic=ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

List of solutions found:

pic[[1, 1]][[1]]
(*{{0.8, 0.0416416},{0.804065,0.0422764}, 
{0.807732,0.0428571}, {0.8375,0.0477974}, 
{0.838872, 0.0480182},...}*)

Change the integration variable x to xdh=x/h , in this way parameter h dissappears

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

pic=ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

List of solutions found:

pic[[1, 1]][[1]]
(*{{0.8, 0.0416416},{0.804065,0.0422764}, 
{0.807732,0.0428571}, {0.8375,0.0477974}, 
{0.838872, 0.0480182},...}*)

Alternatively try direct solution NSolve

sol[Vi_] := {Vi, b} /. NSolve[{0 ==eq[[1]], 0 < b < 1}, b][[1]]   
Table[sol[Vi], {Vi, .8, 5, .1}]  
added 186 characters in body
Source Link
Ulrich Neumann
  • 56.8k
  • 2
  • 26
  • 60

Change the integration variable x to xdh=x/h , in this way parameter h dissappears

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

ContourPlot[0pic=ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

List of solutions found:

pic[[1, 1]][[1]]
(*{{0.8, 0.0416416},{0.804065,0.0422764}, 
{0.807732,0.0428571}, {0.8375,0.0477974}, 
{0.838872, 0.0480182},...}*)

Change the integration variable x to xdh=x/h

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

Change the integration variable x to xdh=x/h , in this way parameter h dissappears

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

pic=ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here

List of solutions found:

pic[[1, 1]][[1]]
(*{{0.8, 0.0416416},{0.804065,0.0422764}, 
{0.807732,0.0428571}, {0.8375,0.0477974}, 
{0.838872, 0.0480182},...}*)
Source Link
Ulrich Neumann
  • 56.8k
  • 2
  • 26
  • 60

Change the integration variable x to xdh=x/h

eq = ArcTan[1/d*Sqrt[(a + b)/(1 - b)]] - (Vi/h)*Integrate[Sqrt[(1 + (xdh)^2)^(-1)], {xdh,(1/2)*Log[b], 0}] h + (Pi/4) +m*Pi

ConturPlot shows possible real solutions {Vi,b}

ContourPlot[0 == eq, {Vi, .8, 5}, {b, 0.0, .6}, FrameLabel -> {Vi, b}]

enter image description here