Following is the set of commands I performed.

> points := plots[implicitplot](t+sin(8*x^3*t)-x^3 = 0, t = -0.50e-1 .. 0.5e-1, x = -10/9 .. 10/9);
print(??); # input placeholder
> ptslist := convert([1, 1], points, listlist);

In an output I need the data at the specified values of t. It will be great favor, if someone can help me out.

Thanks

Dear mapleprimes users,

I have a problem concerning this function:

Naief := proc(A::integer,B::integer, p::posint)
local x, a:= A mod p, b:= B mod p;

 for x from 1 to p-1 do
    if a^x mod p = b then return x end if
end do;
print(¨No solution.¨)

end proc

It works fine for what it should do, finding x for a^x = b mod p by inserting x from 1 to p-1 until it finds an
apprioperate x.

My problem is concerning its computation time, which I like to calculate with:

Codetools:- Usage(Naief(a,b,p), output = [cputime], quiet, iterations= 2^12)

The problem is that keeps repeating

alot even when a,b and p are immensely huge. I don`t know how to fix it, because I need real CPUtimes which increase
a,b and p increase. I get the idea that the values of the CPUtime are not realistic.

Thanks for the help!

I'm trying to verify whether the result obtained from a recursive Legendre function that I defined is the same as that obtained from the function available in Maple.

The arguments for LegendreP() mentioned seem to be v,x where v is the algebraic expression's degree and x is the expression. I tried calling the function as LegendreP(6,x) since I am to obtain the value at n=6 but this statement is not giving me an output in return that I can use to compare.

Would I have to plot both functions in order to be able to compare them?

I've got a pair of equations :

x^3-4x=y and y^3-4y=x

 which I've defined as eqns:={,}

and 9 solutions as solns1:={,,..}

and being stored as s1,s2,..s9

when I run a command such as testeq(subs(s1,eqns[1])=subs(s1,eqns[2])

I get an error of passing invalid arguments into testeq. What I essentially need to show is that on substituting for x,y from each s1,..s9; both equations get the same result. What am I doing incorrectly?

I've also noticed that just subs(s1,eqns[1]) returns an equality; I don't quite understand why

Dear mapleprimes users,

I have a problem with this function:

B is the length of my array

A := Array(1..B,0);

for i from 1 to B do
p = nextprime(i)

a = primroot(p)

A[i] := A[i] + convert((Usage(mlog(b,a,p,method=indcalc), output=[cputime,output],quiet)[1],decimal,15))
end do;

A;

My idea is to create an array A with all the CPUtimes from the 1st prime til the Bth prime, and then plot it.

But my problem is that I get an error; Illegal use of an object as a name.

I don`t know what to do could you help me?

Thanks!

I'm given the following two equations:

x^3-4x=y, y^3-4y=x

to solve the system, I've just used

eqns:={x^3-4x=y,y^3-4y=x};

vars:={x,y};

solns:=solve(eqns,vars);

and have obtained only four solutions when I should instead get 9. Is there a mistake in my approach?

Here's my problem,

I'm this line of code

result:=solve({x+y=1,x+2y=4},{x,y})
and the value of result is
{x=-2,y=3} and the value of result[1] is x=-2.

I would like the value of result to be {-2,3} for I need to use those output later. Is there any way to put those output in a list?

The following works as intended:

restart;
diff(x^7, x$i):
seq(%, i=0..7);

Combining the two commands into one, however, does not work:

restart;
seq(diff(x^7, x$i), i=0..7);

          Error, invalid input: diff expects 2 or more arguments, but received 1

How does one explane this?  I was unable to find the reason by looking at the help page for seq().

I want to get a step-by-step solution of following integral:

Int(.5*exp(-s*t)*exp(-4.5*t)*sin(2*Pi*t), t)

how I can achieve this? Is there a Task for this? I tried Student Package but failed.

how to demonstrate the first quantization and second quantization of maxwell equation with physics package

How can this error be corrected '' error, (in fsolve) fsolve cannot solve on 0=0 ''. See the worksheet p4.mw

Thanks.

we always have subscript variable in the math book, but how could this be natral done in maple

I want to get a seq a₁ a₂ a₃

seq(a[i],i=1..3)

but how could I get a  a_ij

seq(a[i_j],i=1..3);

and

seq(a[ij],i=1..3);  both was not right

Hello, I need a person who can help me with task from numerical methods.

I'm trying to solve some ODE analitically. But Maple gives me an incorrect solution. What am I doing wrong? Thank you.

Hello Guys,

i´m quite new to Maple and found a guide to draw an Ossannacircle with Maple.
The problem is, that it doesn´t work and i don´t get what´s wrong.

Could someone help?

That´s the code:

restart: with(plots):
interface(imaginaryunit=j, displayprecision=4):
plotsetup("inline",plotoutput=terminal,plotoptions="colour=cmyk,resolution=2000"):
U[n]:= 500: # V
U[2,n]:= 83.2: # V
I[1,n]:= 28.7: # A
phi[n]:= arccos(0.844):
evalf(180/Pi*phi[n]); # Grad
n[0]:= 1000: # 1/min
R1:= 0.465: R2:= 0.00985: # Ohm
ü:= U[n]/U[2,n]; # Übersetzungsverhältnis
I[1,0]:= 8.5: # A
phi[0]:= arccos(0.15):
evalf(180/Pi*phi[0]); # phi[0] in Grad
U[k]:= 170.5: # V
I[k]:= 28.7: # A
phi[k]:= arccos(0.277):
evalf(180/Pi*phi[k]); # phi[k] in Grad
I_[1,n]:= ?I[1,n]*exp(j*phi[n])?;
Pn_:= I_[1,n]:
I_[1,0]:= ?I[1,0]*exp(j*phi[0])?; P0_:= I_[1,0]:
I[1,k]:= ?I[1,n]*U[n]/U[k]?;
I_[1,k]:= 'I[1,k]*exp(j*phi[k])';
Pk_:= I_[1,k]:
G1:= '(Im(P0_)-xm)^2+(Re(P0_)-ym)^2 = r^2';
G2:= '(Im(Pk_)-xm)^2+(Re(Pk_)-ym)^2 = r^2';
G3:= '(Im(Pn_)-xm)^2+(Re(Pn_)-ym)^2 = r^2';
Loe:= solve({G1,G2,G3}, [xm,ym,r]);
assign(Loe);
linie:= proc(ap, ep, opt)
# ap... Anfangspunkt, ep...Endpunkt
# opt...Optionen (z.B. Farbe, Dicke); als Liste vorzugeben
plot([[Im(ap),Re(ap)],[Im(ep),Re(ep)]],op(opt));
end proc:
punkt:= proc(pp)
# pp... Zeiger auf den Punkt
pointplot([[Im(pp), Re(pp)]], symbol=solidcircle):
end proc:
kreis:= plottools[circle]([xm,ym], r, color=black):
P0:= punkt(P0_): Pk:= punkt(Pk_): Pn:= punkt(Pn_):
D_:= 'Re(P0_)+j*Im(Pk_)';
DD:= punkt(D_): # D ist für den Differentialoperator reserviert
P0D:= linie(P0_, D_, [color=black]):
DPk:= linie(D_, Pk_, [color=black]):
I_[2,k]:= 'I_[1,k]-I_[1,0]';
D_E:= 'Re(I_[2,k])*R[1]/(R[1]+R[2]*ü^2)';
E_:= 'D_ + D_E';
EE_:= 'E_ - P0_';
E:= punkt(E_):
ML:= linie(P0_, E_, [legend="Drehmomentlinie ", color=blue]):
LL:= linie(P0_, Pk_, [linestyle=dash, legend="Leistungslinie ",color=blue]):
I_10:= linie(0, P0_, [color=blue]):
d:= 0.05*Im(Pk_): # Textabstand
text:= textplot([[Im(P0_)-0.7*d, Re(P0_)+d/2,"P0"],[Im(Pn_)-0.7*d, Re(Pn_)+d/2,"Pn"], [Im(Pk_)+d, Re(Pk_),"Pk"],[Im(D_)+d/2, Re(D_)+1,"D"], [Im(E_)+d/2, Re(E_)+1,"E"]]):
KD:= display(kreis, P0, Pk, Pn, DD, E, I_10, LL, ML, P0D, DPk, text,scaling=constrained, title=typeset("Asynchronmotor\n U =",U[n], " V, I = ",I[1,n]," A\n"), titlefont=[TIMES,12,BOLD],labels=[typeset(Im('I[1]/A')),typeset(Re('I[1]/A'))]): KD;
G4:= ''(x-xm)*(xk-xm) + (y-ym)*(yk-ym) = r^2'';
G5:= isolate(G4,y);
G6:= collect(G5, x);
AA:= op(rhs(G6));
A := AA[1]/x;
G7:= isolate(A = Re(EE_)/Im(EE_), yk);
G8:= (xk-xm)^2 + (yk-ym)^2 = r^2;
Loe2:= solve({G7,G8}, [xk,yk]);
I_[1,kipp]:= subs(Loe2[1], j*xk + yk); Pkipp_:= I_[1,kipp]:
Pkipp:= punkt(Pkipp_):
LMS:= sqrt(3)*U[n]/1000; # Leistungsmaßstab in kW/A
DMS:= sqrt(3)*U[n]/(2*Pi*n[0]/60); # Drehmomentmaßstab in Nm/A
moment:= proc(P_)
local S; global P0_, E_, D_, DMS;
# Ermittlung der Strecke von P bis zur Drehmomentlinie
S:= Re(P_)-Re(P0_)-(Im(P_)-Im(P0_))*(Re(E_)-Re(D_))/
(Im(D_)- Im(P0_));
evalf(S*DMS);
end proc:
leistung:= proc(P_)
local S; global P0_, Pk_, D_, LMS;
# Ermittlung der Strecke von P bis zur Leistungslinie
S:= Re(P_)-Re(P0_)-(Im(P_)-Im(P0_))*(Re(Pk_)-Re(D_))/(Im(D_)-Im(P0_));
evalf(S*LMS);
end proc:
M[n]:= moment(Pn_);
Pmech[n]:= leistung(Pn_);
M[kipp]:= moment(Pkipp_);
Pmech[kipp]:= leistung(Pkipp_);
text2:= textplot([[Im(Pkipp_), Re(Pkipp_)+d, "Pkipp"]]):
plotsetup(ps, plotoutput="print", plotoptions="noborder,resolution=2000"):
display(KD, Pkipp, text2, labels=[" "," "], tickmarks=[0,0],caption=" ");
ZL:= readstat("Zeigerlänge vom Punkt (0,0) bis Pk in mm:");
IMSd:= abs(Pk_)/ZL; # A/mm
LMSd:= evalf(IMSd*U[n]*sqrt(3)/1000); # kW/mm
DMSd:= evalf(LMSd*1000/(2*Pi*n[0]/60)); # Nm/mm

Regards
Alex