| Author |
Comment/Response |
Antonio
|
 05/05/12 12:00pm
Hello, I am new to mathematica, I have started couple of days ago, but I succeded in doing lot of interesing thigs which I had problems to do with other software.
However I am still stuck on programming and variable values assignement.
I have a small piece of code
alfa = 10.; beta = 10.; i = 3.;
(*Evaluation of polynomial of order i*)
df = Expand[
Pochhammer[(alfa + 1), i]/i!*
Hypergeometric2F1[-i, i + alfa + beta + 1, alfa + 1, 0.5 - x/2]];
(*Evaluation of first derivative of polynomial of order i*)
df3 = 1/N[Expand[D[df, x]^2]];
(*Evaluation of the roots of polynomial order i*)
root = Solve[
Pochhammer[(alfa + 1), i]/i!*
Hypergeometric2F1[-i, i + alfa + beta + 1, alfa + 1,
0.5 - x/2] == 0, x];
(*Save in a the values of the roots of the polynomial*)
a = x /. root[[ ;; ]];
(*Defintion of the integrand - This doesn't work, it should make a \
vector integ with the different values a, but doesn't work!*)
integ = (Gamma[alfa + beta + 2]*(1 - x)^
alfa*(1 + x)^
beta/(2^(alfa + beta + 1)*Gamma[alfa + 1]*
Gamma[beta + 1]))*((df)/(x - a))^2;
(*Numerical integration of the integral - This doesn't work since the \
integrand is not properly done*)
integr = NIntegrate[integ, {x, -1, 1}];
(*Give in df3Num in place of x the value of a - This doesn't work I \
guess I made a mistake in defining *)
df3Num = df3 /. a
(*Evaluation of the weight as product of two vector of numbers*)
weight = df3Num*integr
I am using Mathematica 8.
The next step will be the make the code at different values of i, at the moment I am using only one i.
However, I am stuck in giving the variables the values, it would be appreaciated some help on this.
Thanks a lot
Antonio
URL: , |
|
