Mathematica 9 is now available
Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2012

[Date Index] [Thread Index] [Author Index]

Search the Archive

numerical solution to a set of equations with hypergeometric function

  • To: mathgroup at smc.vnet.net
  • Subject: [mg126120] numerical solution to a set of equations with hypergeometric function
  • From: Joseph Döner <doenermachtscheisse at googlemail.com>
  • Date: Thu, 19 Apr 2012 03:54:04 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com

Hi,

i ran into some trouble solving the following set of equations
http://imageshack.us/photo/my-images/337/gleichungen.jpg/ for N and M,
all other variables are known. I think the problem lies within the F-
function. It's supposed to be the hypergeometric function.
I used FindRoot and the hypergeometric2f1regularized command at first,
but all i got were some errors. When i tried to solve it with the
integral representation and the power series of the function I still
didn't get any results. If i insert all the values it does look like
this:

FindRoot[{((0.5*0.5 + 0.5*0.5)*(2.5*Hypergeometric2F1Regularized[1,
1.894427191, 2.894427191, -G/M] +
1.180339887*Hypergeometric2F1Regularized[1, 0.894427191,1.894427191, -
M/G]))/(1 - (0.5 - G)*(2.5*Hypergeometric2F1Regularized[1,
1.894427191, 2.894427191, -G/M] +
1.180339887*Hypergeometric2F1Regularized[1, 0.894427191,  1.894427191,
-M/G]) - (0.5 - M)*(1.180339887*Hypergeometric2F1Regularized[1,
1.894427191, 2.894427191, -M/G] -
(-2.5)*Hypergeometric2F1Regularized[1, 0.894427191, 1.894427191, -G/
M])) == (((1 - 0.5)*0.5 + (1 -
0.5)*0.5)*(50*Hypergeometric2F1Regularized[1, 1.2, 2.2, -(1 - G)/(1 -
M)] + 8.3333333*Hypergeometric2F1Regularized[1, 0.2, 1.2, -(1 - M)/(1
- G)]))/(1 - (G - 0.5)*(50*Hypergeometric2F1Regularized[1, 1.2, 2.2, -
(1 - G)/(1 - M)] + 8.3333333*Hypergeometric2F1Regularized[1, 0.2, 1.2,
-(1 - M)/(1 - G)]) - (M -
0.5)*(8.333333*Hypergeometric2F1Regularized[1, 1.2, 2.2, -(1 - M)/(1 -
G)] - (-50)*Hypergeometric2F1Regularized[1, 0.2, 1.2, -(1 - G)/(1 -
M)])), ((0.5*0.5 +
0.5*0.5)*(1.180339887*Hypergeometric2F1Regularized[1, 1.894427191,
2.894427191, -G/M] + 8.33333*Hypergeometric2F1Regularized[1, 0.2, 1.2,
-M/G]))/(1 - (0.5 -
G)*(1.180339887*0.5/0.5*Hypergeometric2F1Regularized[1, 1.1894427191,
2.894427191, -M/G] - (-2.5)*Hypergeometric2F1Regularized[1,
0.894427191, 1.894427191, -G/M])) == (((1 - 0.5)*0.5 + (1 -
0.5)*0.5)*(8.3333333*Hypergeometric2F1Regularized[1, 0.2, 1.2, (1 - M)/
(1 - G)]))/(1 - (G -  0.5)*(50*Hypergeometric2F1Regularized[1, 1.2,
2.2, -(1 - G)/(1 - M)] + 8.3333333*Hypergeometric2F1[1, 0.2, 1.2, -(1
- M)/(1 - G)]) - (M - 0.5)*(8.333333*Hypergeometric2F1Regularized[1,
1.2, 2.2, -(1 - M)/(1 - G)] - (-50)*Hypergeometric2F1Regularized[1,
0.2, 1.2, -(1 - G)/(1 - M)]))}, {{G, 0.5}, {M, 0.5}}]

If you have any idea why it does not work, let me now :)


Many thanks in advance!



  • Prev by Date: Re: NIntegrate about singular point
  • Next by Date: Re: convolution involving UnitStep
  • Previous by thread: Re: Extract coefficients of a trig polynomial
  • Next by thread: Encode and Get