MaplePrimes Questions

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!!

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 

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

 

 

 

how to use Appell hypergeometric function in maple? Is there any command for that like hypergeom ???

I have the PDE u_{xx}+u_{yy} = 1 with BC: u|_{x^2+y^2=1} =0 ;

 

how to write down the command of the BC in solving this PDE?, btw can I make maple show me how to solve this PDE analytically?

 

Thanks in advance.

 

Here are the lines that I wrote so far:

pde := diff(u(x, y), x, x)+diff(u(x, y), y, y) = 1;

ans := pdsolve(pde)

 

how to add the BC correctly to pdsolve? I am not sure how to write the condition x^2+y^2=1 and that u will get a value on this boundary.

 

First 1112 1113 1114 1115 1116 1117 1118 Last Page 1114 of 2434