Advanced Search

Dkunb

70 Reputation

5 Badges

3 years, 204 days
Contact Dkunb

MaplePrimes Activity


These are questions asked by Dkunb

How do I integrate in Maple?...

Dkunb 70 Maple 2022
integration
April 30 2024
2 0

How can I  get the result for the integration when m =n or m is not equal to n and How can I add assumption that m, or n can be even or odd?

Thanks in advance for your help.

 
> 

restart;
phi := proc(k,x,L)
  if (type(k,even)) then sqrt(2)*sin(Pi*k*x/L)/sqrt(L)
  else sqrt(2)*cos(Pi*k*x/L)/sqrt(L)
  end if;
end proc:
> 

Int(phi(m,x,L)*_h^2/m2*diff(phi(n,x,L),x,x),x=-L/2..L/2);

Int(-2*cos(Pi*m*x/L)*_h^2*Pi^2*n^2*cos(Pi*n*x/L)/(L^3*m2), x = -(1/2)*L .. (1/2)*L)

 (1)
 

 

Download test1.mw

 

Factor in matrix form...

Dkunb 70 Maple 2022
matrix factorization simplification
April 16 2024
2 3

ras5_v1.mw 

How can I  change eq13 to eq14 without using op command manually?

Thank you for your help in advance,

Integration of Hermite to Laguerre...

Dkunb 70 Maple 2022
February 27 2024
1 3

How  to make the integration 

 

into the following form:



Thanks in advance,Q1.mw

How to find the normalization factor ?...

Dkunb 70 Maple
February 17 2024
2 0

SolnforOS_dsolve.mw 
I would like to solve eq6 in terms of sigma and n. I would appreciate for your help.

How to find the solutions matrix(array)...

Dkunb 70 Maple 2021
array ode syntax
June 13 2023
1 11

How to find Z(t) matrix(array) =[ [z[1,0](t),z[1,1](t)],[z[2,0](t),z[2,1](t)]]

> 

restart;
> 

with(IntegrationTools):with(Physics):
> 

 

> 

 

> 

 

> 

h1 := (m,n)->(n+1/2)*KroneckerDelta[n,m];

proc (m, n) options operator, arrow; Physics:-`*`(n+Physics:-`*`(1, Physics:-`^`(2, -1)), Physics:-KroneckerDelta[n, m]) end proc

 (1)
> 

h2 := (mu,nu,m2,l)->-(nu*Pi/l)^2/(2*m2)*KroneckerDelta[mu,nu];

proc (mu, nu, m2, l) options operator, arrow; Physics:-`*`(-1, Physics:-`*`(Physics:-`^`(Physics:-`*`(nu, Pi, Physics:-`^`(l, -1)), 2), Physics:-`^`(Physics:-`*`(2, m2), -1), Physics:-KroneckerDelta[mu, nu])) end proc

 (2)
> 

v1 := (m,n,m1)->sqrt(min(n,m)!/max(n,m)!)*(2*m1)^(-abs(n-m)/2)*exp(-1/(4*m1))*LaguerreL(min(n,m),abs(n-m),-1/(2*m1));

proc (m, n, m1) options operator, arrow; Physics:-`*`(sqrt(Physics:-`*`(factorial(min(n, m)), Physics:-`^`(factorial(max(n, m)), -1))), Physics:-`^`(Physics:-`*`(2, m1), Physics:-`*`(-1, Physics:-`*`(abs(n-m), Physics:-`^`(2, -1)))), exp(Physics:-`*`(-1, Physics:-`*`(1, Physics:-`^`(Physics:-`*`(4, m1), -1)))), LaguerreL(min(n, m), abs(n-m), Physics:-`*`(-1, Physics:-`*`(1, Physics:-`^`(Physics:-`*`(2, m1), -1))))) end proc

 (3)
> 

v2 := (mu,nu,l)->4*Pi^2*l*mu*nu*(exp(l/2)-(-1)^(mu+nu)*exp(-l/2))/((Pi*(mu+nu))^2+l^2)/((Pi*(mu-nu))^2+l^2);

proc (mu, nu, l) options operator, arrow; Physics:-`*`(4, Physics:-`^`(Pi, 2), l, mu, nu, exp(Physics:-`*`(l, Physics:-`^`(2, -1)))-Physics:-`*`(Physics:-`^`(-1, mu+nu), exp(Physics:-`*`(-1, Physics:-`*`(l, Physics:-`^`(2, -1))))), Physics:-`^`(Physics:-`^`(Physics:-`*`(Pi, mu+nu), 2)+Physics:-`^`(l, 2), -1), Physics:-`^`(Physics:-`^`(Physics:-`*`(Pi, mu-nu), 2)+Physics:-`^`(l, 2), -1)) end proc

 (4)
> 

h:=(m,n,mu,nu,m1,m2,l)->evalf(h1(m,n)+h2(mu,nu,m2,l)+v1(m,n,m1)+v2(mu,nu,l));

proc (m, n, mu, nu, m1, m2, l) options operator, arrow; evalf(h1(m, n)+h2(mu, nu, m2, l)+v1(m, n, m1)+v2(mu, nu, l)) end proc

 (5)
> 

m1:=1:m2:=1:l:=1:
> 

H:= (m,n,mu,nu)->h(m,n,mu,nu,m1,m2,l);

proc (m, n, mu, nu) options operator, arrow; h(m, n, mu, nu, m1, m2, l) end proc

 (6)
> 

H(1,1,0,0);

2.668201175

 (7)
> 

eq1:= diff(z(m,nu,t),t)=-I*Sum(Sum(H(m,n,mu,nu)*z(n,mu,t),n=1..N),mu=0..M);

diff(z(m, nu, t), t) = -I*(Sum(Sum(((n+.5000000000)*Physics:-KroneckerDelta[m, n]-4.934802202*nu^2*Physics:-KroneckerDelta[mu, nu]+.7788007831*(factorial(min(m, n))/factorial(max(m, n)))^(1/2)*2.^(-.5000000000*abs(-1.*n+m))*LaguerreL(min(m, n), abs(-1.*n+m), -.5000000000)+39.47841762*mu*nu*(1.648721271-.6065306597*(-1.)^(mu+nu))/((9.869604404*(mu+nu)^2+1.)*(9.869604404*(mu-1.*nu)^2+1.)))*z(n, mu, t), n = 1 .. N), mu = 0 .. M))

 (8)
> 

 

> 

zint := Array([[1,0],[0,1]]);

Matrix(2, 2, {(1, 1) = 1, (1, 2) = 0, (2, 1) = 0, (2, 2) = 1})

 (9)
> 

Z:= Array(1..2,1..2);
 for i from 1to 2 do
    for j from 1 to 2 do
        Z[i,j]:= dsolve({eq1,zint[i,j]},numeric,output=listprocedure);
    end do;
end do;
> 

 

Matrix(2, 2, {(1, 1) = 0, (1, 2) = 0, (2, 1) = 0, (2, 2) = 0})

  

Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations

  
> 

 

Download test3.mw

1 2 3 4 5 6 Page 2 of 6