Where is the problem within this procedure which I can not get the numerical values of V[1] , V[2], ...

Download ProcPropl.mw

fig([L1(color = blue), L2(color = blue), L3(color = green), C1(color = black), C2(color = black), C3(color = orange), C4(color = orange)]);
point(oo, [0, 0]);
                               oo

rotation(fig1, fig, oo, Pi/2, 'counterclockwise');
Error, (in geometry:-rotation) wrong type of arguments why thos error ? Thank you;
 

q := (f)->local ff := `if`(type(f,procedure),f,`x`->f); ff(z + (ff(k))):

Trying to create something so that I can do q(x) or q(x->x^2) or whatever.

Trying to understand the difference between an expression and a function as represented by maple. I thought I could just turn whatever is passed to q, such as x*cos(x) in to a function by x->x*cos(x) as one would do normally but this does not work.

How_to_collect_Coefficent.mw

Kindly help me to find the coefficients from the equation. I have attached the MAPLE file here.

Thank You 

Is there any easy way to do subs([a,b,c]=x, [2,4,5,1,a,4,5,a,b,434,a,b,c,5,5,5,1]) = [2,4,5,1,a,4,5,a,b,434,x,5,5,5,1]

and another way to get [2,4,5,1,a,4,5,a,b,434,x,b,c,5,5,5,1]

The first subsitutes the entire match and the second substitutes when matches but only for first element.

I know I can write routines to do this and will after posting this but is there a native way that is simple?

Hi,

I try to enrich the layout of my documents by tables like that of trigonometric numbers. Is there a way to make it simpler? without going through the encoding of all values.Thanks

TableauTest.mw

The worksheet below displays a sphere with a surface patch defined by the arcs of three intersecting circles on the sphere's surface.

How can the patch be colored differently from the non-patch sphere's surface?

Patch_on_a_Sphere.mw

I am looping over a list of functions to calculate various properties of the said functions. some of these calculations take too long. I was wondering if there was a way to skip to the next element of the loop if a calculation takes more than a specific time.

 I looking to contact a maple sales representative. I am from india I am not getting any contact number. Can anyone help.

Kind help anyone if possible. I do research in mathematics  field but at a personal level. I need help on which product to buy and a sales represitive number or someone to contact me so that i can first understand certain doubts from my personal end then I can opt for.

I may have asked this before, but if so lost answer.

I'm trying to put a relatively elaborate peace of code in a loop, but can't figure out how to have multiple lines. Each command separated by ';' is awkward and seems to cause me to make even more typos then normal and it's hard to read the code.

Is there a way to have command on separate lines with in a loop?

-Arthur 

I used Maple to factor (using ifactor) the number 10^128 + 1.  Maple produced the factors in a matter of seconds.

One of the factors had 116 digits.  I then restarted Maple and executed the command ifactor(X) where X was the

factor with 116 digits.  Maple found the number was prime instantly.  I do not understand how one can determine

that a number  with 116 digits is prime in less than a second.  My guess is that Maple has stored some large primes.

Also I used  Maple to factor (10^128+3) and I turned the program off because it was taking forever.

Download HP-Numeric.mw

The matrix data[gridji,4] in my code has float undefined. Please, how do I handle it?

In part of my program, I have a loop that the iterative calculations are performed there. During running the program when it comes to the loop the memory usage become full in a short period of time and I face the memory shortage. Since my program is so long and huge, I only can show the loop part of it here, you may look at where is the problem of high memory usage in loop and how to deal with it. Any comment will be appreciated.

I again note that, this segment can not be run independently since it is a part of big program. I put it just to see and mmake comment if any.

Download HighMemUseProb.mw

Hello everyone,

I want to use "proc" in a program. Although the original code is good for me, I modify it by "proc" to be convenient. Thank you very much!

With regards

clearall;
restart;
Digits := 15;
## Initial polygon
x[1] := [-1, 1/2, 1, 0, -sqrt(2)/2, -1, 1/2, 1, 0, -sqrt(2)/2, -1, 1/2, 1, 0, -sqrt(2)/2];
y[1] := [0, sqrt(3)/2, 0, -1, -sqrt(2)/2, 0, sqrt(3)/2, 0, -1, -sqrt(2)/2, 0, sqrt(3)/2, 0, -1, -sqrt(2)/2];
s2 := [-1, 1/2, 1, 0, -sqrt(2)/2, -1];
t2 := [0, sqrt(3)/2, 0, -1, -sqrt(2)/2, 0];
assign(a[-6] = 13/1296, a[-5] = -11/648, a[-4] = -1/16, a[-3] = -107/1296, a[-2] = 179/1296, a[-1] = 9/16, a[0] = 137/144, a[1] = 137/144, a[2] = 9/16, a[3] = 179/1296, a[4] = -107/1296, a[5] = -1/16, a[6] = -11/648, a[7] = 13/1296);

L := 3;
N := numelems(x[1]);
## Perform k-1 iterative steps
# L--the number of iterations; nk--the number of refined points after k steps;
for k to L do
    nk := 3^(k - 1)*(N - 6) + 6;
    for i from 3 to nk - 2 do
        x[k + 1][3*i - 8] := a[-6]*x[k][i + 2] + a[-3]*x[k][i + 1] + a[0]*x[k][i] + a[3]*x[k][i - 1] + a[6]*x[k][i - 2];
        y[k + 1][3*i - 8] := a[-6]*y[k][i + 2] + a[-3]*y[k][i + 1] + a[0]*y[k][i] + a[3]*y[k][i - 1] + a[6]*y[k][i - 2];
        x[k + 1][3*i - 7] := a[-5]*x[k][i + 2] + a[-2]*x[k][i + 1] + a[1]*x[k][i] + a[4]*x[k][i - 1] + a[7]*x[k][i - 2];
        y[k + 1][3*i - 7] := a[-5]*y[k][i + 2] + a[-2]*y[k][i + 1] + a[1]*y[k][i] + a[4]*y[k][i - 1] + a[7]*y[k][i - 2];
        x[k + 1][3*i - 6] := a[-4]*x[k][i + 2] + a[-1]*x[k][i + 1] + a[2]*x[k][i] + a[5]*x[k][i - 1];
        y[k + 1][3*i - 6] := a[-4]*y[k][i + 2] + a[-1]*y[k][i + 1] + a[2]*y[k][i] + a[5]*y[k][i - 1];
    end do;
end do;
## Plot the result
s1 := evalf(simplify(convert(x[L+1], list))):t1 := evalf(simplify(convert(y[L+1], list))):
f1 := <<s1> | <t1>>: f2 := <<s2> | <t2>>:
plot([f1, f2], linestyle = [1, 3], color = [black, red]);

T5Scheme:=proc(x[1],y[1],L)          
local i,k,nk,result,N,a[-1],a[7],a[6],a[5],a[4],a[3],a[2],a[1],a[0],a[-2],a[-3],a[-4],a[-5];          N:=numelems(x[1]):         
   for k from 1 to L do              
 nk := 3^(k - 1)*(N - 6) + 6;         
     for i from 3 to nk - 2 do            
      x[k + 1][3*i - 8] := a[-6]*x[k][i + 2] + a[-3]*x[k][i + 1] + a[0]*x[k][i] + a[3]*x[k][i - 1] + a[6]*x[k][i - 2];                       y[k + 1][3*i - 8] := a[-6]*y[k][i + 2] + a[-3]*y[k][i + 1] + a[0]*y[k][i] + a[3]*y[k][i - 1] + a[6]*y[k][i - 2];                  x[k + 1][3*i - 7] := a[-5]*x[k][i + 2] + a[-2]*x[k][i + 1] + a[1]*x[k][i] + a[4]*x[k][i - 1] + a[7]*x[k][i - 2];                  y[k + 1][3*i - 7] := a[-5]*y[k][i + 2] + a[-2]*y[k][i + 1] + a[1]*y[k][i] + a[4]*y[k][i - 1] + a[7]*y[k][i - 2];                   x[k + 1][3*i - 6] := a[-4]*x[k][i + 2] + a[-1]*x[k][i + 1] + a[2]*x[k][i] + a[5]*x[k][i - 1];               
    y[k + 1][3*i - 6] := a[-4]*y[k][i + 2] + a[-1]*y[k][i + 1] + a[2]*y[k][i] + a[5]*y[k][i - 1];                   result[1]:=x[L+1]:           result[2]:=y[L+1]            
  end do;          
     end do;  
     return result  
  end proc: