Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Maple's isprime is not a definitive primality test. The input has to pass a "strong pseudo-primality test" and "one Lucas test". This is well documented. I thought I remembered that there is also a way to get Maple to perform a true primality test, but I don't remember how and don't see anything about this in the Maple help system.

Is my memory faulty, or is there no definitive primality test in Maple?

Thanks in advance,

Doug

---------------------------------------------------------------------
Douglas B. Meade  <><
Math, USC, Columbia, SC 29208  E-mail: mailto:meade@math.sc.edu
Phone:  (803) 777-6183         URL:    http://www.math.sc.edu

Hello,

How to change the length underline in latex ?

For example in this environment :

$\underline{\mathcal{G}}$

 

Thank you,

Gérard.

So I am trying to apply the Optimization function to a rather complicated problem I am currently working on, and having some trouble getting maple to cooporate. This is the part of the code that is giving me the error, with error included. All of the variables imputtet when calling OptimizeSpring are constants. If there is any other details I should add, please say so and I will add them promptly :) Thanks alot!
Maple problem

I suspect that it might be related to this suggestion, but I am not sure how to apply it http://www.mapleprimes.com/posts/95113-Optimization-Workaround-For-Operator-Form

I'm trying to solve a beam vibration problem. I am relatively new to Maple and trying to solve PDE's. I feel like I'm pretty close, but I don't understand what's wrong here.

This is my maple file.

Thesis_Pde2_attempt.mw

as a note, I made some simplifying assumptions. I am just using part of this equation since it is less complex. And my boundary conditions are simply supported on both ends so no displacement and no 2nd derivative at the ends.

 

Any help would be greatly appreciated!!

Hello everyone, 

I have some problems with the "isolve" command on Maple. I am trying to solve for integer a very easy system of equations. When I type the commands

 

restart; 

n := 2;
isolve({sum(a[k], k = 1 .. n)-1 = 1}, d)


I get the expected {a[1] = 2-d, a[2] = d}. However, if I add conditions a[1],a[2] >= 0, that is the commands



restart;
n := 2;
isolve({ge(a[1], 0), ge(a[2], 0), sum(a[k], k = 1 .. n)-1 = 1}, d)


I get the warning "Warning, solutions may have been lost". What am I doing wrong? Is there a way to get Maple to give me the possible values?

 

Thank you in advance,

David

I'm trying to solve a system of equations thats expressed as a summation (the original has the summation symbol as opposed to 'sum'):

d_actual := solve(W_actual = sum(W_guess(def-asp_rad_inverse[i], E, asp_rad[i]), i = 1 .. n), def);

When n<5, I get an answer after a few seconds, but when n is higher, the program sits and 'evaluates' forever...I've waited up to 30 min. 

Background:

I have a plate with a number of hemispheres on the surface (# of impacted hemispheres given a force = n). Each has a unique radius and they're listed from largest to smallest in 'asp_rad[]'. I have an equation for the deflection of a single hemisphere as a function of Force and material properties that I have rearranged with respect to Force (W_guess(deflection, E, r)). There's an opposing plate that stays parallel to the original plate while pushing down on the asperities with Force W_actual. The total deflection of the opposing plate is the sum of: 1) the difference between the tallest radius and the impacted radius in question (asp_rad_inverse[]), and 2) the deflection of the impacted radius in question. 

I'm attempting to solve for the total deflection of the opposing plate via solving for the 'def' in the summation above, but when I run it, the program is not able to compute a solution.

Thanks in advance.

Hello,

I'm a very new beginner, and I have a question about getting the data of a curve from a JPG, without using manual measures an too getting more accurate datas.

I wish to fit the curve with a model. The JPG joined, from EU 2015, is about human thermoregulation. I still try to found a model with Mathematica but am uncertain seeing the result as illogic.

Did anyone write something to read curves from a JPG and get coordinates from these curves ? In the JPG, the curve is above, as an 'U'.

Thanks for your advice,

Milos.

Hello,

I have an error in my code. I don't know where it came from. Earlier today I loaded this and it worked fine and now an error comes up. Any help is greatly appreciated! Thank you in advance!

Kind regards,

Gambia Man

with(plots):

``

HyperionOrbit := proc (`&theta;IC`, `&omega;IC`) local a, Mh, Msat, G, e, beta, M, Eqns, Th, ICs, soln, `&omega;H`; a := 1.501*10^9; Mh := 5.5855*10^18; Msat := 5.6832*10^26; G := 6.67259/10^11; e := .232; beta := .89; M := Mh+Msat; Eqns := diff(theta(t), t) = omega(t), diff(omega(t), t) = -G*Msat*beta^2*(xH(t)*sin(theta(t))-yH(t)*cos(theta(t)))*(xH(t)*cos(theta(t))+yH(t)*sin(theta(t)))/(xH(t)^2+yH(t)^2)^2.5, diff(xH(t), t) = vxH(t), diff(vxH(t), t) = -G*M*xH(t)/(xH(t)^2+yH(t)^2)^(3/2), diff(yH(t), t) = vyH(t), diff(vyH(t), t) = -G*M*yH(t)/(xH(t)^2+yH(t)^2)^(3/2); Th := sqrt(4*Pi^2*a^3/(G*(Mh+Msat))); `&omega;H` := 2*Pi/Th; ICs := xH(0) = a*(1+e), yH(0) = 0, vxH(0) = 0, vyH(0) = sqrt(G*M*(1-e)/(a*(1+e))), theta(0) = `&theta;IC`, omega(0) = `&omega;IC`; soln := dsolve({Eqns, ICs}, numeric); odeplot(soln, [theta(t), omega(t)/`&omega;H`], 0 .. 5*Th, numpoints = 2000, labels = ["&theta;(t)","&omega;(t)/&omega;H"], axes = boxed) end proc

HyperionOrbit(.5, 1.8*`&omega;H`)

Warning, The use of global variables in numerical ODE problems is deprecated, and will be removed in a future release. Use the 'parameters' argument instead (see ?dsolve,numeric,parameters)

 

Error, (in dsolve/numeric/SC/IVPpoints) unable to store '11.3097335529232/Th' when datatype=float[8]

 

``

``

``

 

Download Poincare_section.mw

Hi everyone,

 

Consider this maple 18 doc: Euler18.mw

 

The code is regular code for Julia sets of the exponential.

 

To see how the Julia set behaves at infinity, I apply the transform mu(z)=1/z.

 

The plot3d command correctly plots the Julia set at an appropriate neighborhood of infinity, but:

1) Axes are not transformed

2) Saving as .eps produces an empty plot, sans the axes (plot is saved correctly, when not applying mu(z))

 

Is there any trick to force the axes to also show correctly with the transformed ranges?

Seems that this misalignment is bothering the .eps renderer, which probably plots the graph in twilight zone, given the false ranges of the untransformed axes.

 

Any ideas on how to force the saveas .eps to work in this case?

 

Many thanks,

Yiannis

Hi

If I have this equation

and I need to use different Z, it's like

how can I solve it

Many Thanks

Amal

Dear Maple users

Physical experiment: I dropped a ball with low mass from a height of approximately 7 meters and wanted to test if the air resistance was proportional to the square of the velocity. I filmed the fall and used the program Logger Pro to collect data: a number of datapoints (time,height) was collected. I copy/pasted the datapoints into MS Excel, from where I could import data into Maple via Tools > Assistants > Import Data ... Then I wanted to make a fit with the theoretical solution, given by a function having just one parameter: the Drag coefficient. Unfortunately I received an error "complex values encountered" (see below). I can solve the problem manual by making a number of guesses for the drag coefficient, until the theoretical curve approximates the data points well. I wanted to make Maple do the fitting job, though. I will appreciate if someone could give an idea how to fit the data properly.

NB! Mass m and g is defined above in the Maple document. The Statistics and plots package is called too.

I would like to plot the following singular double integral, but I cannot due to singularities...

 

where x>0, t=0.2 and m=0.2.

I defined f(y) function as f:=y->exp(-(y-4.68)^2/0.4):

I attached my file:
1st_try.mw

Thank you !

There is an error in my implementation as follow:

"Error, (in unknown) incorrect syntax in parse: `*` unexpected (near 1st character of parsed string)"

What I have to do to remove this error?

Hello,

I've got this error in my code and I don't know why as I didn't get it when I used a different xn function. Any help would be greatly appreciated! Thank you in advance!

Kind regards,

Gambia Man

restart

with(plots):

boxcount := proc (data, N) local n, xmax, xmin, xrange, ymax, ymin, yrange, dx, dy, i, j, ix, iy, sum, res; n := (1/2)*ArrayNumElems(data); xmax := max(seq(data[i, 1], i = 1 .. n)); xmin := min(seq(data[i, 1], i = 1 .. n)); ymax := max(seq(data[i, 2], i = 1 .. n)); ymin := min(seq(data[i, 2], i = 1 .. n)); xrange := xmax-xmin; xmin := xmin+(-1)*0.1e-1*xrange; xmax := xmax+0.1e-1*xrange; xrange := 1.02*xrange; yrange := ymax-ymin; ymin := ymin+(-1)*0.1e-1*yrange; ymax := ymax+0.1e-1*yrange; yrange := 1.02*yrange; dx := xrange/N; dy := yrange/N; res := Array(0 .. N-1, 0 .. N-1, 0); for i to n do ix := trunc((data[i, 1]-xmin)/dx); iy := trunc((data[i, 2]-ymin)/dy); res[ix, iy] := 1 end do; add(add(res[i, j], i = 0 .. N-1), j = 0 .. N-1) end proc:

``

bicationplot := proc (N) local Nr, Nt, x0, rmin, rmax, bif, k, ir, r, xn, i; global pts; Nr := 100; Nt := 200; x0 := .1; rmin := .75; rmax := 3.5; bif := Array(1 .. Nr*N, 1 .. 2); k := 1; for ir to Nr do r := rmin+ir*(rmax-rmin)/Nr; xn := x0; for i to Nt do xn := xn^2-r end do; for i to N do xn := xn^2-r; bif[k, 1] := r; bif[k, 2] := xn; k := k+1 end do end do; pts := bif end proc:

bif

bif

(1)

fractaldimension := proc (Noofitterations::integer, Npoints::integer, Nmax::integer) local res, xv, yv, line, stderrors, avgstderrors, i, avgline; avgstderrors := 0; avgline := 0; for i to Noofitterations do bicationplot(Npoints); res := [seq([1.0/n, boxcount(pts, n)], n = 1 .. Nmax, 10)]; xv := [seq(log(res[i][1]), i = 1 .. nops(res))]; yv := [seq(log(res[i][2]), i = 1 .. nops(res))]; line[i] := Fit(m*x+const, xv, yv, x, output = [leastsquaresfunction]); stderrors[i] := Fit(m*x+const, xv, yv, x, output = standarderrors) end do; for i to Noofitterations do avgstderrors := avgstderrors+stderrors[i] end do; avgstderrors := avgstderrors/Noofitterations; for i to Noofitterations do avgline := avgline+line[i] end do; avgline := avgline/Noofitterations; return FD = -(diff(avgline, x)), avgline, avgstderrors, loglogplot(res) end proc:

fractaldimension(10, 100, 100)

Error, (in boxcount) bad index into Array

 

``

 

Download First_part_fractal_determination_.mw

Explore the values of km digit(n,m) using km list for all m, 0 ≤ m ≤ 8.
Look at the output until you can make a conjecture that concerns the pattern
obtained for each fixed m, 0 ≤ m ≤ 8 using 

km := proc (n::posint, m::nonnegint)

local k,

mySum := 0;

for k to n do

mySum := mySum+k^m

end do;

return mySum

end proc

using a list km list(m,6,20) when m is not a multiple of 4, and km list(m,6,50) when m is a multiple of 4.

 

any help appreciated..THank you

 

First 160 161 162 163 164 165 166 Last Page 162 of 2097