Hi all,

I'm new to Maple. Probably this is trivial, but how can I solve an inequality as e.g., the following

-2 <=  a/x <= 0

knowing that x > 0?

Thanks

S.

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; g ← f′;
  if g ≠ 0 then {
     c ← gcd(f, g);
     w ← f/c;
     while w ≠ 1 do {
           y ← gcd(w, c); z ← w/y;
           R ← R·zi; i ← i+1; 

           w ← y; c ← c/y }
     if c ≠ 1 then {
           c ← c1/p;

           Output(R·SFF(c)p) }

     else  Output(R)
  else {
           f ← f1/p;

           Output(SFF(f)p) }

  end.

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

Help.mw

With my best regards and sincerely.

I want convert the integro differential equation to volterra integral equation by integration both side of integro-diff then subs in  initial condition 

intE := diff(u(x), x) = 1+x-x^2+int((x-t)*u(t), t = 0 .. x), u(0) = 3

How do it?

I need to replace the "MISSING" in MultiDiff := proc(f,n,x) ---MISSING--- end proc; so that MultiDiff(f,n,x) evaluates as f (n) (x) if f is a suitably differentiable function, n is a positive integer, and x is an undefined variable..F(n) being an nth derivitave..Not sure where to start..Any help appreciated

I was hoping to browse color results using a slider and text box. I see the "use DocumentTools" command, but was not able to get "use ColorTools" to allow execution of the "Color()" function. 

My slider "n" has the following code;

use DocumentTools in

Do(%outn=%n);
Do(%outrgb=my_rgb(%n));
use ColorTools in
Do(%outcolor=color_my_rgb(%n));
end use;
end use;

How equating the coefficient in the same power of x  

diff(u(x), x) = 1+int(u(t), t = 0 .. x), u(0) = 0

intE1 := subs({u(t) = sum(a[i]*t^i, i = 0 .. 3), u(x) = sum(a[i]*x^i, i = 0 .. 3)}, intE);
intE2 := eval(intE1);
intE3 := collect(intE2, x);
ans := zip(`=`, [coeffs(lhs(intE3)-rhs(intE3), [x, x^2, x^3, x^4], 't')], [seq(0, i = 1 .. nops([t]))]);

Shows me wrong I also want to find values of a[1] and a[2] and a[3] and a[4] without using my hand and from the initial condition a[0]=0?

Assume that I have r:= -6x+3y+23x²-4xyz+7z² . By using coeffs(r,x,'k') I can find the coefficients of 1, x, and x². 

What should I write to get the conditions that make the coefficient of x zero? How can I just pick the coefficient of x, and solve it. 

I have several coefficients but the answer of this question will help me. 

Hello,

I use a map() command to get the values of function for each element of vector. Here is an example of a simple task:

restart:

A:=[1,2,8];

f1:=x->2*x;
f1_table_A:=Vector([map(f1,A)]);

I would like to do the same with piecewise function but it doesn't work. Here is an example:

restart:
with(plots):
f:=2*x:
g:=x^2:

h:=x->piecewise( 0<=x and x<= 5, f, 5<x and x <= 10, g ) ;
A:=[1,2,8]:
h_table_A:=map(h,A):

Is it possible to work with piecewise function? Maybe someone has an idea how to do it in different way?

Kind regards,

Iza

Hi all,

I need to solve det[M]=0 for omega.

M is:

M := Matrix(8, 8, {(1, 1) = BesselJ(n, tp*a), (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = -BesselJ(n, tg*a), (1, 6) = -BesselY(n, tg*a), (1, 7) = 0, (1, 8) = 0, (2, 1) = k*n*BesselJ(n, tp*a)/(tp^2*a), (2, 2) = I*k*n*mu0*omega*(diff(BesselJ(n, tp*a), a))/(tp^2*a), (2, 3) = 0, (2, 4) = 0, (2, 5) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (2, 6) = -k*n*BesselY(n, tg*a)/(tg^2*a), (2, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (2, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (3, 1) = 0, (3, 2) = BesselJ(n, tp*a), (3, 3) = 0, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (3, 7) = -BesselJ(n, tg*a), (3, 8) = -BesselY(n, tg*a), (4, 1) = -I*omega*`&varepsilon;p`*(diff(BesselJ(n, tp*a), a))/tp^2, (4, 2) = k*n*BesselJ(n, tp*a)/(tp^2*a), (4, 3) = 0, (4, 4) = 0, (4, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*a), a))/tg^2, (4, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*a), a))/tg^2, (4, 7) = -k*n*BesselJ(n, tg*a)/(tg^2*a), (4, 8) = -k*n*BesselY(n, tg*a)/(tg^2*a), (5, 1) = 0, (5, 2) = 0, (5, 3) = k*n*BesselY(n, t0*b)/(t0^2*b), (5, 4) = I*mu0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (5, 5) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (5, 6) = -k*n*BesselY(n, tg*b)/(tg^2*b), (5, 7) = -I*`&mu;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (5, 8) = -I*`&mu;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (6, 1) = 0, (6, 2) = 0, (6, 3) = BesselY(n, t0*b), (6, 4) = 0, (6, 5) = -BesselJ(n, tg*b), (6, 6) = -BesselY(n, tg*b), (6, 7) = 0, (6, 8) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = BesselY(n, t0*b), (7, 5) = 0, (7, 6) = 0, (7, 7) = -BesselJ(n, tg*b), (7, 8) = -BesselY(n, tg*b), (8, 1) = 0, (8, 2) = 0, (8, 3) = -I*epsilon0*omega*(diff(BesselY(n, t0*b), b))/t0^2, (8, 4) = -k*n*BesselY(n, t0*b)/(t0^2*b), (8, 5) = I*`&varepsilon;g`*omega*(diff(BesselJ(n, tg*b), b))/tg^2, (8, 6) = I*`&varepsilon;g`*omega*(diff(BesselY(n, tg*b), b))/tg^2, (8, 7) = -k*n*BesselJ(n, tg*b)/(tg^2*b), (8, 8) = -k*n*BesselY(n, tg*b)/(tg^2*b)});

and n=0,1,2.

Except of omega and k ,other parameters is canstant.

After using
with(LinearAlgebra):
detM := Determinant(M):

I used solve(detM=0,omega) and fsolve() but it dosnt work. how can i solve it?

Thanks alot.

can anyone please help me to find roots of (2*cos(0.5*x)*sin(0.5*x)*cos(3.775*x)+2.2075*((cos(0.5*x))^2)*sin(3.775*x)-0.453*((sin(0.5*x))^2)*sin(3.775*x)=0 ?

i type on maple like this:

solve(2*cos(0.5*x)*sin(0.5*x)*cos(3.775*x)+2.2075*((cos(0.5*x))^2)*sin(3.775*x)-0.453*((sin(0.5*x))^2)*sin(3.775*x)=0);

but it said the solutions may have been lost

thanks before

Working on a code to create random lines. Here is the code: 

M:=1; N:=1500; R3:=rand(1..3): # M=lines, N=steps

for i from 1 to M do
X[i,0]:=0;
Y[i,0]:=0;
   for j from 0 to N do
   if j=1 then X[i,1]:=1; Y[i,1]:=0;
    elif j>1 then
      r:=R3();
        if r=1 then X[i,j]:=2*X[i,j-1]-X[i,j-2]; Y[i,j]:=2*Y[i,j-1]-Y[i,j-2]
        elif r=2 then X[i,j]:=X[i,j-1]+Y[i,j-1]-Y[i,j-2]; Y[i,j]:=Y[i,j-1]-X[i,j-1]+X[i,j-2];
        elif r=3 then X[i,j]:=X[i,j-1]-Y[i,j-1]+Y[i,j-2]; Y[i,j]:=Y[i,j-1]+X[i,j-1]-X[i,j-2];
      end if;
    end if;
   R[i,j]:=[X[i,j],Y[i,j]];
   K[i,j]:=[X[i,j],Y[i,j]];
  end do:
end do:

Now the code works fine. But I don't want the lines to cross its own path. I want it to stop if R[i,j]=R[i,k] for any 0<k<i-2 for j=0..N. What I have so far is:

for i from 1 to M do
  for j from 1 to 40 do
    R[i,j];
    for k from 1 to j-2 do
     if R[i,j]=K[i,k] then x:=j; print(b[i]=j); break; end if;
     L[i,j]:=R[i,j];
     end do;
   end do;
end do;

Now it works aswell, but it doesn't stop when it hits the crossing point. I want it to stop so i get one b[i]=j per M. So I need it to break two for loops so that it goes through i M times and stops everytime it reaches a crossing point

Hello,

I am trying to differentiate a matrix containing four variables, alpha, alphaB (representing alpha bar), beta, and betaB (for beta bar). They are are variables with respect to t. I then need to let t=0 and then simplify the result with some initial conitions i have. Could you tell me any useful maple functions which i can use to do this? If you need anymore information let me know and thank you for helping me.

Robbie

How to use The  value of lambda in other operation ?

solve({alpha = (1/2)*lambda*alpha*Pi, beta = (1/2)*lambda*beta*Pi}, lambda)

for exampel 

lambda^6;

we have a curve (x^2+y^2)^2=x^2-y^2

parametric representation of the curve is given by

x(t)=cos(t)/1+sin^2(t)

y(t)=sin(t)cos(t)/1+sin^2(t)

and we are asked to plot the curve using a parametric plot for -Pi<t<Pi

Dear

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.