MaplePrimes Questions

Hi,

See the code mli.mws

What is the mistake in this command ?

chfonc:=cat('piecewise(' , seq(cat(convert(t<lf[k],string) , ' , 1 , ' ,convert(t<lf[k],string), ' , 1 , ' ) , k=1..p) , '1) ') :

I am trying to solve three non-linear coupled initial value problems. The equations are a bit complicated. The MAPLE code is 

restart

with*plots

eq[1] := (1/6)*Ra*v(x)^3*((2*v(x)*w(x)^2*(v(x)^2/w(x)^2-3)*(diff(w(x), x))-6*w(x)^3*(1-v(x)/w(x))^2*(diff(v(x), x)))*A+v(x)^2*w(x)^2*(2*v(x)/w(x)-3)*(diff(w(x), x))-3*w(x)^4*(1-v(x)/w(x))^2*(diff(v(x), x)))/(w(x)^2*(w(x)+2*A)^2*(2*v(x)+3*B))+u(x)

eq[2] := (1/105)*(diff(u(x)^2*(v(x)^2+11*B*v(x)+39*B^2...

The following simple code

myproc1:=proc({parlist::list:=[]})
return parlist[1];
end proc:

myproc2:=proc(parlist::list)
return myproc1('parlist'=parlist);
end proc:

myproc2([1,2]);

fails with


Error, (in myproc1) invalid subscript selector

The reason is that myproc1 has the optional argument parlist, which has the same name as myproc2. myproc2 seems to make...

Hi all,

I want to compare two signals where intersections between the two signals represent the instant of command (like boolean pulse) of a thyristor. The two signals are volatges and command of thyristor is boolean. A conversion of voltage to boolean signal should also done, I think ?

Thanks in advance

See image below.

First of, thank you commenters, you have helped me so much in the past!

I am doing a matrix inversion with random numbers and I would like to loop the process, and have each inverstion go into a matrix, such that the matrix will be a bunch of rows ( i number of rows) and three columns (the matrix inversion yeieds three values)

 

> for i to 2 do;
>
> a1h := Generate(float(range = 0.1316881e-2 .. 0.1359706e-2)); a1f := Generate(float(range = .317149811 .. .327463626...

I'm am new to maple and am not very familiar with it.

I have a state space represnetaiton (Xdot = AX+Bu  y+CX+Du)of a 3rd order non linear system.

X=[X1; X2; X3]

Xdot=[X1dot; X2dot; X3dot]

A=[0 -1 -1; 1 0.2 0; X3 0 -2.4]

B=[0; 0; 2]

C=[1 0 0]

D=0

With initial conditons:

X1(0)=X2(0)=x3(0)=An arbitrary concept of my choice.

I need to plot x1(t) Vs time.

I also need to plot the state trajectories...

Hi,

I am really stuck in one more bit thats the least square problem.  leastsquare1.mw 

Hope someone can help me to solve this bit and get me out of this problem :(.

Least square method is successfuly implemented in MATLAB code (WRITTEN IN THIS MESSAGE AT THE END), but I couldn't do it in Maple.  I am having problem in the for loop section,  I am attaching Matlab code and  maple code ...

I have a differential that I'm solving using rk4 on maple and I want to plot the resaulting points in an interval. How do I do that?

here is what I've already done:

> f := {y(0) = 0, diff(y(x), x) = cos(5*y(x))-x};

> p := dsolve(f, y(x), numeric, method = classical[rk4], stepsize = .1);
proc(x_classical)  ...  end;

> for i from 0 by .1 to 3 do p(i) end do;
           ...

scatter2d  quantile2 scatter1d is removed?


with(stats): with(statplots): with(transform):
with(random): with(describe): with(plots):
datax:=[normald[5,1](100)]:
datay:=[normald[7,3](100)]:
scatter2d(datax,datay);
quantile2(datax,datay);
scatter1d(datax);
scatter1d[jittered](datax);

who(things) replace it?

i have a systems of differential equation:

J[n]=2*A[n]

S[n+1]=J[n]

A[n+1]=S[n] + 0.5*A[n]

with initial condition A[0]=10, J[0]=30, S[0]=5

i have a had time writing a procedure for this systems. help please

do i need a for command for this procedure

Can you divide 2 polynomials with Real coeffs.   eg 1.87x3-2x2+5.789x- 3,2 by    1.5 x2+2x-1.4 to get the quotient and remainder? If so how?

hello
my duty variational iteration method for solving nonlinear equations with maple code need finding.
I looked but I could not google
please help me

Hello,

Suppose that S1,S2,S3 are 3 blocks of a matrix S := Matrix(...). I would like to do

 

out := ArrayTools:-Concatenate(2,convert(S1,vector),convert(S2,vector),convert(S3,vector)):

but in the case when I have n blocks.

 

Thanks,

Dimitar

given : A(n+1)=2*A(n) , A(0)=70

write a procedure and find A(5)=?

> dsol1 := dsolve({-(diff(y[1](t), t, t))*sqrt(y[1](t)^2+y[2](t)^2)^3/y[1](t) = 0, -(diff(y[2](t), t, t))*sqrt(y[1](t)^2+y[2](t)^2)^3/y[1](t) = 0, y[1](0) = 1, y[2](0) = 0, (D(y[1]))(0) = 0, (D(y[2]))(0) = 1});                             dsol1:= >

Why the result can't execute?

First 1745 1746 1747 1748 1749 1750 1751 Last Page 1747 of 2428