Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi guys,

I'm trying to apply boundary condition on the seris i generated using Adomian decomposition. Could some one please check it for me?
 

NULL

u[0] := a1+a2*y:

NULL

NULL

w[0] := a3+a4*y:

theta[0] := a6*y+a5

phi[0] := a8*y+a7

NULL

NULL

``

``

``

``

``

A[1] := R*(diff(u[0], y))+A-Gr*(B*phi[0]+theta[0])/Ree+Ha^2*(alpha*u[0]+beta*w[0])/(alpha^2+beta^2)

R*a2+A-Gr*(a5+a6*y+B*(a7+a8*y))/Ree+Ha^2*(alpha*(a1+a2*y)+beta*(a3+a4*y))/(alpha^2+beta^2)

(1)

u[1] := int(A[1], y = 0 .. y)

(1/2)*(-Gr*(a6+B*a8)/Ree+Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^2+R*a2*y+A*y-Gr*(a5+B*a7)*y/Ree+Ha^2*(alpha*a1+beta*a3)*y/(alpha^2+beta^2)

(2)

u[11] := int(u[1], y = 0 .. y)

(1/3)*(-(1/2)*Gr*(a6+B*a8)/Ree+(1/2)*Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^3+(1/2)*(R*a2+A-Gr*(a5+B*a7)/Ree+Ha^2*(alpha*a1+beta*a3)/(alpha^2+beta^2))*y^2

(3)

u = u[0]+u[11]

u = a1+a2*y+(1/3)*(-(1/2)*Gr*(a6+B*a8)/Ree+(1/2)*Ha^2*(alpha*a2+beta*a4)/(alpha^2+beta^2))*y^3+(1/2)*(R*a2+A-Gr*(a5+B*a7)/Ree+Ha^2*(alpha*a1+beta*a3)/(alpha^2+beta^2))*y^2

(4)

A[2] := R*(diff(w[0], y))-Ha^2*(beta*u[0]-alpha*w[0])/(alpha^2+beta^2)

R*a4-Ha^2*(beta*(a1+a2*y)-alpha*(a3+a4*y))/(alpha^2+beta^2)

(5)

w[1] := int(A[2], y = 0 .. y)

-(1/2)*Ha^2*(beta*a2-alpha*a4)*y^2/(alpha^2+beta^2)+R*a4*y-Ha^2*(beta*a1-alpha*a3)*y/(alpha^2+beta^2)

(6)

w[11] := int(w[1], y = 0 .. y)

-(1/6)*Ha^2*(beta*a2-alpha*a4)*y^3/(alpha^2+beta^2)+(1/2)*(R*a4-Ha^2*(beta*a1-alpha*a3)/(alpha^2+beta^2))*y^2

(7)

``

w = w[0]+w[11]

w = a3+a4*y-(1/6)*Ha^2*(beta*a2-alpha*a4)*y^3/(alpha^2+beta^2)+(1/2)*(R*a4-Ha^2*(beta*a1-alpha*a3)/(alpha^2+beta^2))*y^2

(8)

A[3] := R*Pr*(diff(theta[0], y))-2*Br*((diff(u[0], y))^2+(diff(w[0], y))^2+M^2*C*(u[0]^2+w[0]^2))-gamma*R*Pr*theta[0]

R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*((a1+a2*y)^2+(a3+a4*y)^2))-gamma*R*Pr*(a5+a6*y)

(9)

theta[1] := int(A[3], y = 0 .. y)

-(2/3)*Br*M^2*C*(a2^2+a4^2)*y^3+(1/2)*(-2*Br*M^2*C*(2*a1*a2+2*a3*a4)-gamma*R*Pr*a6)*y^2+R*Pr*a6*y-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))*y-gamma*R*Pr*a5*y

(10)

theta[11] := int(theta[1], y = 0 .. y)

-(1/6)*Br*M^2*C*(a2^2+a4^2)*y^4+(1/3)*(-Br*M^2*C*(2*a1*a2+2*a3*a4)-(1/2)*gamma*R*Pr*a6)*y^3+(1/2)*(R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))-gamma*R*Pr*a5)*y^2

(11)

theta = theta[0]+theta[11]

theta = a5+a6*y-(1/6)*Br*M^2*C*(a2^2+a4^2)*y^4+(1/3)*(-Br*M^2*C*(2*a1*a2+2*a3*a4)-(1/2)*gamma*R*Pr*a6)*y^3+(1/2)*(R*Pr*a6-2*Br*(a2^2+a4^2+M^2*C*(a1^2+a3^2))-gamma*R*Pr*a5)*y^2

(12)

``

 

NULL

A[4] := R*Sc*(diff(phi[0], y))-K*Sc*phi[0]

R*Sc*a8-K*Sc*(a7+a8*y)

(13)

phi[1] := int(A[4], y = 0 .. y)

-(1/2)*K*Sc*a8*y^2+R*Sc*a8*y-K*Sc*a7*y

(14)

phi[11] := int(phi[1], y = 0 .. y)

-(1/6)*K*Sc*a8*y^3+(1/2)*(R*Sc*a8-K*Sc*a7)*y^2

(15)

phi = phi[0]+phi[11]

phi = a7+a8*y-(1/6)*K*Sc*a8*y^3+(1/2)*(R*Sc*a8-K*Sc*a7)*y^2

(16)

``


 

Download second_problem.mw

Hi everyone,

Please, I really need your expertise advice(s) on what i am not doing right in the code attached below. I was actually writing the code on Multi-step DTM, but, instead of continuing from the last point, it starting all over again. Please, your expertise will save a soul here. Regard

Thanks

MsDTM3.mw

I'm pretty annoyed with maple because when I try to create proc's it dooes not use tabbing or smart tabbing. This makes writing readable code a chore. Basically it's always screwing up alignment. (it wants to left justify things a lot)...

 

What's worse is that using the tab key takes one outside the proc rather than adding a tab, making one having to use spaces.

Any way to fix this?

Hello, I try this integral, but Maple not solutions, the answer is the same integral.

int(N*exp((2*(-(1/2)*x^2-2*a^2*ln(x)*x/(-2*x^2+2)))/sigma^2)/(-2*x^2+2), x = -1 .. 1)

N,a, sigma are constants.

Regards

I'm having trouble connecting from Maple on Windows 10 to MSQL Server. I tried Microsoft recommended drivers such as sqljdbc_6.4.0.0, did (as I thought) all required steps. The only invariable result I get is "Cannot load driver". I was wandering if anyone had implemented such a construction. Driver name & version , connection string and Java version would be greatly appreciated. Another option is to have any driver, which connects to any of standard databases (Oracle, MySQL).  The only limitation is- it must be from Windows 7 or 10.
          Thanks.
           A.B.

I remember there was a command that used to start an independent java session each time. I thought it's "xmaple -singlemode" but it's not. So could anybody please remind me what was the command. I just killed a second java session after losing kernel connection and java running wild. I'd be very happy to do my calculation in another java session so that if anything goes bad II won't have to open all my 5 tabs again.

Integrating a positive definite function (normal distribution) and a >= 0 function (Heaviside) should not return a negative value.

 

 

with(Statistics);
X := RandomVariable(Normal(1, sqrt(2.25)));
int(PDF(X, x)*Heaviside(x^7-5*x^4-3*x+1), x = -infinity .. infinity);
                         -0.08507120131

 

 

Bug.mw

I'm using dsolve command to solve a differential equation. Using infolevel to 3 will tell me the classification of said DE. However, how can I see the step by step solution? I'm using Maple as a study tool so I do solve manually a DE then I'd like to compare my answer with Maple's. How can I acomplish this? Thanks in advance. 

Please may I know if you can offer ma student's discount ob the seleted version.

Thank you.

Fred.

I open a completely new document, type in a few simple expressions and get this weird complex number out. Can someone pleas help me fix this?

Is there any way we can get maple to tell us which variable we are recursively assigning? 

solve({-infinity < a , a < -1, -1 < b ,  b < 0});

 

f(x) := piecewise(0 < x, x^(3/2)*sin(1/x), x = 0, 0, undefined);
plot(f(x));

gives me the following error:

Error, (in plot) incorrect first argument piecewise(0 < x, (HFloat(2.739493386336394e-116)+HFloat(2.739493386336394e-116)*I)*x^(3/2), x = 0, 0, undefined)

I just want to see the function plot. With Wolfram Alpha this is no deal at all!
 

Maple doesn't completely check the condition on the number of trials "n" for Binomial and NegativeBinomial distributions (package Statistics).
The attribute "Conditions" explicitely says that n must be a strictly positive integer but no strictly positive real valuereturna an error (ok, it would be stupid to set n to a non integer value !!!).

I think it is a default that ought to be corrected in future releases (this default still exists in Maple 2018)

 

restart

kernelopts(version)

`Maple 2015.2, APPLE UNIVERSAL OSX, Dec 20 2015, Build ID 1097895`

(1)

with(Statistics):


BINOMIAL DISTRIBUTION

X := RandomVariable(Binomial(n, p)):
L := [attributes(X)][3]:
A := exports(L)

Conditions, ParentName, Parameters, CDF, CharacteristicFunction, Kurtosis, Mean, Median, Mode, MGF, ProbabilityFunction, Skewness, Support, Variance, VariationCoefficient, CDFNumeric, QuantileNumeric, RandomSample, RandomSampleSetup, RandomVariate

(2)

L:-Conditions

[0 <= p, p <= 1, n::posint]

(3)

# Maple should return an error for N is not of type posint
#
# It seems that Sample uses floor(N)

N := 10.49; type(N::posint);
P := 1/2:
X := RandomVariable(Binomial(N, P)):
Mean(X), N*P;
ProbabilityFunction(X, k);
S := Sample(X, 10^6):
Mean(S);


# A non consistent result (only non negative values of k should be accepted)

eval(ProbabilityFunction(X, k), k=evalf(Pi));

N := 10.49

 

false

 

5.245000000, 5.245000000

 

piecewise(k < 0, 0, binomial(10.49, k)*(1/2)^k*(1/2)^(10.49-k))

 

HFloat(4.998903)

 

.1096019539

(4)


NEGATIVE BINOMIAL DISTRIBUTION

X := RandomVariable(NegativeBinomial(n, p)):
L := [attributes(X)][3]:
A := exports(L):
L:-Conditions

[0 < n, 0 < p, p <= 1]

(5)

N := 10.49:
P := 1/2:
X := RandomVariable(NegativeBinomial(N, P)):

Mean(X)


 

Download BinomialLaw.mw

 

First 77 78 79 80 81 82 83 Last Page 79 of 2097