Hi, i want to solve the Ax=b with jacobi method and 4 Iteration, A=Matrix(3, 3, [1, 2, 1, 2, 1, 3, -1, -3, 4]) and b=[8, 13, 5]. 

i want to write a code with for or if or while. tnx

Dear sir / madam

Hi
I'm a student of mechanical engineering and looking for some HPM ( homotopy perturbation method) maple codes with the related article to practice more in this field.
can I ask you for help?
unfortunately, I couldn't find any articles with maple solution about HPM.

sincerely
Aidin

When running this code in Maple 13 and 16 this works just fine,

GenMs:= (k::posint)-> assign(m||(1..2*k) =~ seq([dx||i, dy||i][], i= 1..k)):
GenMs(k), m || (1 .. 2*k);
gln := evalDG(l1 &s n1+`&s`(l2, n2)+sum('epsilon || i'*(cat(m, 2*i-1) &s cat(m, 2*i-1)+cat(m, 2*i) &s cat(m, 2*i)), i = 1 .. k));

But running the same code in Maple 2018 I get the following error:

Error, (in DifferentialGeometry:-Tensor:-SymmetrizeIndices) expected 1st argument to be a tensor. Received: `m2*i-1`^2

I don't know if it's related to Maple 2018 screwing with my tensor formatting:
https://gyazo.com/88eaeceda2e36cf411df44dbc4aa3ab6

(compare https://gyazo.com/b33396131165d66bbcf16e45d20cc579)

For whatever reason I just noticed that my inline &t's have been turned into `&t` but that's a separate issue.

For the most part, when I have come across bugs (or what appears to be bugs from my naive view, maybe it is left unevaluated for some means of diplomacy that I am not privy to) they are very difficult or impossible for me to find a way to get around.

But I have found alot recently that I was able to fix very quickly using functions that are already inbuilt in maple. Most have to do  with evaluating functions of trancendental arguements or composites of trancendental functions, but still, it leaves me feeling a little disturbed just how readily it was resolved.

I don't know exactly what to say, perhaps a 'dispatch' needed to be included in an inbuilt procedure to inform the 'engine' that a user has come across a bug of some kind at that address, so that maplesoft can keep the process of applying patches in house, avoiding having to enable the evolution/personal development the inept user that clearly stumbled upon the bug by sheer coincidence whist entering formula and observing output in a relatively sub human cognitive state?

I only just noticed the Start.mw file being always number 1 in my Recent document list.  I don't recall it being there, ever, and I don't think it should.  Why is it there now?  Anyone else have this?

Ok the error message I originally recieved in requesting 'expression' rather than 'value' was telling me the input is ambigous, isn't exactly helpful.

With the uploaded worksheet as last attempted, it gives the error of something about an invalid sequence, but still and I honestly do not see how when MathML is an XML application, surely whatever is put into the math container is retrievable considering XML applications are based on just handling the rendering of symbols.

Anyway I can only post here seeings that following the link for error messages rarely opens an actual help page in the brower, it simply tells you there is no help page for that error.

Download 26052018.mw

I'm trying to define a monomial order - Monomials are first compared by their weighted degree, with ties broken by pure lexiocraphic order (plex). I want the monomial order just to find the LeadingTerm of polynomials.

Maple has wdeg(W,V) as a built in monomial order, which compares monomials first by their weights, but breaks ties using reverse lexicographic order (tdeg).

One way I could plausibly do this is using the 'matrix' constructor, and defining the plex order by a succession of weights, but this seems like overkill, and I would have a n*n matrix, where I'd like n to reach over 1000.

Is anyone aware of a simple way of defining such a monomial order?

I just noticed something funny. If I use the left side tool of Maple to define a matrix with zero entries and then change some entries and take determinant, then Maple treats zeros as symbold!!!! If I have to define my matrix completely myself, then what is the use of this tool?

Dear Friends, 

I would appreciate your help in resolving some issues. Let me describe my dummy code and the issues I am having. 

I want to know the proper value of two parameters beta and `ΔA` in a tripe integral function  by the NonlinearFit command.The triple integral function is complex.

The code of function is below:

int(tan(beta)^2*exp(-Pi*tan(beta)^2*((x-varepsilon)^2+(0-varsigma)^2)/eta^2)/eta^2, [eta = 22.83-sqrt((5.83+`ΔA`)^2-varepsilon^2) .. 22.83+sqrt((5.83+`ΔA`)^2-varepsilon^2), varepsilon = -5.83-`ΔA` .. 5.83+`ΔA`, varsigma = -1 .. 1])-(int(tan(beta)^2*exp(-Pi*tan(beta)^2*((x-varepsilon)^2+(0-varsigma)^2)/eta^2)/eta^2, [eta = 22.83-sqrt(5.83^2-varepsilon^2) .. 22.83+sqrt(5.83^2-varepsilon^2), varepsilon = -5.83 .. 5.83, varsigma = -1 .. 1])) 

And this is the complete code:

datax := [-8, -4.5, -.5, 4.5, 8, 11.5, 14.5];

datay := [0.287e-2, 0.266e-2, 0.259e-2, 0.199e-2, 0.164e-2, 0.113e-2, 0.78e-3];

f := NonlinearFit(int(tan(beta)^2*exp(-Pi*tan(beta)^2*((x-varepsilon)^2+(0-varsigma)^2)/eta^2)/eta^2, [eta = 22.83-sqrt((5.83+`ΔA`)^2-varepsilon^2) .. 22.83+sqrt((5.83+`ΔA`)^2-varepsilon^2), varepsilon = -5.83-`ΔA` .. 5.83+`ΔA`, varsigma = -1 .. 1])-(int(tan(beta)^2*exp(-Pi*tan(beta)^2*((x-varepsilon)^2+(0-varsigma)^2)/eta^2)/eta^2, [eta = 22.83-sqrt(5.83^2-varepsilon^2) .. 22.83+sqrt(5.83^2-varepsilon^2), varepsilon = -5.83 .. 5.83, varsigma = -1 .. 1])), datax, datay, x)

I try to run it and get the value of parameters beta and `ΔA` ,but I keep getting this error,

Error, (in Statistics:-NonlinearFit) integration range or variable must be specified in the second argument, got HFloat(1.0) = HFloat(16.073603031200726) .. HFloat(29.58639696879927)
Does someone how to deal with this problem?

I would sincerely appreciate any inputs in this regard. 

Ok so as title says that is what i need, or the process of changing the document property listing that is generated when i start a new worksheet. 

The reason being is that of course the obvious solution is for me to have a personalized template for a worksheet with the desired attribute listing, which i have, but the trouble is i have already generated a large number of archive entries in a database or "library" of worksheets that i am building that have the standard attribute listing, and for the new functionality that i have chosen i rely on having the specific custom attribute listing for the new template, and i want to implement this functionality into use for the existing database entries as well as the future ones created, without having to copy and paste things for each of the existing entries for the obvious reason that this is boring and annoying

I using Maple 18 (not Maple 2018) and I'm trying to figure out how to grab earthquake data from earthquakescanada database from here http://www.earthquakescanada.nrcan.gc.ca//stndon/NEDB-BNDS/bull-en.php using the HTTP requests.

First I used the default search within the web browser, and get a new address which I enter as the URL

HTTP:-Get("http://www.earthquakescanada.nrcan.gc.ca//stndon/NEDB-BNDS/bull-en.php?time_start=2018%2F04%2F23+22%3A47%3A00&time_end=2018%2F05%2F23+22%3A47%3A00&depth_min=0&depth_max=100&mag_min=-3&mag_max=9.9&shape_type=region&radius_center_lat=50&radius_center_lon=-95&radius_radius=1000&region_north=90&region_south=41&region_east=-40&region_west=-150&eq_type_L=1&display_list=1&list_sort=date&list_order=a&tpl_output=html&submited=1"

It takes a long time to download the information and would require HTML surgery but changing the option for output to txt or csv, it's faster and in a much more readable form.  However it's not in a table or Array format, it has become a string.

Is there any way to use ImportMatrix, or ImportData to get a better format of the information?  - both give errors in Maple18.  Or am I stuck trying to use string surgery in Maple 18?  The Import command isn't available until Maple 2016 (I don't mean the Import command within ExcelTools) and I believe that works in Maple 2018 however I'm at a loss for trying to use it in Maple 18. 

I need some help. I'm trying to solve this system of equations, but maple says the solutions may have been lost.

I don't know why. Here are the equations:

I have four equations,very unknown variables in the equation.

I am trying to solve for any 4 unknowns ,not must had be Zoo1.Zoo2.th1.th2,it can be Za1.Za2.Zb1.Zb2.  

Any help  would be greatly appreciated.

    /  /                                                    0.5
    |  |/         2                                        \   
 int\x \\D(c[2]) x  + 2 c[2] D(x) x + D(c[1]) x + c[1] D(x)/   

                                              \            \
            2             2                1.5|            |
    + c[2] x  + c[1] x - x  - 1.504505556 x   /, x = 0 .. 1/
   /   /                                                    0.5
   | 2 |/         2                                        \   
int\x  \\D(c[2]) x  + 2 c[2] D(x) x + D(c[1]) x + c[1] D(x)/   

                                             \            \
           2             2                1.5|            |
   + c[2] x  + c[1] x - x  - 1.504505556 x   /, x = 0 .. 1/
 

Hi,

I naively thought that
f := t -> Some expression containing t
and
f := unapply(Some expression containing t, t)
where two different ways to define the same function f.

Recently I defined f and g this way
X := RandomVariable(Normal(m, s)):
U := Sample(X, 30000) ;  
f := t -> PDF(X, t);
g := unapply(PDF(X, t), t);

Could you explain me why f~(U) exectues in about 10 seconds and g~(U) in less than 0.1 s ?

Thanks in advance

 

Hellow , Any body can help me to write the codes for the differential equation to solve by Galerkin finite element method by four elements by taking step size h=pi/2 and compare with homotropy method by plotting the graph of Galerkin method, homotropy method and exact, here is my codes

ode := diff(y(x), x, x) + y(x)-1=0:
ics := y(0) = 0, y(2*pi) = 0:
dsolve({ics, ode}):