I am trying to solve the following integral

beta^alpha*x^(alpha-1)*exp(-beta*x)/((1+x)*GAMMA(alpha)) 

with respect to x, between zero and infinity.

If i just integrate it the result that I obtain is

(alpha^2*beta^(alpha+1)*exp(beta)*GAMMA(-alpha-1, beta)+alpha*beta^(alpha+1)*exp(beta)*GAMMA(-alpha-1, beta)+beta-alpha)/beta

My problem is that the alpha is positive and consequently the GAMMA functions do not exist, since the argument...

Optimal_Portfolio_Ch.mw

I was using maple as a documentation software and was preparing some notes for my exam. I had a sudden computer shutdown and the file didn't saved properly.

The file now opens but most of the content is gone. Even the one which was properly saved.

I see Joe Riel 

Hello,

I need to solve bond orders in Maple the formula is

Bond order = [number of bonding electrons – number of anti-bonding electrons]/2

I need to implement it in Maple, I dk how

Thanks

when calculate the intersection of I2 and J2, i exactly follow this picture
https://skydrive.live.com/redir?resid=E0ED7271C68BE47C!311

however, it is G list which is not the stated in this picture

restart;
with(Groebner):
I2 := {x^2+y^3-1, x-y*x+3};
J2 := {y*x^2-1};
G1 := Basis({expand(w*I2[1]),expand(w*I2[2]), expand((1-w)*J2[1])}, 'tord', deglex(x,y));
G2 := Basis({expand(w*I2[1]),expand(w*I2[2]), expand((1-w)*J2[1])}, 'tord', deglex(w,x,y));

Dear All

How do I solve the followinh set of non-linear equations for the unknowns as given below:
    

    Zth1=R1(1-e^(-t1/tau1))+R2(1-e^(-t1/tau2))+R3(1-e^(-t1/tau3))+R4(1-e^(-t1/tau4))+R5(1-e^(-t1/tau5))+R6(1-e^(-t1/tau6))

    Zth2=R1(1-e^(-t2/tau1))+R2(1-e^(-t2/tau2))+R3(1-e^(-t2/tau3))+R4(1-e^(-t2/tau4))+R5(1-e^(-t2/tau5))+R6(1-e^(-t2/tau6))

    Zth3=R1(1-e^(-t3/tau1...

I'm trying to graph a contour map of this function: e^(y/x)

>with(plots)

>contourplot(e^(y/x),x=-3..3,y=-0.5..0.5)

I thought the contour map for this function would looks more like several lines intersects at the origin. But when I graph it on Maple, it doesn't looks like the graph I'm looking for.

Did I put in a wrong range?

Hello everybody,

I' m trying to run completely The optmization template from a batch file (.bat) and wonder how to :

1.  Make the buttons working automatically.
2.  What  command should be used in a batch file to execute utterly a maple program
3. What kind of script file can be perfect for Matlab / Maple.

The point is to run in parallel Matlab and Maple, using a script file(.bat) or another external script format.

Hi everyone, I'm new to the forum and also new to Maple.

I am trying to do the following:

"Write a procedure, marks, to take a student’s mark and output the appropriate grade, based on the following table:"

Hello,

I am trying to make a three dimensional cube in maple by using matrices. I also have to make it rotate about the x,y,z axes, project it with three different origins, and translate it (I'm assuming this means translate it across an axis?). I have already made the cube, but when I try to rotate it across one of the axes by multiplying my matrix by one of the rotation matrices (this is what I put un maple: P1:= matrix([[1,0,0,0],[0,cos(theta),-sin(theta),0],[0,sin(theta),cos(theta),0...

Hi.

I need to plot some function f(x,y) on a right hexagon pllaced on (x,y)-plane. I can define my hexagon in following way:

if -(4*pi)/3/sqrt(3)<x<(2*pi)/3/sqrt(3) then -sqrt(3)*(x-4*pi/3/sqrt(3))<y<sqrt(3)*(x-4*pi/3/sqrt(3))

if -(2*pi)/3/sqrt(3)<x<(2*pi)/3/sqrt(3) then -2*pi/3<y<2*pi/3

if (2*pi)/3/sqrt(3)<x<(4*pi)/3/sqrt(3) then -sqrt(3)*(4*pi/3/sqrt(3)-x)<y<sqrt(3)*(4*pi/3/sqrt(3)-x)

My haxagon looks like

hey friends..
Am trying to solve 2d poission equation with finite difference method and using different scheme is got

u_i-1,j -8u_i,j + u_i+1,j + 3u_i,j-1 +3u_i,j+1 = -h *16

and i have made code for this which is

> restart;
> with(LinearAlgebra);
> printf(`The lefthand endpoint on the X-axis is 0.\n`);
> OK := FALSE;
> while OK = FALSE do printf(`Input...

OK, I want to plot 10-15 curves on the same coordinate system P. I generate and put the equations of said curves in an array with a procedure, but when plotting the array, I receive an array of plots, instead of the 10 curves on one coordinate system. Any ideas how can I do that?

The obvious solution, something like P:=plot (something) Q:=plot(other thing) display({P,Q}) obviously doesn't work for me.

the goal is write a for loop to iterate and get the power of a variable with degree function

i use nops and convert(F, list) can access each terms and variables

however y^2 can be convert to list [y,2] 

hope to add if condition to stop convert this.

how to do better?

i just guess...

r-x^4
u-(x^3)*y
v-x*y^3
w-y^4

Step 1
with(Groebner):
K := {r-x^4,u-(x^3)*y,v-x*y^3,w-y^4};
G := Basis(K, 'tord', degrevlex(r,u,v,w));

Step 2
G intersect k[r,u,v,w]

G2 := Basis({w*G[1],w*G[2],w*G[3],w*G[4],w*G[5],w*G[6],w*G[7],w*G[8],w*G[9],w*G[10],w*G[11],w*G[12],w*G[13],w*G[14], (1-w)*K[1], (1-w)*K[2], (1-w)*K[3], (1-w)*K[4]}, 'tord', degrevlex(r,u,v,w));

G2 has too much though include the correct solution

correct solution...