I need help to create a program that will find all the positive integers n, where n < 1000, such that
(n 􏰀-1)!= 􏰁 􏰀-1 (mod n^2 ) . program has to be in full and state the values of n obtained. 

Try solve 

restart; with(PDEtools);
U := diff_table(u(x, y, z));

pde[1] := x*y*U[z]+x*U[x]+2*y*U[y] = 0;

bc[1] := eval(U[], z = 0) = x^2+y^2;
sys[1] := [pde[1], bc[1]];
pdsolve(sys[1])

Maple nothing returned. Where my mistake ?? Thank you.

restart;
eq1 := T - M*g = 0;
eq2 := T1*sin(theta1) + T2*sin(theta2) - T = 0;
eq3 := -T1*cos(theta1) + T2*cos(theta2) = 0;
soln := solve({eq1,eq2,eq3},{T,T1,T2});

limit(subs(theta2=theta1,soln[3]),theta1=0);

This returns a limit -> ... expression.  Since there is a sin(theta1) in the denominator, this should be undefined or infinity.

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

Hey everyone!

My question is pretty straightforward: I have a second degree tensor and I need a fourth degree order one from it.

In Python, there is a command that performs the action, would there be one in Maple? Or any kind of process to go through to get the result?

Thanks

My code :

> restart; with*PDEtools;

> PDE := diff(u(x, t), `$`(t, 2)) = 4*(diff(u(x, t), `$`(x, 2)));

> pdsolve({PDE, u(x, 0) = cos((1/2)*x), (D[1](u))(x, 0) = 0}, u(x, t));
 

And i get Error, (in pdsolve/BC) invalid input: indets expects 1 or 2 arguments, but received 3.Please help me find my Error.Thank you!!!

I would like to see the list of metrics recognize by Maple with their acornym.  For example,:

>Setup(coordinates = spherical, metric = kerr)

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:

Dear friends

I have a long-running code that sometimes takes an hour to complete. I wonder is there a way to find out what line of code is currently running in a long-running maple code or to show an update of the variables. I have some print commands but are shown just after completion not before. 

Thank you in advance for your time.

Ok, I'm  a confident Matlab user, but for some reason Maple befuddles me.

Let's say I have two functions, x=f(t) and y=g(t) and I want to fill both two vectors with values, with t=-n:n.  How do I construct a loop to do that?

There is a horse a buggy ride around a small village which takes roughly 30 minutes.  Here is an example timing for 12 consecutive rides [34, 29, 32, 32, 28, 28, 27, 28, 39, 24, 27, 27].

How can I create a monte carlo simulation graph that would estimate the future times based on given data?  Do I randomly pick numbers from the given list for a simulation or generate random numbers based on mean and standard deviation generated from the data?

When would the best possible time to come back after 4 rides be?

How could i show wilsons theorom on maple?

(p-1)!=-1(modp) if and only if p is prime.

Hello,

I have a question in factoring a matrix of complex polynomials into its Smith Normal Form. It seems that Maple is not giving me correct answers.

I tried a simple matrix:

A:=Matrix([[I*z-1, z+I], [z+I, z+I]]);

Maple gives me:

SmithForm( A, z, method='rational'  ,output=['S'] );

which is:

Matrix(2, 2, {(1, 1) = 1, (1, 2) = 0, (2, 1) = 0, (2, 2) = -4*z^2-(3*I)*z^3+z^4+1+(3*I)*z})

Since the gcd of the 4 elements is (z+i), it seems the Smith Normal Form should not be in this form. Also, I got errors in solving the other 2 matrices:

SmithForm( A, z, method='rational'  ,output=['U'] );
Error, (in convert/radical) numeric exception: division by zero

How to make it work for complex polynomials?

Sorry that I am not an expert in Algebra~

Thanks a lot!

William

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?

Hi I have the question where i have to create a program in Maple

to find all the solutions to x^2 = -1(mod p) where 0 <= x < p . 

The progam has to be tested with different p values.