Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

The following is the log file after running a code file,

Extrema_Network := table([])

Extrema_Network[Han] := [20.0, 385.61]

Extrema_Network[Liv] := [20.0, 385.61]

Extrema_Network[Vir] := [20.0, 385.61]

values := [[2.8274333874308139146163, 2 Pi]]

theta_step := 0

phi_step := 0

counter_theta := 0

counter_phi := 0


distance eff distance_eff
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
im in has not
Im in has par
Im in has par
the time offset for Han was .71848e-2*cos(lat_begin)*cos(long_begin)+.12746e-1*cos(lat_begin)*sin(long_begin)-.15394e-1*sin(lat_begin)
the time offset for Liv was .24690e-3*cos(lat_begin)*cos(long_begin)+.18302e-1*cos(lat_begin)*sin(long_begin)-.10809e-1*sin(lat_begin)
the time offset for Vir was -.15117e-1*cos(lat_begin)*cos(long_begin)-.28029e-2*cos(lat_begin)*sin(long_begin)-.14657e-1*sin(lat_begin)
the angles with the coordinate axis are [long_begin-3.1416, 1.5708+lat_begin]
Using logChirp

common factor in Signal() was .43329e-21

Using logChirp

common factor in Signal() was .43329e-21

Using logChirp

common factor in Signal() was .43329e-21

snr_network was (187.00-.48032e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*sin(long_begin)-378.50*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*sin(long_begin)-247.94*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)^3-56.997*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)-31.130*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)^3+32.512*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*sin(long_begin)+196.54*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)+.65598e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)-.13119e-2*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)^3+64.585*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*cos(long_begin)^2*sin(long_begin)-61.911*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*cos(long_begin)^2*sin(long_begin)+110.68*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^3*sin(long_begin)^2*cos(long_begin)-375.04*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)*sin(long_begin)^2*cos(long_begin)+.85041e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*cos(long_begin)^2*sin(long_begin)+.10495e-3*(1.+sin(lat_begin))^(1/2)*(1.-1.*sin(lat_begin))^(1/2)*sin(lat_begin)^2*sin(long_begin)^2*cos(long_begin)-.49247e-3*sin(lat_begin)^3+260.07*cos(long_begin)^4+.16559e-3*cos(long_begin)^4*sin(lat_begin)+.16559e-3*sin(lat_begin)^3*cos(long_begin)^4+.90222e-3*sin(lat_begin)^3*cos(long_begin)^2-.12334e-2*sin(lat_begin)*cos(long_begin)^2-268.44*cos(long_begin)^2+102.34*sin(long_begin)*cos(long_begin)-590.30*sin(lat_begin)^2*cos(long_begin)^2-312.01*sin(lat_begin)^4*cos(long_begin)^2-102.19*sin(lat_begin)^4*sin(long_begin)^2+260.07*sin(lat_begin)^4*cos(long_begin)^4+650.59*sin(lat_begin)^2*cos(long_begin)^4+114.40*cos(long_begin)^3*sin(long_begin)+33.292*sin(long_begin)^2*cos(long_begin)^2-.11326e-2*cos(long_begin)^3*sin(lat_begin)*sin(long_begin)-.85371e-4*sin(long_begin)^2*cos(long_begin)^2*sin(lat_begin)-.11326e-2*sin(lat_begin)^3*cos(long_begin)^3*sin(long_begin)+.49010e-3*sin(lat_begin)^3*sin(long_begin)*cos(long_begin)-.85371e-4*sin(lat_begin)^3*sin(long_begin)^2*cos(long_begin)^2-394.84*sin(lat_begin)^2*sin(long_begin)*cos(long_begin)+287.80*sin(lat_begin)^4*sin(long_begin)*cos(long_begin)+114.40*sin(lat_begin)^4*cos(long_begin)^3*sin(long_begin)-219.38*sin(lat_begin)^2*cos(long_begin)^3*sin(long_begin)+33.292*sin(lat_begin)^4*sin(long_begin)^2*cos(long_begin)^2+1085.5*sin(lat_begin)^2*sin(long_begin)^2*cos(long_begin)^2+.64246e-3*sin(lat_begin)*sin(long_begin)*cos(long_begin)+157.13*sin(long_begin)^2-54.941*sin(lat_begin)^2*sin(long_begin)^2+78.538*sin(lat_begin)^4+27.150*sin(lat_begin)^2+.57530e-3*sin(lat_begin))^(1/2)
low_fisher for detector Han was


Error, (in fprintf) number expected for floating point format
13

32

31

17

hou := 0

mini := 0

seci := 0

memory used=26.3MB, alloc=32.3MB, time=0.47




Not quite sure why I am getting this error. Also this code was written in maple 14 and I am using maple 17. Could that also explain why I am getting all these errors? Is it possible to downgrade maple 17 to maple 14? Sorry if my question above on fprintf is a bit too cryptic as it is research work. I will upload the code file in a reply.

DEAR SIR,

PLEASE HELP ME WITH THAT QUESTION

 

DEAR SIR

ANYONE CAN HELP TO COMPUTE TIME IN DSOLVE COMMAND?

restart

with(plots)

Nb := 0.1e-4; Nt := 0.1e-4; Sc := 3.0; Sb := 15; Pe := 1; Bi := .5; Pr := 6.8; c[4] := 0; c[6] := .3; c[8] := .4; k[1] := 0; k[2] := 1; k[3] := 0; Un := .1; M := .5

Eq1 := (101-100*lambda)*(1+c[2]*phi(eta))*(diff(f(eta), `$`(eta, 3)))+(diff(f(eta), `$`(eta, 2)))*(f(eta)+g(eta)+c[2]*(diff(phi(eta), eta)))-(diff(f(eta), eta))^2-M*(diff(f(eta), eta)-k[2])+k[2]+Un*(k[2]-(1/2)*eta*(diff(f(eta), `$`(eta, 2)))-(diff(f(eta), eta))) = 0

(101-100*lambda)*(1+c[2]*phi(eta))*(diff(diff(diff(f(eta), eta), eta), eta))+(diff(diff(f(eta), eta), eta))*(f(eta)+g(eta)+c[2]*(diff(phi(eta), eta)))-(diff(f(eta), eta))^2-.6*(diff(f(eta), eta))+1.6-0.5000000000e-1*eta*(diff(diff(f(eta), eta), eta)) = 0

(1)

Eq2 := (101-100*lambda)*(1+c[2]*phi(eta))*(diff(g(eta), `$`(eta, 3)))+(diff(g(eta), `$`(eta, 2)))*(f(eta)+g(eta)+c[2]*(diff(phi(eta), eta)))-(diff(g(eta), eta))^2-M*(diff(g(eta), eta)-k[2])+k[2]+Un*(k[2]-(1/2)*eta*(diff(g(eta), `$`(eta, 2)))-(diff(g(eta), eta))) = 0

(101-100*lambda)*(1+c[2]*phi(eta))*(diff(diff(diff(g(eta), eta), eta), eta))+(diff(diff(g(eta), eta), eta))*(f(eta)+g(eta)+c[2]*(diff(phi(eta), eta)))-(diff(g(eta), eta))^2-.6*(diff(g(eta), eta))+1.6-0.5000000000e-1*eta*(diff(diff(g(eta), eta), eta)) = 0

(2)

Eq3 := (1+c[4]*phi(eta))*(diff(theta(eta), `$`(eta, 2)))+Pr*(diff(theta(eta), eta))*(f(eta)+g(eta))+Nb*Pr*(diff(theta(eta), eta))*(diff(phi(eta), eta))*(1+c[6]*(2*phi(eta)+1))+Nt*Pr*(diff(theta(eta), eta))^2+c[4]*(diff(theta(eta), eta))*(diff(phi(eta), eta))-(1/2)*Pr*eta*Un*(diff(theta(eta), eta)) = 0

diff(diff(theta(eta), eta), eta)+6.8*(diff(theta(eta), eta))*(f(eta)+g(eta))+0.68e-4*(diff(theta(eta), eta))*(diff(phi(eta), eta))*(1.3+.6*phi(eta))+0.68e-4*(diff(theta(eta), eta))^2-.3400000000*eta*(diff(theta(eta), eta)) = 0

(3)

Eq4 := (1+c[6]*phi(eta))*(diff(phi(eta), `$`(eta, 2)))+Sc*(f(eta)+g(eta))*(diff(phi(eta), eta))+c[6]*(diff(phi(eta), eta))^2+Nt*(diff(theta(eta), `$`(eta, 2)))/Nb-(1/2)*Sc*eta*Un*(diff(phi(eta), eta)) = 0

(1+.3*phi(eta))*(diff(diff(phi(eta), eta), eta))+3.0*(f(eta)+g(eta))*(diff(phi(eta), eta))+.3*(diff(phi(eta), eta))^2+1.000000000*(diff(diff(theta(eta), eta), eta))-.1500000000*eta*(diff(phi(eta), eta)) = 0

(4)

Eq5 := (1+c[8]*phi(eta))*(diff(chi(eta), `$`(eta, 2)))+c[8]*(diff(phi(eta), eta))*(diff(chi(eta), eta))-Pe*(chi(eta)*(diff(phi(eta), `$`(eta, 2)))+(diff(phi(eta), eta))*(diff(chi(eta), eta)))+Sb*(diff(chi(eta), eta))*(f(eta)+g(eta))-(1/2)*Sb*eta*Un*(diff(chi(eta), eta)) = 0

(1+.4*phi(eta))*(diff(diff(chi(eta), eta), eta))-.6*(diff(phi(eta), eta))*(diff(chi(eta), eta))-chi(eta)*(diff(diff(phi(eta), eta), eta))+15*(diff(chi(eta), eta))*(f(eta)+g(eta))-.7500000000*eta*(diff(chi(eta), eta)) = 0

(5)

Vc[2] := [.2, .3, .4]

etainf := 1.85

bcs := (D(f))(0) = k[1], (D(g))(0) = k[3], f(0) = 0, g(0) = 0, (D(theta))(0) = -Bi*(1-theta(0))/(1+c[4]*phi(0)), Nb*(D(phi))(0)*(1+c[6]*(2*phi(0)+1))+Nt*(D(theta))(0) = 0, chi(0) = 1, (D(f))(etainf) = k[2], (D(g))(etainf) = k[2], theta(etainf) = 0, phi(etainf) = 0, chi(etainf) = 0

(D(f))(0) = 0, (D(g))(0) = 0, f(0) = 0, g(0) = 0, (D(theta))(0) = -.5+.5*theta(0), 0.1e-4*(D(phi))(0)*(1.3+.6*phi(0))+0.1e-4*(D(theta))(0) = 0, chi(0) = 1, (D(f))(1.85) = 1, (D(g))(1.85) = 1, theta(1.85) = 0, phi(1.85) = 0, chi(1.85) = 0

(6)

dsys := {Eq1, Eq2, Eq3, Eq4, Eq5, bcs}

for i to 3 do c[2] := Vc[2][i]; dsol[i] := dsolve(dsys, numeric, continuation = lambda); print(c[2]); print(dsol[i](0)) end do

.2

 

[eta = 0., chi(eta) = HFloat(1.0), diff(chi(eta), eta) = HFloat(3.888290578689045), f(eta) = HFloat(0.0), diff(f(eta), eta) = HFloat(0.0), diff(diff(f(eta), eta), eta) = HFloat(2.244199282192492), g(eta) = HFloat(0.0), diff(g(eta), eta) = HFloat(0.0), diff(diff(g(eta), eta), eta) = HFloat(2.244199282192492), phi(eta) = HFloat(-2.044191234673432), diff(phi(eta), eta) = HFloat(5.227515304629519), theta(eta) = HFloat(0.2317093657771352), diff(theta(eta), eta) = HFloat(-0.38414531711143246)]

 

.3

 

[eta = 0., chi(eta) = HFloat(1.0), diff(chi(eta), eta) = HFloat(4.148187853914835), f(eta) = HFloat(0.0), diff(f(eta), eta) = HFloat(0.0), diff(diff(f(eta), eta), eta) = HFloat(3.1086884419918364), g(eta) = HFloat(0.0), diff(g(eta), eta) = HFloat(0.0), diff(diff(g(eta), eta), eta) = HFloat(3.1086884419918364), phi(eta) = HFloat(-2.049332060722701), diff(phi(eta), eta) = HFloat(5.527786294980874), theta(eta) = HFloat(0.2216792480031605), diff(theta(eta), eta) = HFloat(-0.38916037599841974)]

 

.4

 

[eta = 0., chi(eta) = HFloat(0.9999999999999998), diff(chi(eta), eta) = HFloat(4.580796631072469), f(eta) = HFloat(0.0), diff(f(eta), eta) = HFloat(0.0), diff(diff(f(eta), eta), eta) = HFloat(5.687607599246298), g(eta) = HFloat(0.0), diff(g(eta), eta) = HFloat(0.0), diff(diff(g(eta), eta), eta) = HFloat(5.687607599246298), phi(eta) = HFloat(-2.0568809171520708), diff(phi(eta), eta) = HFloat(6.0203396482123575), theta(eta) = HFloat(0.20686299926628238), diff(theta(eta), eta) = HFloat(-0.39656850036685887)]

(7)

NULL

 

 

Download 3DAKc2w_-_Copy.mw

I have a gensym routine (gensym = generate symbol) which appends Chinese characters to some base symbols to create 'new' symbols.  This works well on Windows, Mac and tty maple.  But I get all square boxes in the linux GUI.

I know the fonts exist as, on the same machine, it works in tty.  But I have no idea how to tell the GUI to go and use those fonts.  Nor why it isn't.

Amusingly, if I try to paste something from the GUI as MapleMath here on Primes, in the 'paste' box, it shows just fine, but then it says it is invalid Maple.

For example, from the GUI I see

but the same thing in TTY is

(x三0不 + x三1下)^(-M)

and I just did cut-and-paste of the same thing.  By this I mean that I took that output with boxes above, put it into a terminal window, re-copied that, pasted it in here, and voila!

So clearly the issue is with the GUI, and only on linux (Ubuntu 16.04).

I am curious about the numerical method(s) used by Maple to calculate fsolve. I've looked at the documentation (https://www.maplesoft.com/support/help/maple/view.aspx?path=fsolve%2fdetails), but the method(s) used are not stated on there. Does anyone know which method Maple uses in fsolve? Additionally, does Maple use multiple methods and if so, how does it determine which one to use? Thanks!

My task is Second Order Implicit Differentiation (if my english is correct)

My paper's photoI tried to accomplish it usig maple:

ScreenShoot of my maplebut my  mixed partials  are different but they should be equal if i used the correct expression in mapple.. Why is that or simply how do I lolve my problem using mapple?

Help pleease!!

We've added a collection of thermal engineering applications to the Application Center. You could think of it as an e-book.

This collection has a few features that I think are pretty neat

  • The applications are collected together in a Workbook; a single file gives you access to 30 applications
  • You can navigate the contents using the Navigator or a hyperlinked table of contents
  • You can change working fluids and operating conditions, while still using accurate thermophysical data

If you don't have Maple 2016, you can view and navigate the applications (and interactive with a few) using the free Player.

The collection includes these applications.

  • Psychrometric Modeling
    • Swamp Cooler
    • Adiabatic Mixing of Air
    • Human Comfort Zone
    • Dew Point and Wet Bulb Temperature
    • Interactive Psychrometric Chart
  • Thermodynamic Cycles
    • Ideal Brayton Cycle
    • Optimize a Rankine Cycle
    • Efficiency of a Rankine Cycle
    • Turbine Analysis
    • Organic Rankine Cycle
    • Isothermal Compression of Methane
    • Adiabatic Compression of Methane
  • Refrigeration
    • COP of a Refrigeration Cycle
    • Flow Through an Expansion Valve
    • Food Refrigeration
    • Rate of Refrigerant Boiling
    • Refrigeration Cycle Analysis 1
    • Refrigeration Cycle Analysis 2
  • Miscellaneous
    • Measurement Error in a Manometer
    • Particle Falling Through Air
    • Saturation Temperature of Fluids
    • Water Fountain
    • Water in Piston
  • Heat Transfer
    • Dittus-Boelter Correlation
    • Double Pipe Heat Exchanger
    • Energy Needed to Vaporize Ethanol
    • Heat Transfer Coefficient Across a Flat Plate
  • Vapor-Liquid Equilibria
    • Water-Ethanol

I have a few ideas for more themed Maple application collections. Data analysis, anyone?

I have worksheet that generates very large matrices and it seems to save all results, so it requires large amount of RAM to load and storage to save it.  

There used to be a pop-up prompt not to save the results of the worksheet - to save space, but that no longer occurs.  Can someone help me to switch on that prompt not to save results or otherwise prevent saving results, when I open the spreadsheet.  It maybe that I agreed not to be prompted again, but now I am regretting it!!

MRB 

equidistant_curve_MP.mw  Equidistant curves to the curves on the surface. (Without any sense, but real.)







hello

I have a equation of 4th degree with two variable (x,y) that i want to solve my equation according to one of variable (y).

myeq:= y^4 +(-(5/96)*x^2 -(5/6)*x -2) *y^3 +(-(55/288)*x^2 -(7/12)*x +1) *y^2 +((65/288)*x^2 +(2/3)*x) *y +(5/288)*x^2 -(1/4)*x;

i use "solve(myeq,y)" but can't solve. how can i solve this equation?

I watched a webinar training:   Maple Training for Engineers, Researchers, and Scientists.

 

The file demonstrated topics and used collapsible regions with a colored band across the document and a + / - sign on the left to open/close the section or region.  Different than the document gray down arrow.

 

I have searched all the documentation and can not find how to do this.   I have done this in other programs like MathCAD to close and lock a region from viewing or editing of specific formulas.

 

Any help always appreciated.

 

Regards,

Bill

I'm currently working on building a Grid Layout for a project, and I'm having trouble coding in the RunWindow and GetFile elements into buttons under the grid layout. I've gone through the overviews and examples for them, but had no luck. I'm using Maple 2016.1 for OS X.

Additionally, the structure of the code is slightly different as to how many of the example worksheets structure their Grid Layout code, since the code originated from a Maplet Builder file. I.e. in the example worksheets they would follow as:

maplet := Maplet('onstartup' = 'Action1', 'reference' = 'Maplet1',
         BoxLayout('background' = "#D6D3CE", 'border' = 'false', 'halign' = 'center', 'inset' = '5', 'reference' = 'BoxLayout1', 'valign' = 'center', 'vertical' = 'false', 'visible' = 'true',
                       BoxColumn( BoxCell('hscroll' = 'never', 'value' = 'Button1', 'vscroll' = 'never'),
         GridLayout('background' = "#D6D3CE", 'border' = 'false','halign'='center','inset'='5', 'reference' = 'GridLayout1', 'valign' = 'center', 'visible' = 'true',
                   GridRow('valign' = 'top', GridCell('height' = '1', 'hscroll' = 'never', 'value' = 'BoxLayout1', 'vscroll' = 'never', 'width' = '1' ))),
         Window('layout'= 'GridLayout1', 'reference' = 'W1', 'resizable' = 'true', 'title' = "Maplet"),
          Action('reference' = 'Action1', RunWindow('window'= 'W1'))

However the structure for the code I am working with has action at the very start of the code, follwed by the the code for the buttons then layouts/window.  E.g. (the code has been shortened)

with (Maplets[Elements]):
maplet :=
Maplet('onstartup'='Action1','reference'='Maplet1',
Action('reference'='clickButton1'),
Action('reference'='clickButton9',
Evaluate('function'='plot3d(x^2*cos(y),x = -1 .. 1,y = -2*Pi .. 2*Pi)','target'='Plotter1','waitforresult'='true')),
Action('reference'='clickButton11'),
Action('reference'='clickButton12'),
Action('reference'='clickButton10'),
Button('background'="#D6D3CE",'caption'="Insert Molecular Geometry",'enabled'='true','foreground'="#000000",'onclick'='clickButton1','reference'='Button1','visible'='true'),

....

BoxLayout('background'="#D6D3CE",'border'='false','halign'='center','inset'='5','reference'='BoxLayout1','valign'='center','vertical'='false','visible'='true',
BoxColumn(
BoxCell('hscroll'='never','value'='Button1','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout2','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout3','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout9','vscroll'='never'),
BoxCell('hscroll'='never','value'='BoxLayout14','vscroll'='never')),
BoxColumn(
BoxCell('hscroll'='never','value'='Label3','vscroll'='never'),
BoxCell('hscroll'='never','value'='Plotter1','vscroll'='never'),
BoxCell('hscroll'='never','value'='Slider1','vscroll'='never'))),
GridLayout('background'="#D6D3CE",'border'='false','halign'='center','inset'='5','reference'='GridLayout1','valign'='center','visible'='true',
GridRow('valign'='top',
GridCell('height'='1','hscroll'='never','value'='BoxLayout1','vscroll'='never','width'='1'))),
Window('layout'='GridLayout1','reference'='Window1','resizable'='true','title'="Maplet"),
Action('reference'='Action1',
RunWindow('window'='Window1'))):

Maplets[Display](maplet);

 

If anyone would be able to provide an example of code or some guidance I could follow that would be greatly appreciated! 

hello

i have some implicit plot that i want intersect between them.

for example:


implicitplot(x^2+y^2<1,x=-3..3,y=-3..3);
implicitplot((x-1)^2+y^2<1,x=-3..3,y=-3..3);

how can i do intersect between them!?

thanks for guidance

Hi there,

Recently, I started to learn maple. I came across a problem. 

I entered following codes:

with(plots);

ball := proc (x, y) fieldplot([x, y], x = -1 .. 1, y = -1 .. 1) end proc;

ball(1,1);      


However, it showed me Error, (in plots/animate) bad range arguments -10. = -10 .. 10, -10. = -10 .. 10.

could someone tell me how to put fieldplot into procedure? Very much appreciated.

I am trying to use y* to label a point on the axis of a graph made with DEplot, and am currently unable to.


with(DEtools);
NLC := diff(y(t), t) = k*(Am-y(t));
Am := 20; k := .1;
ivs := [y(0) = 10, y(0) = 30, y(0) = 50];
DEplot(NLC, y(t), t = 0 .. 20, ivs, tickmarks = [default, [20 = y^`*`]], font = [default, default, 30]);

makes y`*` apear as the label, as does the code

tickmarks = [default, [20 = y^`&ast;`]]

wheras if i remove the `` marks I get an error

 

 

 

First 1087 1088 1089 1090 1091 1092 1093 Last Page 1089 of 2224