Advanced Search

Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015
View unanswered posts
Maple 2015 Questions and Posts Feed

Question about Igusa invariant ...

Asked by: LeeHoYeung 535
invariant
June 04 2016
0  3

how to compute Igusa invariant in maple?

why it is not included in maple?

Roots of the characteristic polynomial ...

Asked by: iman 35
dsolve differential-equations characteristic_polynomial
June 02 2016
0  1 
how i can calculate roots of the characteristic polynomial equations {dsys and dsys2}
and dsolve them with arbitrary initial condition for differennt amont of m and n?
thanks
Kr.mw

restart; a := 1; b := 2; Number := 10; q := 1; omega := 0.2e-1
``

Q1 := besselj(0, xi*b)*(eval(diff(bessely(0, xi*r), r), r = a))-(eval(diff(besselj(0, xi*r), r), r = a))*bessely(0, xi*b):

J := 0:

m := 0:

U1 := (int(r*K1[m]*(diff(K_01[m], r)+K_01[m]/r), r = a .. b))/(int(r*K1[m]^2, r = a .. b)); -1; U2 := -(int(r*K_01[m]*(diff(K1[m], r)), r = a .. b))/(int(r*K_01[m]^2, r = a .. b)); -1; U3 := (int(r^2*omega^2*K_01[m], r = a .. b))/(int(r*K_01[m]^2, r = a .. b))

0.6222222222e-3/K_01[12]

 (1)

Q2 := besselj(1, eta*b)*(eval(diff(bessely(1, eta*r), r), r = a))-(eval(diff(besselj(1, eta*r), r), r = a))*bessely(1, eta*b):

E2 := unapply(Q2, eta):

m := 0:

 
dsys := {diff(S_mn(t), t, t, t)+xi[m]^2*(diff(S_mn(t), t, t))+(-U1*U2+eta__n^2)*(diff(S_mn(t), t))+xi[m]^2*eta__n^2*S_mn(t) = -(2*U2*b_m/(Pi*xi[m])*(-besselj(0, xi[m]*b)/besselj(1, xi[m]*a)))*q+xi[m]^2*U3}; 1; dsolve(dsys)

{S_mn(t) = (3111111111/5000000000000)/(K_01[12]*eta__n^2)+_C1*cos(eta__n*t)+_C2*sin(eta__n*t)+_C3*exp(-xi[12]^2*t)}

 (2)

dsys2 := {diff(Q_mn(t), t, t, t)+xi[m]^2*(diff(Q_mn(t), t, t))+(-U1*U2+eta__n^2)*(diff(Q_mn(t), t))+xi[m]^2*eta__n^2*Q_mn(t) = -2*besselj(0, xi[m]*b)*U1*U2*b_m*(1-exp(-xi[m]^2*t))/(besselj(1, xi[m]*a)*Pi*xi[m]^3)}; 1; dsolve(dsys2)

{Q_mn(t) = _C1*exp(-xi[12]^2*t)+_C2*sin(eta__n*t)+_C3*cos(eta__n*t)}

 (3)

``

 

``

Download Kr.mw

 

A problem in using Grid[Launch]...

Asked by: sand15 867
parallel grid
June 02 2016
0  10

Hi all,

I start working with the Grid package.
To familiarize myself with it I ran the "primeChecker" example, which of course has worked perfectly well.

Next I did this (a priori harmless) simple modifications :

  1. within the primeChecker procedure :
    replace  myVal := userData[thisNode+1] :
    by         myVal := userData[thisNode+1, 1] :
  2. before launching the procedure with Grid[Launch] :
    replace  userData := [ .... ] :  #which is a list
    by         userData := [ .... ] :  # the same thing
                 userData := convert(userData, matrix): 


I get the following error message :
error, (in unknown) Matrix index out of range


What does it mean and how can I fix this ?

Thanks in advance

Integral of complex function...

Asked by: Alger 524
integration
June 01 2016
1  8

Dear all,

I have somme difficulties to calculate this integral.

int(I*sqrt((R*exp(I*theta)+1)/(R*exp(I*theta)-a)), theta = 0 .. Pi);

Thanks

Parametric solution of polynomial system...

Asked by: iman 35
solve parameter
May 31 2016
0  12

hi.

how i can solve two equation with respect to parameter sigma1

SOLL.mw

> 

restart; pprime11 := -16395.36603*q1+5.811117425*q1*sigma1^2+3526.724044*p1-1.250000000*p1*sigma1^2+4.999870968*10^11*p1^3+4.999870970*10^11*p1*q1^2+7.967307034*10^14*p1^2*q1+4.999870966*10^12*sigma1*p1^2*q1-2.655769012*10^14*q1^3+4.999870968*10^12*sigma1*q1^3-17633.62022*q1*sigma1+6.250000000*q1*sigma1^3
> 

qprime11 := 2.655769012*10^14*p1^3-7.967307034*10^14*p1*q1^2+4.999870970*10^11*p1^2*q1-4.999870968*10^12*sigma1*p1^3-4.999870966*10^12*p1*sigma1*q1^2+3526.724044*q1-1.250000000*q1*sigma1^2+16395.36603*p1-5.811117425*p1*sigma1^2+4.999870968*10^11*q1^3+17633.62022*p1*sigma1-6.250000000*p1*sigma1^3:
> 

-50 < sigma1 and sigma1 < 50:
> 

sigma1 <> 53.11665685, -53.11665685:
> 

SOLL := solve({pprime11, qprime11}, real)

Warning,  computation interrupted

  
> 

``

 

Download SOLL.mw

Combining two plots...

Asked by: redgazelle 20
statistics plot3d
May 31 2016
4  2

Hello,

does anyone know a way to combine two plots in one where one is created with ScatterPlot3D (Package:Statistics) and the other one with plot3D (Package:plots)?

 

Normally you would write something like this:

but that only works if the plots are from the same package...

Problem with solution of linear algebraic system ...

Asked by: Muhammad Usman 240
solve linear-algebra
May 30 2016
0  11

Hello dear!

Hope everyone fine, I am facing to solve the attached problem please find and fix it. I am waiting your kind response. 

naveed.mw

Composition of functions...

Asked by: MrNsh 10
composition
May 30 2016
2  7

I have entered three functions into Maple and I would like to create a set of all possible two-function and three-function compositions involving the three functions. For instance, for my three functions f(x), g(x) and h(x), the set would contain f(g(x)), g(f(x)), f(g(h(x))), f(f(f(x))), etc. 

I'm also looking for a method that will be generalisable to larger numbers of initial functions than just three.

Labeling solutions...

Asked by: MrNsh 10
solve label
May 30 2016
0  3

I have used solve to find the solution to an equation that has two solutions, and I want to give each solution a label so that I can use each individually in subsequent manipulations. How do I label each solution separately?

How to find this integral?...

Asked by: taro 505
integration int
May 29 2016
0  6

Hello people in mapleprimes,

 

I cannot obtain a proper result from the following code.

a:=int(((beta/beta[1,2])^(-theta/(1-theta))-kappa[1]^(-theta/(1-theta)))*m*beta^(m-1),beta=0 .. kappa[1]*beta[1,2]);

 

Please tell me if you know how to have maple calculate it.

 

Thanks in advance.

 

taro

Want to solve system of linear equations...

Asked by: Muhammad Usman 240
linear
May 28 2016
0  1

Hello Dear!

I want to solve the system of linear equation but facing some problem please see the attachmen. I am waiting your positive response 

1_(1).mw

How to solve polynomial system?...

Asked by: farahnaz 15
polynomial solve
May 27 2016
0  1

how i can solve equations and gain{p1,p2,q1,q2}?

thanks

solve.mw

Q1 := -mu1*p1-upsilon1*q1+gamma11*q1*(p1^2+q1^2)+gamma12*q1*(p2^2+q2^2)-delta1*(2*p1*q1*p2-q2*(p1^2-q1^2))-s11*F*q1+s12*F*q2 = 0

-mu1*p1-upsilon1*q1+gamma11*q1*(p1^2+q1^2)+gamma12*q1*(p2^2+q2^2)-delta1*(2*p1*q1*p2-q2*(p1^2-q1^2))-s11*F*q1+s12*F*q2 = 0

 (1)

Q2 := -mu1*q1+upsilon1*p1-gamma11*p1*(p1^2+q1^2)-gamma12*p1*(p2^2+q2^2)-delta1*(2*p1*q1*q2+p2*(p1^2-q1^2))-s11*F*p1-s12*F*p2 = 0

-mu1*q1+upsilon1*p1-gamma11*p1*(p1^2+q1^2)-gamma12*p1*(p2^2+q2^2)-delta1*(2*p1*q1*q2+p2*(p1^2-q1^2))-s11*F*p1-s12*F*p2 = 0

 (2)

Q3 := -mu2*p2-upsilon2*q2+gamma21*q2*(p1^2+q1^2)+gamma22*q2*(p2^2+q2^2)+delta2*(3*p1^2*q1-q1^3)+s21*F*q1 = 0

-mu2*p2-upsilon2*q2+gamma21*q2*(p1^2+q1^2)+gamma22*q2*(p2^2+q2^2)+delta2*(3*p1^2*q1-q1^3)+s21*F*q1 = 0

 (3)

Q4 := -mu2*q2+upsilon2*p2-gamma21*p2*(p1^2+q1^2)-gamma22*p2*(p2^2+q2^2)+delta2*(-p1^3+3*p1*q1^2)-s21*F*p1 = 0

-mu2*q2+upsilon2*p2-gamma21*p2*(p1^2+q1^2)-gamma22*p2*(p2^2+q2^2)+delta2*(-p1^3+3*p1*q1^2)-s21*F*p1 = 0

 (4)

SOLL := solve({Q1, Q2, Q3, Q4}, real)

Error, (in PolynomialSystemSolvers:-PseudoResultant) numeric exception: division by zero

  

SOLL := solve({Q1, Q2, Q3, Q4}, {p1, p2, q1, q2})

Warning,  computation interrupted

  

``

 

Download solve.mw

 

How to differentiate with respect to y?...

Asked by: iman 35
differentiation
May 27 2016
1  1

how i can differential with respect to ''y'' variable

thanks...help_(2).mw

T(r, theta) := .5*alpha+sum([gamma*n*(r^n+R^(2*n)+r^(-n))*sin(n*theta)], n = 1 .. infinity)

.5*alpha+[-gamma*r*(r^2-1)*sin(theta)/(4*cos(theta)^2*r^2-4*cos(theta)*r^3+r^4-4*cos(theta)*r+2*r^2+1)-((1/2)*I)*R^2*gamma*(exp(I*theta)*R^4-exp((3*I)*theta)*R^4-exp(I*theta)+exp((3*I)*theta))/((R^2*exp(I*theta)-1)^2*(-R^2+exp(I*theta))^2)+((1/2)*I)*r*gamma*(exp(I*theta)*r^2-exp((3*I)*theta)*r^2-exp(I*theta)+exp((3*I)*theta))/((exp(I*theta)-r)^2*(exp(I*theta)*r-1)^2)]

 (1)

``

y := r, proc (theta) options operator, arrow; r*sin(theta) end proc

r, proc (theta) options operator, arrow; r*sin(theta) end proc

 (2)

 

 

diff(T, y)

Error, invalid input: diff received proc (theta) options operator, arrow; r*sin(theta) end proc, which is not valid for its 3rd argument

  

``

 

Download help_(2).mw

How to export undirected graph to C#?...

Asked by: longqnh 10
May 26 2016
2  2

Hi all,

I drew a undirected graph using Maple

restart;
with(GraphTheory);
with(SpecialGraphs);
with(RandomGraphs);
G := Graph(undirected, {[{1, 2}, 4], {1, 4}, {2, 3}, {2, 5}, {3, 4}});
DrawGraph(G, style = circle);

I want to import this graph to picturebox or panel on C#, please help me!

Replica of MATLAB code in Maple?...

Asked by: iman 35
May 22 2016
0  0

hi

if possible to convert matlab file in to maple fie program??convert_to_maple_program.txt

thanks

 

First 44 45 46 47 48 49 50 Last Page 46 of 73