Maple 17 Questions and Posts

These are Posts and Questions associated with the product, Maple 17

Dear all,

I developed a program to solve f(x, y) = 0 and g(x, y) = 0, I obtained as results (x=2.726, y=2.126) . running the same program another time it gives (x=2.762, y=1.992). how to explain this?

> fsolve({f(x, y) = 0, g(x, y) = 0}, {x = 0 .. infinity, y= 0 .. infinity});

Thanks in advance.

So have to calculate frequency response of RLC-circuit
Long story short - I can't get how to solve something like this
solve(abs(I+x) = 1)
Obv - the answer is 0 or -I, but i get 1-I, -1-I instead

It's been a couple of years since I last used the Excel add-in for Maple 17. It was working before, but now it isn't.

When I click any of the Maple icons in the Excel Add-Ins -> Custom Toolbars tab, I get a Visual Basic run-time error 453, "Can't find DLL entry point vbLoadMapleCursor in WMIMPLEX". The only choice in the dialogue is End, after which none of the Maple icons respond.

When I quit Excel after trying this, I get the same error dialogue except the message refers to vbMapleQuit. I select End, and then Excel exits. (If I don't touch the Maple icons during an Excel session, then of course no VB error dialogue referring to Maple pops up.)

An answer to a ten year old MaplePrimes question on the Excel add-in suggests editing the Excel 2010 32-bit path registry entry. Indeed my HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Windows\CurrentVersion\App Paths\Excel.exe Path key has no reference to Maple. However that info is so old, for Maple 9.5, that I hesitate to use it. I tried though. As suggested in that question, updated for Maple 17, I tried adding ";C:\Program Files (x86)\Maple 17\Excel\" to the end of the existing Path key entry, but I still get the same VB error. I also tried adding ";C:\Program Files (x86)\Maple 17\bin.win\", and both: "C:\Program Files (x86)\Maple 17\bin.win\;C:\Program Files (x86)\Maple 17\Excel\". I tried disabling and reloading the add-in from the Excel Add-Ins menu. Interestingly, even when I disable the add-in, the toolbar icons still show.

Anyway I'm at a loss how to proceed. Any help is appreciated. --Graham

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

``

Hello everyone, how can this error be corrected

pde.mw

Thanks.

I'm not really an expert in maple. I'm stuck with a problem that could be related to the precision of the computational engine...

I defined a function which is relatively complicated as follows:

m := proc (x, l, T) options operator, arrow; (exp(l*(x-T)/(1-exp(-l*T)))*exp(l*T/(1-exp(-l*T)))*(1-exp(-l*T))+(2-exp(-l*T)-exp(l*T))*exp(l*(x-T)/(1-exp(-l*T))))/(exp(l*T)-1)-(l*(x-T)*exp(l*(x-T)/(1-exp(-l*T)))-2*exp(-l*T)+4-2*exp(l*T))/(exp(l*T)-1) end proc

whenever i simplify the equation or factor some of the terms, I get completely different results.

As an example, you can define this function as follows

m := proc (x, l, T) options operator, arrow; (exp(x*l/(1-exp(-l*T)))*(1-exp(-l*T))+(2-exp(-l*T)-exp(l*T))*exp(l*(x-T)/(1-exp(-l*T))))/(exp(l*T)-1)-(l*(x-T)*exp(l*(x-T)/(1-exp(-l*T)))-2*exp(-l*T)+4-2*exp(l*T))/(exp(l*T)-1) end proc

 

By drawing the function, you can see how things get messy.

plot(m(450, (1/250)*x, 250), x = 0 .. 100)

 

I know the problem occurs as a results of computation round-off. How could I ameliorate that? I have already increase the precision in the option menu, but that didn't help. By the way, how can I be sure that the anwers maple gives me is correct, except using my intuition?

 

Hello everyone!

i want to make a triangle which height's is a. Sorry, my English is not very good!

And I writed it by C language, but It is not right in maple. and I don't know. Can you help me? thank you very much! triangle.mw

with sin(x) and cos(x). how to draw  this graph:

 

question.mw

``


thanks you!!!

Download question.mw

why the the software can't plot the function like x^(4/3)*sin(1/x) or x^(1/3)

it could only plot where x>0,but the value is does exist where x<0.

Thanks in advance for your help.

Hi Everyone,

I want to draw a 3D cylindrical object with follwoing constraint:

(1) For z=0 to 5, it should be a cylinder with radius 2 units

(2) For z = 10 to 20, it should be a cylinder with radius 8 units

(3) Connecting piece between z=5 to 10 is a type of cylinder such that its radius increses from 2 to 8 as z moves from 5 to 10.

 

I know how to draw (1) and (2) individually, but can't make a connecting object.

Can anyone help me drawing that object?

Thanks in advance.

hello

I have recently install Maple17 on my computer (Windows10) and I need to use some Greece alphabet such as ß but I look everywhere in maple's icon and I just could find capital Greece alphabet.

does anybody know how can I find those?

 

 

The plot command outputs a graph of symbols where the line is made up of "H"s, vertical lines on the axes are "+"s, horizontal lines on the axes are "-"s, and intercepts are "*"s.

 

How do I fix this problem so that the graph displays a graph normally?

I have several *.mw files which I have built seperately. Now, for readability reasons, I would like to create a Master document, which will ues the first files, and then generate a PDF document.
I tried the read command, but it does not recognize the files.

Any help is appreciated.

 

Hello every one. I writed:

M:=-2,-1,0,1,2:
f:=x->2^x
with(LinearAlgebra);
for i from 2 by 1 to nops(M) do
for j from 2 by 1 to i
Q(i,j)->((x-M[i-j])*Q(i,j-1)-(x-M[i])*Q(i-1,j-1))/(M[i]-M[j-1]):
A:=Matrix(Q,nops(M),nops(M)+1):

end do;
end do;
print(A)
But I don't know wrong in my code. That is Neville algorithm

can you help me? sorry, my English is not very good

can someone tell me how can I derive and solve d(g1(x))=0 a complicated function as

I tried with the instructions simplify and combine but without results.

Thanks in advance.

 

 

First 23 24 25 26 27 28 29 Last Page 25 of 61