Maple 2017 Questions and Posts

These are Posts and Questions associated with the product, Maple 2017

How can I draw Steiner trees of hypercubes in graph theory? 

hi

I have a linear system with varibles trying to plot 3d the solutions x, y, z

here is my code: linear_var.mw

please any comment might help.

*******************************

restart;

Omega:=10:N:=0.5:M:=sqrt(N(N+1)):

a11:=0.5*(1+2*N)+M*cos(phi):
a12:=-0.5*((1+theta)^3+(1-theta)^3):a13:=-0.5*(N+M*cos(phi))*((1+theta)^3-(1-theta)^3): a21:=M*sin(phi): a22:=(-(1+2*N)+0.5*M*cos(phi))*((1+theta)^3+(1-theta)^3): a23:=-(Omega+0.5*((1+theta)^3-(1-theta)^3)*M*sin(phi)): a31:=0.25*((1+theta)^3-(1-theta)^3): a32:=Omega: a33:=-0.5-(N+0.25)*((1+theta)^3+(1-theta)^3): b1:=-0.5*a31: b2:=0: b3:=0.25+((1+theta)^3+(1-theta)^3)/8:

slove([a11*x+a12*y+a13*z=b1,a21*x+a22*y+a23*z=b2,a31*x+a32*y+a33*z=b3[,[x,y,z]);
Error, unable to match delimiters
Typesetting:-mambiguous(Typesetting:-mambiguous(slovelparlsqba11

  sdotx + a12sdoty + a13sdotzequalsb1commaa21sdotx + a22sdoty + 

  a23sdotzequalsb2commaa31sdotx + a32sdoty + a33sdotzequalsb3lsqb

  comma(xyz)rparsemi, 

  Typesetting:-merror("unable to match delimiters")))

plot3d(x, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed);
plot3d(y, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed); plot3d(z, theta = .1 .. 5, phi = 0 .. 2*Pi, axes = boxed);

I have a guess about the set of the zeros of the following polynomial

y(1-x^{m+1}z)+(1-x^{n+1}z), (here m,n are positive integers and z is a primitive d-root of unity)

which are located on the complex 2-dimensional torus. The set of solutions is finite (I think the system is zero-dimensional). My goal is to verify my guess numerically using Maple for some small values of m and n and a fixed value of z. I think if (x,y) is a solution, then x is either a (n-m) root of unity or a (n-m) root of 1/z^2 (where n>m).

You can find my code for n=3 and m = 1 attached (I was not able to load the mw format so I put the zip version). I consider z to be a third root of unity but actually, I am interested in putting z= exp(2pi/3*i) and even the real third root z=1 is not interesting for me, but since the exponential representation led to an error, I changed it to z and mentioned that z^3=1. Still, it has an error and  I would be grateful if you could let me know how I can correct this code.  

Question.maple.zip

Is there a way to delete frames arround the legend in a plot?

Why does this not return a quick solution?

sol := solve(.2894606222 = .4090662520/((1+0.1481289280e-2*x)^5.034901366*(1/(1+0.1481289280e-2*x)^5.034901366)^.1986136226), x)

mw.mw

Are Maple settings responsible for this error:

HHow do i correct "error, illegal use of an object as a name"

Hello everyone,

I am trying to solve a system of 4 ODEs that describe a physical system, however, I get the message:

Error, (in dsolve/numeric/bvp) matrix is singular

Any help is much appreciated. The initial conditions are such that to find solution in steady state, that means, v(0)=v(period) and so on with all variables.

Download LLCcircuit_SteadyState_numerical.mw

Hello everyone,

I want to solve a simple ODE numerically having a piecewise periodic function in the ODE.

The periodic piecewise function definition is taken from this  post.

When solving using "dsolve" numerical, I get this error:

Error, (in f) unable to store '100000*square_wave_periodic' when datatype=float[8]

restart;
d := .5;
fs := 100*10^3;
Vin := 10;
C := 10*10^(-6);
R := 10;
L := 10*10^(-6);

DEFINE PERIODIC FUNCTION
square_wave := proc (t) options operator, arrow; piecewise(0 <= t and t <= d/fs, Vin, d/fs < t and t < 1/fs, 0) end proc

square_wave_periodic := proc (t) options operator, arrow; square_wave(t-floor(t*fs)/fs) end proc

ode1 := L*(diff(i(t), t))+v(t) = square_wave_periodic
ode2 := C*(diff(v(t), t))+v(t)/R = i(t)
IC := v(0) = 0, i(0) = 0
sol := dsolve({IC, ode1, ode2}, {i(t), v(t)}, numeric)

LCRcircuit_transient_numerical.mw

Any help is much appreciated ! Thank you !

I tryning to to plot 3d and then contour plot of the function N vs alpha

This is my try please any comments might help

N_vs_alpha.mw

Hello everyone,

I am trying to solve a simple nonlinear 4 variable function with NLPsolve giving real values ranges for the variables but the NLPsolve gives "non numerical" result. Find attached the worksheet.CoreOptimization.mw

Am I doing something wrong?

Regards,

Cata

Hhow can i solve R it says here root of

Calculate.mws

How do you stop Maple during this code execution:

proc()
x:
print(1);
goto(x);
end():

Hello everyone,

I am trying to solve an equation using the "solve(...)" function but the coefficients of the equation are changing.

I am using the coefficients as list but the "solve(...)" function gives error: "Error, (in assuming) when calling 'Engine:-Dispatch'. Received: 'badly formed input to solve: not fully algebraic'"

I am attaching the Maple document.

 

However, what I want to achieve, is the solution for each value in the list.

Here is an easy example of what I am saying:

a=3*x+k where k=[list of values], a=[another list of values].

I want to find "x" for each pair of values (a,k) in the list.

 

Thank you !
 

``

``

LCC Resonant Converter - Frequency Finder

 

 

restart``

with(Student[Calculus1])

V__out := 50.0

`&omega;__line` := 2.0*Pi*50

t := [seq(0 .. 10*10^(-3), .1*10^(-3))]

theta := `&omega;__line`*t

v__line := 85*sqrt(2.0)*`~`[sin](theta)

m := 1

Q__s := 3.2*`~`[`^`](`~`[sin](theta), 2)

``

f__n := `assuming`([solve((1/120)*V__out = 1/sqrt((-m*x^2+m+1)^2+(Q__s*(x-1/x))^2), x, useassumptions)], [x::positive])

Error, (in assuming) when calling 'Engine:-Dispatch'. Received: 'badly formed input to solve: not fully algebraic'

 

``


 

Download LCC_gain_freqFinder.mwLCC_gain_freqFinder.mw

H all expert

first this Expression 

ni := diff(Q(x, t), t)+a*Q*(x, t)*(diff(Q(x, t), x))+b*(diff(Q(x, t), t$3))+d*(diff(Q(x, t), t$5)) = 0

then I want to solve  diff(Q(x, t), t$5) 

diff(Q(x, t), t$5) = solve(ni, diff(Q(x, t), t$5))#eq2

p (x, t) :=H1(t)*Q(x*H2(t), H3(t)) #assumption 

k := diff(p(x, t), t)+a*p(x, t)*(diff(p(x, t), x))+diff(p(x, t), t$3)+d*(diff(p(x, t), t$5))+c*p(x, t) #eq3

r := subs(diff(Q(x, t), t$5), k) #subs eq2 in eq3

i recived error

SUBS1.mw

1 2 3 4 5 6 7 Last Page 1 of 40