I don't know where to ask this.  Sorry if this is the wrong forum.

I use maplev-mode

;; Version:    2.155
;; X-RCS:      $Id: maplev.el,v 1.14 2006-06-02 14:02:38 joe Exp $
 

I am attempting to modify it so
  f1 invokes maplev-cmaple-send-line
  f2 invokes maplev-cmaple-send-region
  f3 invokes maplev-cmaple-send-buffer
 

tried adding
(define-key map [(control c) (control c) ?l] 'maplev-cmaple-send-line)
(define-key map [(f1)]  'maplev-cmaple-send-line)

but, emacs seems to not recognize this.  I can change the global keymap to get this behavior.
Is there a way to do this in maplev.el?
 

Tom Dean

What can i do in this case? I,m traying solve a system of non linear equation like this, i want to know Rs and Rsh:

I have the following fuction in Laplace domain,

restart:with(plots):with(inttrans):

u:=Pi^4*s3^(alpha-1)/((s1^2+Pi^2)*(s2^2+Pi^2)*(-s1^2+Pi^2*s3^alpha-s2^2))-Pi*s1*s2^(alpha-1)/(s3*(1+s2^alpha)*(-s1^2+Pi^2*s3^alpha-s2^2))-Pi*s1*s2^(alpha-1)/(-s1^2+Pi^2*s3^alpha-s2^2);

Where, s1, s2, s3 are the Laplace variables.

x1:=invlaplace(u, s1, x);

This worked. But the next two doesn't work.

y1:=invlaplace(x1, s2, y);

uu:=invlaplace(y1, s3,t);

Even, I tried to plot the unevaluated invlaplace but no luck.

alpha:=1:t:=1:
plot3d(uu,x =1..2, y=1..2);
 

Am I missing something?

Dear all,

I have the following question, this code:

restart:
with(DifferentialGeometry):
DGsetup([w1,w2],N):
eq1 := ExteriorDerivative(w1);  
eq2 := ExteriorDerivative(w1) &wedge ExteriorDerivative(w2);
eq1 &wedge eq2;

Gives the error:
Error, (in DifferentialGeometry:-Tools:-DGzero)  given degree, 3, exceeds that of frame dimension, 2

Unfortunately, I am not so familiar with differential geometry but as far as I know dw1  (dw1  dw2) = 0 should be correct.

Thank you for your help
best
baustamm1

I am looking for persons with experience in setting up maple for the blind. The setting is a 21-year old blind
 student who wants to do math. He uses Mac OS/X and the built-in VoiceOver screen reader, plus refreshable Braille display (one-line), a Braille printer, and Nemeth, the Braille extension for math symbols.

Initially, we are working from no experience whatsoever, trying ideas, mostly failure. Getting some advice from persons who are blind and using maple would be wonderful. If anyone knows of an instance, a person who is blind and regularly uses maple, then we would both would be delighted to receive advice on what to install, how to set it up, and how to use maple on a daily basis.

One trouble for the VoiceOver screen reader is verbosity of maple output. Any methods to reduce the amount of output would be welcome. To uderstand the issue, look at the output of these two commands when interface(prettyprint=0) has been set to help out the VoiceOver screen reader:

matrix([[1,2],6,7]]); Matrix([[1,2],[6,7]]);

The first outputs the same text, which VoiceOver will read accurately. The second dumps out two lines of extraneous and confusing information, only 20% of which is useful. The verbosity is not restricted to matrix packages linalg or LinearAlgebra. It is a problem wtih nearly all packages.

-Grant Gustafson, Salt Lake City, Univ of Utah Math Dept

Dear friends:

I am facing two problems

1. one is to get solution of the below system of ODE for L=100 (highlited as red) and

2. the other is I want the graph in the form of solid line not poit, asterisk etc.

restart; epsilon := .1; Pr := 1; beta := .1; Sc := 1; S := 1; L := 20;
for i from -L while i <= L do;
a[i] := 1.0*i/L;
end do;
for i2 from -L while i2 <= L do;

fw := a[i2]; 

Eq1[i2] := eval(diff(F(eta), eta, eta, eta)+F(eta)*(diff(F(eta), eta, eta))-(diff(F(eta), eta))^2+S*(epsilon-(diff(F(eta), eta)))+epsilon^2);
Eq2[i2] := eval((diff(G(eta), eta, eta))/Pr-G(eta)*(diff(F(eta), eta))+F(eta)*(diff(G(eta), eta))); 
Eq3[i2] := eval(diff(H(eta), eta, eta)+Sc*(F(eta)*(diff(H(eta), eta))-beta*H(eta)));
IC[i2] := F(0) = a[i2], (D(F))(0) = 1, (D(F))(L) = epsilon, G(0) = 1, G(L) = 0, H(0) = 1, H(L) = 0;
dsys1[i2] := {Eq1[i2], Eq2[i2], Eq3[i2], IC[i2]};
dsol1[i2] := dsolve(dsys1[i2], numeric, output = listprocedure, range = 0 .. L);
dsol1x[i2] := subs(dsol1[i2], diff(F(eta), eta, eta));
dsol1y[i2] := subs(dsol1[i2], G(eta));
dsol1z[i2] := subs(dsol1[i2], H(eta)) end do;

for j from -L while j <= L do; 
g[j] := eval(-dsol1x[j](0)) end do;
with(plots); 

g6 := pointplot({seq([n/L, g[n]], n = -L .. L)}, symbol = asterisk, symbolsize = 15, color = red);
display(g6);

Please see the problem and correct as soon as possible. I am waiting your positive respone.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

It is well known that fsolve usually increases (internally) Digits in order to obtain the desired accuracy.

But in the following example, it seems that fsolve highly exaggerates :-)   

restart;
N:=40:
Digits:=100:
F:=expand(mul(x-k,k=1..N)):
f:=evalf(F):
S:=[fsolve(f,complex)];

Error, (in fsolve) Digits cannot exceed 38654705646

Note that the bug does not appear if e.g. F:=expand(mul(x-k-I, k=1..N)):

How to pdeplot Curl(A(x,y,z))=0 in maple

it has A1(x,y,z) , A2(x,y,z) and A3(x,y,z)

Hi all,

I have a partial differential equation similar to the following:

The solution is that f is identically equal to 0.

However, I am having trouble solving this equation in Maple. I type the following:

pde := diff(f(x, y), x)+diff(f(x, y), y) = f(x, y)+f(x, 0);

bv1 := f(x, 10) = 0;

bv2 := f(10, y) = 0;

solution := pdsolve(pde, {bv1, bv2}, numeric, time = x, range = 0 .. 10);

When Maple tries to evaluate the last expression, I get the error

Error, (in pdsolve/numeric/process_PDEs) PDEs can only contain dependent variables with direct dependence on the independent variables of the problem, got {f(x, 0)}

It seems to have difficulties with the expression "f(x,0)". Is there some trick to typing this in a way that makes Maple interpret it correctly?

Edit: I encounter the same problem, when I try to solve the ODE f'(x) = f(x) + f(0), where f(10) = 0.

Best regards.

Is there a method to relate groebner bases with monomials ideals

I attempt to understand is it possible whether to calculate the next expression by means of mathematica:

where

{,}- anti-communicator;

 a=0,1,2,3;

\tau^{0}-unity matrix; \tau^{i} - Pauli matrix;

and:
 

 
finally - Levi-Civita symbol; -Hermitian conjugation.

Thank you for your kind replies .

Is there a way to change the numeric formatting of 'Scientific' to use a cdot instead of a cross to represent multiplication?

I have a need to calculate base 2 math as would be done in an integrated circuit. Math will be done using 15 bit 2's complement mantissa and 8 bit exponent. (the exponent is always assumed to be negative) We need to perform multiplication and addition, where each internal operation is represented by such a number.

The goal is to best represent the errors associated with such calculations and export the resulting code to a hardware description language for implementation on an integrated circuit.

Is there a package in Maple that can do that? Any advice on how te proceed?

Hi all,

I am using Maple 2016.

I have defined 5 polynomials: f1, f2, f3, f4 and f5 with 5 unknowns q1,q2 ,q3, q4 and lamda.

After this, I generated the Gröbner basis. But when I try to find the normal set I got an error.

with(Groebner);

f1 := lamda*q1-(3380075947548081*q1*(1/140737488355328)-259050600068343*q2*(1/140737488355328)-1826834460600733*q3*(1/1125899906842624)+4414049272733425*q4*(1/9007199254740992))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f2 := lamda*q2+(259050600068343*q1*(1/140737488355328)+3380075947548081*q2*(1/140737488355328)-4414049272733425*q3*(1/9007199254740992)-1826834460600733*q4*(1/1125899906842624))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f3 := (1826834460600733*q1*(1/1125899906842624)-4414049272733425*q2*(1/9007199254740992)+843667886835955*q3*(1/35184372088832)-862655592804515*q4*(1/18014398509481984))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)))+lamda*q3;
f4 := lamda*q4-(4414049272733425*q1*(1/9007199254740992)+1826834460600733*q2*(1/1125899906842624)+862655592804515*q3*(1/18014398509481984)+843667886835955*q4*(1/35184372088832))*(q2*(8289619202186977*q1*(1/9007199254740992)+3380075947548081*q2*(1/281474976710656)-4414049272733425*q3*(1/18014398509481984)-1826834460600733*q4*(1/2251799813685248))+q3*(1826834460600733*q1*(1/2251799813685248)-4414049272733425*q2*(1/18014398509481984)+843667886835955*q3*(1/70368744177664)-215663898201129*q4*(1/9007199254740992))-q4*(4414049272733425*q1*(1/18014398509481984)+1826834460600733*q2*(1/2251799813685248)+431327796402257*q3*(1/18014398509481984)+843667886835955*q4*(1/70368744177664))-q1*(3380075947548081*q1*(1/281474976710656)-259050600068343*q2*(1/281474976710656)-1826834460600733*q3*(1/2251799813685248)+4414049272733425*q4*(1/18014398509481984)));
f5 := q1^2+q2^2+q3^2+q4^2-1;
ord := tdeg(q1, q2, q3, q4, lamda);
                  tdeg(q1, q2, q3, q4, lamda)
G := Basis([f1, f2, f3, f4, f5], ord);

IsZeroDimensional(G);
                             false
ns, rv := NormalSet(G, ord);
Error, (in Groebner:-NormalSet) The case of non-zero-dimensional varieties is not handled.

Any help please ?

Thank you.

Hello,

I was wondering if it is possible to use units in Maple so I can always check if the result I have at the end of calculation is the meter.  For example:

The answer is of course 3.10^8 m^3*kg/s^3

I try to do something with the units but I am unable to crreate something that will simplify the m/s ffactor to 1.

Any idea?

Thank you in advance for your help.