Skip to main content

MINIMUM OF GAMMA FUNCTION

## This string finds the minimum of the Gamma function
## integral(x^(alpha-1)*exp(-x))dx) (x=[0,infinity])
## in the interval 0<alpha<4.
dx=1E-2; ## Increment in x
x=0:dx:100; ## x array
Gamma=[]; ## Empty Gamma function
alpha=0.00:dx:4.00; ## The independent variable of Gamma function
for i=1:length(alpha)
## Call the trapezoid function to
## calculate all entries of Gamma function
## corresponding to entries of alpha
f=x.^(alpha(i)-1).*exp(-x); ## integrand of Gamma function
Gamma(i)=trapezoid(x,f);
endfor
plot(alpha,Gamma);
title('Gamma values vs alpha');
legend('Gamma(alpha)');
xlabel('alpha');
ylabel('Gamma');
## The following 'for loop' finds the minimum value of Gamma
## and the corresponding alpha value.
minimum=Gamma(1);
for n=1:length(Gamma)
if(minimum>Gamma(n))
minimum=Gamma(n);
n; ## The array index where the minimum of Gamma is.
m=alpha(n); ## The alpha value corresponding to the min of Gamma
endif
endfor;
printf('Minimum of the Gamma Function in the interval 0<alpha<4 is %f\n for the alpha=%f\n' ,minimum,m)
print('-dpsc',MINGAMMA.ps');
save -text MINGAMMA.dat

Comments

Post a Comment

Popular posts from this blog

FACTORIAL

## Function that calculates the factorial of a number ## Usage : f=factorial(n) function f=factorial(n) ## Initialize the output f=1; ## Check whether the input is correct if ( (n<0) || (rem(n,1)~=0) ) printf("n cannot be a negative number. Exiting...\n"); return endif for num=1:n f*=num; endfor endfunction
Post a Comment
Read more

One Dimensional Harmonic Oscillator-Numerov Method

x=[]; h0=1; M=4; N=M+1; x(1)=0; x(N)=x(1)+h0*M; x=x(1):h0:x(N) A=zeros(N); A(1,1)=-2*(5*(x(N)*h0)^2/12+1); A(N,N)=-2*(5*(x(1)*h0)^2/12+1); A(1,2)=1-(x(M)*h0)^2)/12; A(N,M)=1-(x(2)*h0)^2)/12; B=zeros(N); B(1,1)=B(N,N)=-10*(h0^2)/6; B(1,2)=B(N,M)=-(h0^2)/6; for i=2:M B(i,i)=-10*(h0^2)/6; B(i,i-1)=B(i,i+1)=-(h0^2)/6; A(i,i)=-2*(5*(x(N+1-i)*h0)^2/12+1); A(i,i+1)=1-(x(N-i)*h0)^2)/12; A(i,i-1)=1-(x(N+2-i)*h0)^2)/12; end A B
Post a Comment
Read more

REALISTIC AND NON-REALISTIC STRINGS

## A 'realistic', 'non-elastic' string, which responses to any ## bending and has stifness. This script takes in the previous ## and the present profiles and iterates to find ## the profile in the next time step. The ratio 'r' is not 1 ## like in the 'non-realistic' string since the speed of the wave ## always less than the speed of the string it should be less than 1 ## for best and most stable solution ##constants dx=1e-2 ## Spatial increment (m) L=2 ## Length of the string (m) M=L/dx ## Dimensionless partition E=1e-4 ## Dimensionless stiffnes function ynext=propagate_stiff(ynow,yprev,r) ## Quick and dirty way to fix boundary conditions -- for each step ## they are the same as the previous step. ynext=ynow; ynow(1)=ynow(2)=0; ##Entering the loop for i=3:length(ynow)-1 ## boundary condition ynow(length(ynow)-1)=ynow(length(ynow)-2)=0; ## Divide the ynext with many terms into three parts for easiness ynext(i)=(2−(2*r^2)−(6*E*(r^2)*...
Post a Comment
Read more