Hi Users!

Hope everyone in fine and enjoying good health. I am facing problem to differential the following expression with respect to first variable (mentioned as red). Please help me to fix this query

f(z*sqrt(a/nu[f]), U*t/(2*x))

Thanks in advance

Dear Users!

Hope you would be fine with everying. I want to solve the following 2nd order linear differential equation. 

(1+B)*(diff(theta(eta), eta, eta))+C*A*(diff(theta(eta), eta)) = 0;
where A is given as

A := -(alpha*exp(-sqrt((omega+1)*omega*(M^2+alpha+1))*eta/(omega+1))*omega+alpha*exp(-sqrt((omega+1)*omega*(M^2+alpha+1))*eta/(omega+1))+exp(-sqrt((omega+1)*omega*(M^2+alpha+1))*eta/(omega+1))*omega-alpha*omega+exp(-sqrt((omega+1)*omega*(M^2+alpha+1))*eta/(omega+1))-alpha-omega-1)/sqrt((omega+1)*omega*(M^2+alpha+1));
I want solution for any values of omega, alpha, M, B, C and L. The BCs are below:

BCs := (D(theta))(0) = -1, theta(L) = 0.

I am waiting your response, 

Dear Users!

I want to simple the Gamma function occures in a matrix. I need the simpliest form of this matrix. If there is some thing common in all entries take it common. Thanks

Matrix(6, 6, {(1, 1) = 2*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(alpha*(2*alpha+3)*(1+2*alpha)*GAMMA(alpha)*GAMMA(alpha+1/2)^2), (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (2, 1) = (1/2)*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(alpha*(2*alpha+3)*(1+2*alpha)*GAMMA(alpha)*GAMMA(alpha+1/2)^2), (2, 2) = (1/4)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(GAMMA(alpha+1/2)^3*alpha^3*(alpha+1)^3*GAMMA(alpha)^3)*Pi^(1/4)/((2*alpha+3)*(1+2*alpha)*GAMMA(alpha+1/2)^2*alpha^2*(alpha+1)^2*GAMMA(alpha)^2), (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (3, 1) = (1/16)*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(GAMMA(alpha+1/2)^2*alpha*(alpha+1)*(1+2*alpha)*GAMMA(alpha)), (3, 2) = (1/8)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(GAMMA(alpha+1/2)^3*alpha^3*(alpha+1)^3*GAMMA(alpha)^3)*Pi^(1/4)/((2*alpha+3)*(1+2*alpha)*GAMMA(alpha+1/2)^2*alpha^2*(alpha+1)^2*GAMMA(alpha)^2), (3, 3) = (1/32)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(alpha^3*GAMMA(alpha+3/2)^3*(2+alpha)^3*GAMMA(alpha)^3)*Pi^(1/4)/((alpha+1)*(2*alpha+3)*alpha^2*GAMMA(alpha+3/2)^2*(2+alpha)^2*GAMMA(alpha)^2), (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 2*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(alpha*(2*alpha+3)*(1+2*alpha)*GAMMA(alpha)*GAMMA(alpha+1/2)^2), (4, 5) = 0, (4, 6) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = (3/2)*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(alpha*(2*alpha+3)*(1+2*alpha)*GAMMA(alpha)*GAMMA(alpha+1/2)^2), (5, 5) = (1/4)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(GAMMA(alpha+1/2)^3*alpha^3*(alpha+1)^3*GAMMA(alpha)^3)*Pi^(1/4)/((2*alpha+3)*(1+2*alpha)*GAMMA(alpha+1/2)^2*alpha^2*(alpha+1)^2*GAMMA(alpha)^2), (5, 6) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = (1/16)*(18*alpha+19)*GAMMA(alpha+5/2)*Pi^(1/4)*sqrt(GAMMA(alpha+1)*GAMMA(alpha+1/2)^3)/(alpha*(alpha+1)*(2*alpha+3)*(1+2*alpha)*GAMMA(alpha)*GAMMA(alpha+1/2)^2), (6, 5) = (3/8)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(GAMMA(alpha+1/2)^3*alpha^3*(alpha+1)^3*GAMMA(alpha)^3)*Pi^(1/4)/((2*alpha+3)*(1+2*alpha)*GAMMA(alpha+1/2)^2*alpha^2*(alpha+1)^2*GAMMA(alpha)^2), (6, 6) = (1/32)*GAMMA(alpha+5/2)*sqrt(2)*sqrt(alpha^3*GAMMA(alpha+3/2)^3*(2+alpha)^3*GAMMA(alpha)^3)*Pi^(1/4)/((alpha+1)*(2*alpha+3)*alpha^2*GAMMA(alpha+3/2)^2*(2+alpha)^2*GAMMA(alpha)^2)})

Dear Users!

Hope you would be fine. In the following maple code, I want to write the derivative of psi in term of psi like it did manually in red portion. For higher M and k it very hard to do it manully. It there any command to fix my problem for any value of k and M.

restart; k := 2; M := 4;

with(linalg); with(LinearAlgebra);

printlevel := 2;

for i while i <= 2^(k-1) do

for j from 0 while j <= M-1 do

psi[M*i+j-M+1] := simplify(2^((1/2)*k)*sqrt(GAMMA(j+1)*(j+alpha)*GAMMA(alpha)^2/(Pi*2^(1-2*alpha)*GAMMA(j+2*alpha)))*(sum((-1)^i1*GAMMA(j-i1+alpha)*(2*(2^k*x-2*i1+1))^(j-2*i1)/(GAMMA(alpha)*factorial(i1)*factorial(j-2*i1)), i1 = 0 .. floor((1/2)*j))));

`&psi;&psi;`[M*i+j-M+1] := simplify(diff(psi[M*i+j-M+1], x))

end do

end do; r := 2^(k-1)*M;

VV := Vector[column](r, proc (i) options operator, arrow; psi[i] end proc);

DV := Vector[column](r, proc (i) options operator, arrow; `&psi;&psi;`[i] end proc);

``&psi;&psi;`[2] := 8*sqrt((alpha+1)*(1/2))*sqrt(2)*sqrt(alpha*GAMMA(alpha)^2*4^alpha/GAMMA(2*alpha))/sqrt(Pi) = 8*sqrt((alpha+1)*(1/2))*psi[1];

`&psi;&psi;`[3] := 16*sqrt((2+alpha)*(alpha+1)/(1+2*alpha))/sqrt(2)*(2*sqrt(2)*sqrt((alpha+1)*GAMMA(alpha)^2*4^alpha/GAMMA(1+2*alpha))*alpha*(4*x+1)/sqrt(Pi)) = 16*sqrt((2+alpha)*(alpha+1)/(1+2*alpha))*psi[2]/sqrt(2)

I am waiting your response. Thanks

Dear Users! 

Hope you would be fine with everything. I want to define a matrix F of M+1 by M+1 order having element of the following form:

I derived the F[r,s] but confuse who to generate matrix now.

restart; M := 5; printlevel := 3; for r from 2 while r <= M+1 do for s while s <= r-1 do if type(r+s, odd) then F[r, s] := 2^(k+1)*sqrt((2*r-1)*(2*s-1)) end if end do end do