Items tagged with field field Tagged Items Feed

I'm trying to calculate flux through a cone but gets the following:

Flux(VectorField(`<,>`(x, y, z), cartesian[x, y, z]), Surface(`<,>`(x,y,z), x^2+y^2<=2*a,z = 2*a-sqrt(x^2+y^2)))
Error, (in VectorCalculus:-Flux) error with surface input

Dear hope you will fine. I am try to make a program of square free factorization over a finite field whose alogrithm is below:

Algorithm: SFF (Square-Free Factorization)
  Input: A monic polynomial f in Fq[x]
Output: Square-free factorization of f

i←1; R ← 1; gf′;
  if g ≠ 0 then {
     cgcd(f, g);
     while w ≠ 1 do {
           ygcd(w, c); zw/y;
           RR·zi; i ← i+1; 
           wy; cc/y }
     if c ≠ 1 then {
           Output(R·SFF(c)p) }
     else  Output(R)
  else {
           Output(SFF(f)p) }

The attached file my try to make this, please find and help me to complete this. I am waiting your kind response.

With my best regards and sincerely.



Hope everything going fine with you. I have question

f := x^11+2*x^9+2*x^8+x^6+x^5+2*x^3+2*x^2+1

g := 2*x^10+x^7+2*x^4+x


if we take gcd of infinite field its answer is x^6+1 

and the GCD(3)[x]=Z_3[x] is given by x^9+2x^6+x^3+2

How we find GCD(3)[x] in maple. 

With my best regards and sincerely.



I am trying to obtain the splitting field of New_polyq. evala@AFactor did not complete. Applying splitting sequentially produced independent extensions from the first 2 (3?) factors. evala@Indep did not complete for the union of all 4 extensions.

What libraries would handle this better?

restart; _EnvExplicit:=false;interface(labelwidth=200);
Rho_polys:=rho[3,1]^3-2, rho[3,2]^2+rho[3,2]*rho[3,1]+rho[3,1]^2, 2*rho[6,1]^3+rho[6,1]^6-2, rho[12,1]^2+rho[6,1]^2-1, 2*rho[12,2]^2-rho[6,1]^2*rho[3,2]*rho[3,1]^2-2*rho[6,1]^2-2;

How to create polynomial ideals over algebraic extensions of the field of  rationals Q with Maple?
The Maple help to PolynomialIdeals
"All package commands support computations over the rational numbers, algebraic number fields, rational function fields, and algebraic function fields, as well as finite fields. Coefficients from algebraic extension fields can be specified using radicals or RootOfs"
is too poor. I also don't find any example on this topic in examples,PolynomialIdeals.

given a scalar field function f(x,y,z), how to plot it in maple?

I can't find a suitable command, if I missed, please tell me.

thank you for any advice.



I was trying to obtain the field equations for the QED Lagrangian and I was not sucessfull.


All my calculations are equal to 0.


Can someone give a hand?


Thanks a lot.

guys, is there any possibility to obtain field equations of einstein-hilbert action?

best regards

Hi there,
I have a set of differential equations whose solution, Jacobian matrix and its eigenvalues, direction field, phase portrait and nullclines, need to be computed.

Each of the equations has a varying parameter.

I know how to get the above for a single parameter value, but when I set a range of values for the parameters, Maple is not able to handle all cases as I would expect: solving the differential equation system:

eq1 := x*(1.6*(1-(1/100)*x)-phi*y)
eq2 := (x/(15+x)-0.3e-1*x-.4)*y+.6+theta
desys := [eq1, eq2];
vars := [x, y];
steadyStates := map2(eval, vars, [solve(desys)])

already yields an error:
Error, (in unknown) invalid input: Utilities:-SetEquations expects its 2nd argument, equations, to be of type set({boolean, algebraic, relation}), but received {-600*y+(Array(1..2, {(1) = 8400, (2) = 15900})), Array(1..5, {(1) = 0, (2) = 0, (3) = 0, (4) = 0, (5) = 0})}

The equations are the following:
de1 := diff(x(t), t) = x(t)*(1.6*(1-(1/100)*x(t))-phi*y(t));
de2 := diff(y(t), t) = (x(t)/(15+x(t))-0.3e-1*x(t)-.4)*y(t)+.6+theta

the parameters being:
phi:=[0 0.5 1 1.5 2]
theta:=[5. 10.]

How can I handle the situation so that Maple computes each of the above for each combination of the parameters?

I would like to avoid using two for loops and having to store all results in increasingly bigger and complicated arrays.

The worksheet at issue is this:


Hello everyone.

I have a vector field in 2d-cartesian coordinates which I would like to convert to a "normal" function, that is f(x, y) where when you put x and y in, you get the magnitude of the vector at that point.

Example vector field:

This one is very hard to handle by hand which is why I want to use Maple for it.

I tried VectorCalculus[Norm] but it gave me this:

The question of how to do matrix operations over a finite field with p^k elements k>1 was raised in this forum in April 27 2010. I wonder if there has been any progress since then for recent versions of Maple.

I am expecially interested in finding the rank of a matrix over a field of order p^k, k > 1.


Hello, I have a problem in plotting the vector field of dx/dt=sinx/sint

It shows me an error:

Error, (in DEtools/dfieldplot) unable to evaluate function `sin(x)` in evalhf

Dunno how to fix it

Here is my code:

DE := [diff(x(t), t) = sin(x)(t)/sin(t)];
dfieldplot(DE, x(t), t = -2 .. 2, x = -1 .. 2, arrows = medium, title = 'Diff', color = DC);

Hi, I am trying to plot a solution curve on a vector field with an initial condition of y(0)=0 and I keep getting error messages. This is what I have so far:


with(DEtools); dfieldplot(diff(y(x), x) = x^2+y(x)^2-1, y(x), x = 0 .. 5, y = -1 .. 5, arrows = line, title = 'Slope*Field');


Thank you. 

I am working over F_p, where p is prime. When p = 2, the number of 2x2x2 arrays (aka hypermatrices or tensors) over this field is 2^(2*2*2) = 256. Of those 256 arrays, I only want the symmetric ones: that is, if x_{ijk} is the ijk-th element of the array X, then x_{ijk} = x_{ikj} = x_{jik} = x_{jki} = x_{kij} = x_{kji}. Is there a quick loop that does this? 

I have plotted a predator/prey system using the following code and would now like to plot the same with the initial value problem of y=200, x=1000, with the end result being a graph of that one trajectory, without the phase field.  How can I do this?  Thanks.

DEtools[DEplot] ([diff(x(t), t) = 2*x(t)*(1-0.1e-3*x(t))-0.1e-1*x(t)*y(t), diff(y(t), t) = -.5*y(t)+0.1e-3*x(t)*y(t)], [x(t), y(t)], t=0 .. 10, x = 0 .. 11500, y = 0 .. 200, [[x(0) = 1000, y(0) = 200...

1 2 Page 1 of 2