MaplePrimes Questions

I am trying to plot several curves with values from matrix columns on one plot.

 

L1 := zip(`[]`, groups, waitArray[ .. , 1]);

L2 := zip(`[]`, groups, waitArray[ .. , 2]);

plot([L1, L2])

plot(L1) and plot(L2) work fine.

But for a joined plot I have an error

 

Error, (in plot) incorrect first argument [Array(1..3, {(1) = [5, HFloat(0.2889190508261187)], (2) = [10, HFloat(0.9017773091123074)], (3) = [20, HFloat(2.4000000155705576)]}), Array(1..3, {(1) = [5, HFloat(0.019538366928174924)], (2) = [10, HFloat(0.06875266075705706)], (3) = [20, HFloat(0.1884304855100058)]})]

f1 := expand((a1*x^2+b1*x+c1)^n);
f2 := a1*x^2+b1*x+c1;

solve([coeff(f1, x, 4) = coeff(f2, x, 4),
coeff(f1, x, 3) = coeff(f2, x, 3),
coeff(f1, x, 2) = coeff(f2, x, 2),
coeff(f1, x, 1) = coeff(f2, x, 1),
coeff(f1, x, 0) = coeff(f2, x, 0)], [a1,b1,c1,d1]);

hi,

i cant solve this 

int(1/1-exp(a*x)*erfc(a*x))

restart;

Vector([a, b]);  <a, b>;  # Identical results

a:=[1, 2]:  b:=[3, 4]:

Vector([a, b]), <a, b>;   # Different results. Why?

 

 

 

I have a vector of lists, and each list is composed of indexed names.  For each list, I would like to assign all the indexed names in the list to the first variable of the list.  Right now, I'm trying to select the entries of the list using the op command and then assign one to the other, and get the following error:

Error, invalid left hand side in assignment

June_29.mw

Thanks!

how to convert decimal number into given decimal number like algebra

for example, convert 191.715 , given a=12.2, b=3.5

how to find this a^2 + b^3

Hello, I am a student typing up my homework assignments with maple, and I am takinng Inferential Statistics.  I can't figure out how to make the bar over a variable, so I can note the average.  It is just a line over a letter.  Kind of when notating a vector without the arrow.  It would be very helpful because I can't figure out how to make that symbol, it is not listed on the symbols list on the left side.

Thanks!

 

hi.i encountered this erroe  [Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system] with solving set of differential equation.please help me.thanks a lot  

dsys3 := {`1`*h1(theta)+`1`*(diff(h1(theta), theta, theta))+`1`*(diff(h2(theta), theta))+`1`*(diff(h2(theta), theta, theta, theta))+`1`*h3(theta)+`1`*(diff(h3(theta), theta, theta))+`1`*(diff(h1(theta), theta, theta, theta, theta)) = 0, `1`*h2(theta)+`1`*(diff(h2(theta), theta, theta, theta, theta))+`1`*(diff(h2(theta), theta, theta))+`1`*(diff(h1(theta), theta))+`1`*(diff(h1(theta), theta, theta, theta))+`1`*(diff(h3(theta), theta))+`1`*(diff(h3(theta), theta, theta, theta)) = 0, h3(theta)^5*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+(diff(h3(theta), theta, theta))*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+(diff(h3(theta), theta, theta, theta, theta))*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+h1(theta)*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+(diff(h1(theta), theta, theta))*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+(diff(h2(theta), theta))*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+(diff(h2(theta), theta, theta, theta))*h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+h3(theta)^4*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)+h3(theta)^4*(diff(h2(theta), theta, theta, theta, theta, theta, theta))*(`1`+ln(h3(theta))^2*`1`+2*ln(h3(theta))*`1`)-beta*h3(theta)^3*`1`-chi*ln(h3(theta))^2*`1`/kappa-chi*`1`/kappa-2*chi*ln(h3(theta))*`1`/kappa = 0, h1(0) = 0, h1(1) = 0, h2(0) = 0, h2(1) = 0, h3(0) = 1, h3(1) = 1, ((D@@1)(h1))(0) = 0, ((D@@1)(h1))(1) = 0, ((D@@1)(h2))(0) = 0, ((D@@1)(h2))(1) = 0, ((D@@1)(h3))(0) = 0, ((D@@1)(h3))(1) = 0, ((D@@2)(h3))(0) = 0, ((D@@2)(h3))(1) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 600, numeric, output = listprocedure);
%;
Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

 

Let be given the complex number z  satisfying condition abs(z+3-2I)=3. I want to find the set of points representing the complex number w, knowing that w - z = 1 +3I. I tried

Restart:
assume(a::real, b::real,x::real, y::real);
z:=x+y*I;
w:=a+b*I;
abs(w-1-3*I+3-2*I)=3;

 

Edit

Restart:
assume(a::real, b::real,x::real, y::real);
Set:=abs(z+3-2*I):
w:=x+y*I;
sol:=solve(w - z =1+3*I,{z});
z1:=subs(sol,Set);
A:=abs(z1);
map(x->x^2,A=3);

I do have some m files and they containing procedures where i am interested in
How can i read those procedures in a m file

Thanks in advance

Hi Maple friends. :)

I would like to graph and solve: 2*sin(2 x)=1 for the interval 0<=x<=360. My textbook gives the answers of x=15, 75, 195 and 225.

But I get an error messages:

plot(2*sin(2*x) = 1, x = 0 .. 360);
Error, invalid input: plot expects its 1st argument, p, to be of type {array, list, rtable, set, algebraic, procedure, And(`module`, appliable)}, but received 2*sin(2*x) = 1

solve(2*sin(2*x) = 1, x = 0 .. 360);
Error, invalid input: too many and/or wrong type of arguments passed to solve; first unused argument is x = 0 .. 360

Any help would be appreciated! Thanks in advance.

 

x: =Matrix([[a1,a2],[a3,a4]])

after some calculation,

assign(%)

a1 etc have value,

how to make a1,a2,a3,a4 back to variable in maple 12?

How to calculate floor(10^(10^(10^(10^(10^(-10^10))))))? My simpleminded try is

Error, numeric exception: underflow
PS. I think that is about 10^(10^10).

m := Matrix([[a1,a2,a3],[a4,a5,a6],[a7,a8,a9]]);

m2 := Determinant(m-Matrix([[1,0,0],[0,1,0],[0,0,1]]));

1. which family of polynomials do m2 belong to?

2. how to analyze m2?

What´s the error here??????

Solução do Sistema de EDO's por Transformada de Laplace

 

restart

with(inttrans):

eq1 := diff(x(t), t)+(2.2*x(t)*y(t)+0.5e-1*x(t)/(5+0.5e-1*t)) = 0;

diff(x(t), t)+2.2*x(t)*y(t)+0.5e-1*x(t)/(5+0.5e-1*t) = 0

(1)

eq2 := diff(y(t), t)+(2.2*x(t)*y(t)-(0.5e-1*(0.25e-1-y(t)))/(5+0.5e-1*t)) = 0;

diff(y(t), t)+2.2*x(t)*y(t)-0.5e-1*(0.25e-1-y(t))/(5+0.5e-1*t) = 0

(2)

eq3 := diff(z(t), t)-2.2*x(t)*y(t)+0.5e-1*z(t)/(5+0.5e-1*t) = 0;

diff(z(t), t)-2.2*x(t)*y(t)+0.5e-1*z(t)/(5+0.5e-1*t) = 0

(3)

EQ := [eq1, eq2, eq3]:

for i to 3 do La[i] := laplace(EQ[i], t, s) end do;

s*laplace(x(t), t, s)-1.*x(0.)+2.200000000*laplace(x(t)*y(t), t, s)+laplace(x(t)/(100.+t)^1., t, s) = 0.

 

-1.*y(0.)-0.2500000000e-1*(exp(100.*s))^1.*Ei(1., 100.*s)^1.+1.*s^1.*laplace(y(t), t, s)^1.+2.200000000*laplace(x(t)^1.*y(t)^1., t, s)+1.*laplace(y(t)^1./(100.+1.*t)^1., t, s) = 0.

 

s*laplace(z(t), t, s)-1.*z(0.)-2.200000000*laplace(x(t)*y(t), t, s)+laplace(z(t)/(100.+t)^1., t, s) = 0.

(4)

LL := subs({laplace(x(t), t, s) = X, laplace(y(t), t, s) = Y, laplace(z(t), t, s) = Z}, [La[1], La[2], La[3]]);

[s*X-1.*x(0.)+2.200000000*laplace(x(t)*y(t), t, s)+laplace(x(t)/(100.+t)^1., t, s) = 0., -1.*y(0.)-0.2500000000e-1*(exp(100.*s))^1.*Ei(1., 100.*s)^1.+1.*s^1.*Y^1.+2.200000000*laplace(x(t)*y(t), t, s)+1.*laplace(y(t)/(100.+t)^1., t, s) = 0., s*Z-1.*z(0.)-2.200000000*laplace(x(t)*y(t), t, s)+laplace(z(t)/(100.+t)^1., t, s) = 0.]

(5)

sol := solve(LL, [X, Y, Z]):

assign(sol):

SOLS[X, Y, Z]:

SOLT := map(invlaplace, [X, Y, Z], s, t);

[-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(x(_U1)/(100.+_U1), _U1 = 0. .. t))+x(0), -1.*(int(y(_U1)/(100.+_U1), _U1 = 0. .. t))-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))+y(0)+0.2500000000e-1*ln(1.+0.1000000000e-1*t), 2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(z(_U1)/(100.+_U1), _U1 = 0. .. t))+z(0)]

(6)

SOLTT := evalf(subs({x(0) = 0.5e-1, y(0) = 0, z(0) = 0}, SOLT));

[-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(x(_U1)/(100.+_U1), _U1 = 0. .. t))+0.5e-1, -1.*(int(y(_U1)/(100.+_U1), _U1 = 0. .. t))-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))+0.2500000000e-1*ln(1.+0.1000000000e-1*t), 2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(z(_U1)/(100.+_U1), _U1 = 0. .. t))]

(7)

xx := evalc(Re(SOLTT[1]));

-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(x(_U1)/(100.+_U1), _U1 = 0. .. t))+0.5e-1

(8)

yy := evalc(Re(SOLTT[2]));

0.2500000000e-1*ln(abs(1.+0.1000000000e-1*t))-1.*(int(y(_U1)/(100.+_U1), _U1 = 0. .. t))-2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))

(9)

zz := evalc(Re(SOLTT[3]));

2.200000000*(int(x(_U1)*y(_U1), _U1 = 0. .. t))-1.*(int(z(_U1)/(100.+_U1), _U1 = 0. .. t))

(10)

plot([xx, yy, zz], t = 0 .. 500, legend = [x, y, z]);

Warning, expecting only range variable t in expression -2.200000000*int(x(_U1)*y(_U1),_U1 = 0. .. t)-1.*int(x(_U1)/(100.+_U1),_U1 = 0. .. t)+.5e-1 to be plotted but found names [_U1, x, y]

 

 

NULL

NULL


Download laplace.mw

First 1413 1414 1415 1416 1417 1418 1419 Last Page 1415 of 2434