MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

I try to make a question with an equation with numbered variables. This works fine when evaluating:

y=x1/x2

I need to have the lefthand variable indexed also like this:

y2=x1/x2

When entering y2=x1/x2 as the answer Maple TA won't evaluate it as a correct answer!?

I have the following code:

restart;
PDE := diff(u(x, t), t) = diff(u(x, t), x, x)-sin(x+t)+cos(x+t);
IBC:= D[1](u)(0,t)=-sin(t),
D[1](u)(1,t)=-(u(1, t))^4+(cos(1+t))^4-sin(1+t),
u(x,0)=cos(x);

pds := pdsolve( PDE, [IBC], numeric, time = t, range = 0 .. 1,

spacestep = 0.1e-1, timestep = 0.1e-1,
errorest=true
)

 

And I want to plot the difference |pds(t,x) - cos(x+t)| in maple for x=1 and t=0..5

 

I thought to use the following piece of commands but I get an error:

P:=unapply(pds, t,x);
Q:=abs(P(t,x)-cos(x+t));
Q:-plot(x=1, t=0..1);

I get an error that Q isn't a module.

I thought that unapply is used for this case, can you help me with this simple task?

 

Thanks in advance.

 

Hi

Im trying to make a small maple document, where i can calculate the moecule mass of any molecule. I want to assign all atoms with a value for its atom mass, but it will not let me unassign I, as this is the standard unit for imaginary numbers. Is there any way i can unassign I and give it a new value just in the current instance of the program? I dont want it to unassign I in other documents.

Thank you in advance

Hi all

I have this repetition with for:

for i from 1 to 10000 do x[i]:=sqrt(i); print(i); end do:

I want to know at the same time when the repetition is working, the value of i where the calculation is. Not at the last time of the calculation of all the repititions.

Example: when i=1, I want to see 1, when i=600, I want to see 600 at the same time of execution the repetition and not the last of the repetition  

Thank you

 

f := x^4-c1*x^3+c2*x^2-c3*x+c4;

restart;
ferrai := -x1^3*x2*x3*x4-x1^2*x2^2*x3^2-x1^2*x2^2*x4^2-x1^2*x3^2*x4^2-x1*x2^3*x3*x4-x1*x2*x3^3*x4-x1*x2*x3*x4^3-x2^2*x3^2*x4^2+x1^2*x2*x3*y+x1^2*x2*x4*y+x1^2*x3*x4*y+x1*x2^2*x3*y+x1*x2^2*x4*y+x1*x2*x3^2*y+x1*x2*x4^2*y+x1*x3^2*x4*y+x1*x3*x4^2*y+x2^2*x3*x4*y+x2*x3^2*x4*y+x2*x3*x4^2*y-x1*x2*y^2-x1*x3*y^2-x1*x4*y^2-x2*x3*y^2-x2*x4*y^2-x3*x4*y^2+y^3;
coeff(ferrai, y^3);
coeff(ferrai, y^2);
coeff(ferrai, y);
res := ferrai - coeff(ferrai, y^3)*y^3 - coeff(ferrai, y^2)*y^2 - coeff(ferrai, y)*y;

c2 := -coeff(ferrai, y^2)/coeff(ferrai, y^3);
sys1 := c1*c3 - 4*c4 = coeff(ferrai, y)/coeff(ferrai, y^3);
sys2 := -c3^2-(c1^2)*c4+4*c2*c4 = res;
solve([sys1, sys2],[c1,c3,c4]);

number of equations is not enough, is it possible to find back c1,c3,c4?

though c2 is easy to know

assume not starting from resolvent of quartic,

how to find [[1,2],[3,4]] and [[1,3],[2,4]] and [[1,4],[2,3]] ?

a1*a2 + a3*a4 if [[1,3,2,4],[1,2]]

a1*a3 + a2*a4 if [2,3][[1,3,2,4],[1,2]]

a1*a4 + a2*a3 if [2,4],[[1,3,2,4],[1,2]]

moreover, i use multiplication of permutation group can not multiply below

though i know how to operate multiplication by hand, i follow the syntax to do multiplication, seems not the way

with(group):
mulperms([[2,3]], [[1,3,2,4],[1,2]]);
mulperms([[1,3,2,4],[1,2]], [[2,3]]);

mulperms([2,3], [[1,3,2,4],[1,2]]);
mulperms([[1,3,2,4],[1,2]], [2,3]);

 

how to find [[1,3,2,4],[1,2]] and [2,3][[1,3,2,4],[1,2]] and [2,4],[[1,3,2,4],[1,2]]?

and why Orbits(G) put permutation group into power, how to display the result of this power group?

i use this or quartic polynomials' resolvent

with(GroupTheory):
G := SpecialUnitaryGroup(3, 1);
Orbits(G);
Orbit(1,G);
Orbit(3,G);

So im trying to solve multiple ODES using dsolve(numeric) but I jsut cant get it to work.

I kep getting this one error: 

Error, (in f) unable to store '-HFloat(0.020918994979034728)-HFloat(0.09184018333019917)*I' when datatype=float[8]

 

This is my uploaded file 

(for some reason the uploaded file didnt show the error at the bottom so i just pasted it in the spot it would appear.

Download CHE504_HW5.mw

Given:

restart

dt := 2.07*0.254e-1

0.52578e-1

(1)

dto := 2.38*0.254e-1

0.60452e-1

(2)

ho := 5000

5000

(3)

``

Ltube := 5

5

(4)

``

Po := 2.5

2.5

(5)

To := 350

350

(6)

dp := 0.3e-2

0.3e-2

(7)

tau := 3

3

(8)

dpore := 5.2*10^(-9)

0.5200000000e-8

(9)

kfluid := 0.485e-1

0.485e-1

(10)

ksolid := 1.67

1.67

(11)

kpipe := 17

17

(12)

Cpair := 1.01

1.01

(13)

`μ_air` := 3.16*10^(-5)

0.3160000000e-4

(14)

rho := P(z)*(28.9*(1/1000))/(Rgas*(T(z)+273.15))

0.2890000000e-1*P(z)/(Rgas*(T(z)+273.15))

(15)

``

Rgas := 8.2057*10^(-5)

0.8205700000e-4

(16)

``

``

density calculation n/v = (P/RT)*mw_air, kg/m^3

``

`ρi` := (2.5*28.966)/(0.82057e-1*(350+273.15))

1.416186012

(17)

NULL

Volumetric flow, m^3/s

((1/60)*((1/300)*((1/1000)*(180*1000000)*.4535)*(1/24))*(1/60))/(1.416)

0.2224085844e-2

(18)

``

NULL

Superficial velocity, m/s

NULL

Vsi := evalf(%/((1/4)*Pi*(2.07*0.254e-1)^2))

1.024362190

(19)

Gs constant, kg/m^2-s

Error, missing operator or `;`

 

Gsi := `ρi`*Vsi

1.450687405

(20)

Gs := 1.45

1.45

(21)

``

``

NULL

Void fraction, epsilon,b

`εb` := .38+0.73e-1*(1+(0.525e-1/(0.3e-2)-2)^2/(0.525e-1/(0.3e-2))^2)

.5102677551

(22)

hi

hi := 3.6*kfluid*(dp*Gs/(`μ_air`*`εb`))^.365/dp

448.9928888

(23)

kinetic parameter

K := ln(19.837-13636/(T(z)+273.15))

ln(19.837-13636/(T(z)+273.15))

(24)

radial disperssion Coeff

Dr := Vs*dp/(9*(1+19.4*(dp/dt)^2))

0.3135309899e-3*Vs

(25)

thermal conductivity calculations

Kbs := kfluid*(`εb`+(1-`εb`)/(2/3*(1+kfluid/ksolid)))

0.5937050269e-1

(26)

``

Kbd := `εb`*Cpair*Gs*dp/(9*(1+19.4*(dp/dt)^2))

0.2342976727e-3

(27)

KB := Kbs+Kbd

0.5960480036e-1

(28)

Heat of reaction (deltaH), Find heat of formation for each reactant and product--> dHrxn = heat of formation(product)-heat of formation(reactant)

Error, missing operator or `;`

 

heat of formation = A + B*T + C*T^2

Hfacrolein := -7.076*10+(-5.59*10^(-2))*(273.15+350)+3.86*10^(-5)*(273.15+350)^2

-90.60509039

(29)

Hfwater := -238.41-0.122e-1*(350+273.15)+2.76*10^(-6)*(273.15+350)^2

-244.9406781

(30)

Hfpropylene := 3.62*10+(-6.49*10^(-2))*(273.15+350)+3.049*10^(-5)*(273.15+350)^2

7.59731748

(31)

`ΔHrxn` := Hfacrolein-Hfwater-Hfpropylene

146.7382702

(32)

Solving nodified Ergun Equation

f := (1-`εb`)*(1.75+(150*(1-`εb`))*`μ_air`/(dp*Gs))/`εb`

2.191735176

(33)

``

Determining U of the wall and Ueffective

``

``

Uwall := 1/(1/hi+ln(dto/dt)/(2*ksolid)+1/ho)

22.61972270

(34)

``

Ueff := 1/(1/Uwall+(1/2)*dt/(4*KB))

6.473624190

(35)

ode1 := Gs*(diff(Cprop(z), z))/rho = -K*Cprop(z)

0.4117046714e-2*(T(z)+273.15)*(diff(Cprop(z), z))/P(z) = -ln(19.837-13636/(T(z)+273.15))*Cprop(z)

(36)

ode2 := Gs*Cpair*(diff(T(z), z)) = -K*Cprop(z)*`ΔHrxn`-4*Ueff*(T(z)-350)/dt

1.4645*(diff(T(z), z)) = -146.7382702*ln(19.837-13636/(T(z)+273.15))*Cprop(z)-492.4968000*T(z)+172373.8800

(37)

ode3 := diff(P(z), z) = -f*Gs^2/(rho*dp)

diff(P(z), z) = -4.361346783*(T(z)+273.15)/P(z)

(38)

Ics1 := Cprop(0) = 0.3e-1

Cprop(0) = 0.3e-1

(39)

Ics2 := T(0) = 350

T(0) = 350

(40)

Ics3 := P(0) = 2.5

P(0) = 2.5

(41)

NULL

dsolve({Ics1, Ics2, Ics3, ode1, ode2, ode3}, {Cprop(z), P(z), T(z)})

Error, (in f) unable to store '-HFloat(0.020918994979034728)-HFloat(0.09184018333019917)*I' when datatype=float[8]

 

 

NULL

NULL

``

hi all

 

i have tow equations and i want plot q vs.A , that A change beetwin 0..phi. vn and a are definite and positive.

plz plz help me

 

(diff(f,y))^2-(f[y]^4)/4+q^2/(f[y]^2)+(1+2*q^2)*(f[y]^2)/2=piecewise(y<-a,1/4+2*q^2,-a<y<a,vn*(f[y]^2)-vn+1/4+2*q^2,y>a,1/4+2*q^2);

phiB[infinity]=q*int(f,y=-infinity..infinity):
A:=phiB[infinity]:

Hi

I have this attched procedure "Alt_lisse1".

Why the result e1 does'nt take into account the known variables

Thanks

ALt_lisse1.mw

Hi

is it possible to remove this warning and not to show it when executing a procedure ? 

Warning, `AzIII` is implicitly declared local to procedure `lisse`

Thanks

I have a question about using embedded component in Maple 2015.

I want to import the MapleSim model to Maple and use the "Combo Box" (which is one of the embedded component in Maple) to load the parameters that were used in the MapleSim model.

 

Embedded Component

 

The picture above shows the Maple worksheet that I made simply.

As you can see, I imported the MapleSim file (which contains three parameters, m, c, and k) to the Maple by using "MapleSim Component (Embedded Component)".

After that, I wanted to load these three parameters to the "Combo Box" automatically like the picture above, but I couldn't do that.

(In the case above, I manually type m, c, and k in the "Combo Box")

 

Isn't there any way to load the parameters to "Combo Box" automatically?

I attached the MapleSim file that I made. It also contains the Maple worksheet.  

[TS]EmbeddedComponents.msim 

 

I appreciate any idea you may have.

with(GroupTheory):
G := Group({[[1, 2]], [[3, 4]]});
Orbit(G);
G := Group({[[1, 2], [3, 4]]});
Orbit(G);

 

GroupTheory:-PermutationGroup(

{module () local cycles, p, d, work; option object; end module,

module () local cycles, p, d, work; option object; end module},

degree = 4)
Error, invalid input: GroupTheory:-Orbit expects its 1st argument, point, to be of type posint, but received module () local labels, minSupp, maxSupp, suppSize, AtkinsonsAlgorithm, IsSimpleGroupOrder, doDerivedSeries, doLowerCentralSeries, Intersection2, RightCosetRepresentatives, LeftCosetRepresentatives, PRA, `Giant?`, `Even?`, doStab1, doStab, CycleIndexMonomial; export generator_list, n, supergroup, Sylows, pCores, ModulePrint, ModuleDeconstruct, Generators, Orbit, Orbits, IsTransitive, Transitivity, IsPrimitive, GroupOrder, Elements, IsAbelian, IsElementary, IsSimple, ConjugacyClass, ConjugacyClasses, CayleyTable, Centre, DerivedSubgrou...
GroupTheory:-PermutationGroup(

{module () local cycles, p, d, work; option object; end module},

degree = 4)
Error, invalid input: GroupTheory:-Orbit expects its 1st argument, point, to be of type posint, but received module () local labels, minSupp, maxSupp, suppSize, AtkinsonsAlgorithm, IsSimpleGroupOrder, doDerivedSeries, doLowerCentralSeries, Intersection2, RightCosetRepresentatives, LeftCosetRepresentatives, PRA, `Giant?`, `Even?`, doStab1, doStab, CycleIndexMonomial; export generator_list, n, supergroup, Sylows, pCores, ModulePrint, ModuleDeconstruct, Generators, Orbit, Orbits, IsTransitive, Transitivity, IsPrimitive, GroupOrder, Elements, IsAbelian, IsElementary, IsSimple, ConjugacyClass, ConjugacyClasses, CayleyTable, Centre, DerivedSubgrou...

how to calculate Ferrari resolvent of x^4-c1*x^3+c2*x^2-c3*x+c4

hi.how i can dsolve couple linear equations with power series solutions or taylor series expantion?

file attached below.

thanks

 

TAYLOR.mw

1 2 3 4 5 6 7 Last Page 3 of 1265