Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hello,

Could you please help me with the following problem? I'm new to Maple and i need some help.

Solve the equation x^3 - a*x + 1 = 0 , in x. Determine the particular solution for a=1,2,... .Graphically represent the polynom that appears in the equation, in a case where the equation has a real root and in a case where the equation has 3 real roots.

Thank you !

Hi

I d like to limit my solution to a real (non complex) solution
There should be a simple solution to my calculation but maple can t process the solutions in one of my "solve" commands.

How can do I tell maple to limit itself to one
 

NULL

restart

with(Student[Calculus1]):


#geometry [mm]

b := 250:

h := 720:

ds := 70:

d := h-ds:

As := 3000:


#concrete [MPa]

fck := 30:

fcm := fck+8;

38

(1)

Ecm := 33000;

33000

(2)

`ϵc1` := 2.2*(1/1000);

0.2200000000e-2

(3)

eta := `ϵc`/`ϵc1`;

454.5454545*`ϵc`

(4)

Ec1 := fck/`ϵc1`;

13636.36364

(5)

k := 1.05*Ecm/Ec1;

2.540999999

(6)

sigma := fcm*(-eta^2+eta*k)/(1+(k-2)*eta);

38*(-206611.5702*`ϵc`^2+1154.999999*`ϵc`)/(1+245.9090904*`ϵc`)

(7)


#steel [MPa]

Es := 200000:

fsy := 400:

fsu := 600:

`ϵy` := fsy/Es;

1/500

(8)

`ϵsh` := 0.9e-2:

`ϵsu` := 0.75e-1:

P := 4:

`ϵs` := `ϵcm`*(d-c)/c;

`ϵcm`*(650-c)/c

(9)

i := 1;

1

(10)

for `ϵcm` from .1*(1/1000) by .1*(1/1000) to 10*(1/1000) do `ϵs` := `ϵcm`*(d-c)/c; T[1] := `ϵs`*Es*As; T[2] := fsy*As; T[3] := (fsu+(fsy-fsu)*((`ϵsu`-`ϵs`)/(`ϵsu`-`ϵsh`))^P)*As; C := b*c*(int(sigma, `ϵc` = 0 .. `ϵcm`))/`ϵcm`; eq[1] := T[1] = C; `ϵl`[1] := `ϵy`; eq[2] := T[2] = C; `ϵl`[2] := `ϵsh`; eq[3] := T[3] = C; cc := max(solve(eq[1], c)); `ϵss` := subs(c = cc, `ϵs`); Ta := subs(c = cc, T[1]); if `ϵss` >= `ϵl`[1] then cc := max(solve(eq[2], c)); Ta := subs(c = cc, T[2]); `ϵss` := subs(c = cc, `ϵs`) end if; if `ϵss` >= `ϵl`[2] then cc := max(`assuming`([solve(eq[3], c, useassumptions)], [c::real])); Ta := subs(c = cc, T[3]); `ϵss` := subs(c = cc, `ϵs`) end if; M[i] := b*cc^2*fcm*(int(sigma*`ϵc`, `ϵc` = 0 .. `ϵcm`))*10^(-6)/`ϵcm`^2+T*(d-cc)*10^(-6); phi[i] := `ϵcm`/cc; cd[i] := cc/d; print(`ϵcm`, `ϵss`, Ta/As); i := i+1 end do

0.1000000000e-3, 0.1955232439e-3, 39.10464877

 

0.2000000000e-3, 0.3845102290e-3, 76.90204580

 

0.3000000000e-3, 0.5671741821e-3, 113.4348364

 

0.4000000000e-3, 0.7437144096e-3, 148.7428819

 

0.5000000000e-3, 0.9143174400e-3, 182.8634880

 

0.6000000000e-3, 0.1079158043e-2, 215.8316087

 

0.7000000000e-3, 0.1238400148e-2, 247.6800297

 

0.8000000000e-3, 0.1392197667e-2, 278.4395334

 

0.9000000000e-3, 0.1540695238e-2, 308.1390476

 

0.1000000000e-2, 0.1684028897e-2, 336.8057793

 

0.1100000000e-2, 0.1822326682e-2, 364.4653363

 

0.1200000000e-2, 0.1955709188e-2, 391.1418377

 

0.1300000000e-2, 0.2226921078e-2, 400

 

0.1400000000e-2, 0.2583745185e-2, 400

 

0.1500000000e-2, 0.2954159196e-2, 400

 

0.1600000000e-2, 0.3336139462e-2, 400

 

0.1700000000e-2, 0.3727768551e-2, 400

 

0.1800000000e-2, 0.4127227927e-2, 400

 

0.1900000000e-2, 0.4532791280e-2, 400

 

0.2000000000e-2, 0.4942818378e-2, 400

 

0.2100000000e-2, 0.5355749454e-2, 400

 

0.2200000000e-2, 0.5770100024e-2, 400

 

0.2300000000e-2, 0.6184456101e-2, 400

 

0.2400000000e-2, 0.6597469817e-2, 400

 

0.2500000000e-2, 0.7007855358e-2, 400

 

0.2600000000e-2, 0.7414385220e-2, 400

 

0.2700000000e-2, 0.7815886744e-2, 400

 

0.2800000000e-2, 0.8211238846e-2, 400

 

0.2900000000e-2, 0.8599369121e-2, 400

 

0.3000000000e-2, 0.8979250971e-2, 400

 

Error, complex argument to max/min

 

`ϵcm` := 0.35e-2; `ϵs` := `ϵcm`*(d-c)/c; T[1] := `ϵs`*Es*As; T[2] := fsy*As; T[3] := (fsu+(fsy-fsu)*((`ϵsu`-`ϵs`)/(`ϵsu`-`ϵsh`))^P)*As; C := b*c*(int(sigma, `ϵc` = 0 .. `ϵcm`))/`ϵcm`; eq[1] := T[1] = C; `ϵl`[1] := `ϵy`; eq[2] := T[2] = C; `ϵl`[2] := `ϵsh`; eq[3] := T[3] = C; cc1 := max(solve(eq[1], c)); `ϵss1` := subs(c = cc, `ϵs`); Ta1 := subs(c = cc, T[1]); cc2 := max(solve(eq[2], c)); `ϵss1` := subs(c = cc, `ϵs`); Ta2 := subs(c = cc, T[2]); cc3 := max(solve(eq[3], c)); `ϵss3` := subs(c = cc, `ϵs`); Ta3 := subs(c = cc, T[3]); M[i] := b*cc^2*fcm*(int(sigma*`ϵc`, `ϵc` = 0 .. `ϵcm`))*10^(-6)/`ϵcm`^2+T*(d-cc)*10^(-6); phi[i] := `ϵcm`/cc; cd[i] := cc/d; print(`ϵcm`, `ϵss`, Ta/As); i := i+1

7501.663386*c

 

308.9753080

 

0.1055633997e-1

 

6333803.979

 

159.9645223

 

0.1055633997e-1

 

1200000

 

Error, complex argument to max/min

 

0.1055633997e-1

 

1254623.652

 

0.35e-2, 0.9349901112e-2, 400

(11)

``

``

``

``

``

NULL

NULL

``

``

``

``

``

``

``

``

``

``

``

``

``

NULL

``

``

``

``

``

``

``

``

``


 

Download mathias.mw

simple non complex solutions ?

Hey mapleprimes

I am trying to get the sum function to give me n variables, well not sure how to explain it, so let me show:

sum(x__i, i = 0 .. 2)

I was expecting to get x__0 + x__1 + x__2

But instead i got:   3*x__i

Thanks for the help in advance :)

Edit:

My problem is gone, but must be some weird character that is playing games with me or something


 

Hi

I m trying to do a nested loop but for some reason maple tells me my for loop is unterminated 

I d love some feedback on my code and on how to make it work
 

NULL

restart

with(Student[Calculus1]):


#geometry [mm]

b := 250:

h := 720:

ds := 70:

d := h-ds:

As := 3000:


#concrete [MPa]

fck := 30:

fcm := fck+8;

38

(1)

Ecm := 33000;

33000

(2)

`ϵc1` := 2.2*(1/1000);

0.2200000000e-2

(3)

eta := `ϵc`/`ϵc1`;

454.5454545*`ϵc`

(4)

Ec1 := fck/`ϵc1`;

13636.36364

(5)

k := 1.05*Ecm/Ec1;

2.540999999

(6)

sigma := fcm*(-eta^2+eta*k)/(1+(k-2)*eta);

38*(-206611.5702*`ϵc`^2+1154.999999*`ϵc`)/(1+245.9090904*`ϵc`)

(7)


#steel [MPa]

Es := 200000:

fsy := 400:

fsu := 600:

`ϵy` := fsy/Es;

1/500

(8)

`ϵsh` := 0.9e-2:

`ϵsu` := 0.75e-1:

P := 4:

`ϵs` := `ϵcm`*(d-c)/c;

`ϵcm`*(650-c)/c

(9)

i := 1;

1

(10)

"for epsiloncm from 0.1/(1000) by 0.1/(1000)to (3.5)/(1000) do  print(epsiloncm);     epsilons:=(epsiloncm)/(c)*(d-c):  T[1]:=epsilons*Es*As:  T[2]:=fsy*As:  T[3]:=fsu+(fsy-fsu)*((epsilonsu-epsilons)/((epsilonsu-epsilonsh)))^(P):  T[4]:=0:  C:=(b*c)/(epsiloncm)*int(sigma,epsilonc=0..epsiloncm):  eq[1]:=T[1]=C:  epsilonl[1]:=epsilony:  eq[2]:=T[2]=C:  epsilonl[2]:=epsilonsh:  eq[3]:=T[3]=C:  a:=1          for j from 1 by 1 to 3 do          if a=1 then          cc:=max(solve(eq[j],c)):          epsilons:=subs(c=cc,epsilons):          if epsilons<=epsilonl[j] then T:=subs(c=cc,T[j]): a:=0:          end if:          end if:          end do:  M[i]:=(b*cc^(2)*fcm)/(epsiloncm^(2))*int(sigma*epsilonc,epsilonc=0..epsiloncm)*10^(-6)+T*(d-cc)*10^(-6):  phi[i]:=(epsiloncm)/(cc):  T[i]:=T:  cd[i]:=(cc)/(d):    print(M[i],epsiloncm,phi[i]);  i:=i+1:  end do:                                        "

Error, unterminated loop

"for epsiloncm from 0.1/1000 by 0.1/1000to 3.5/1000 do  print(epsiloncm);     epsilons:=epsiloncm/c*(d-c):  T[1]:=epsilons*Es*As:  T[2]:=fsy*As:  T[3]:=fsu+(fsy-fsu)*((epsilonsu-epsilons)/(epsilonsu-epsilonsh))^P:  T[4]:=0:  C:=(b*c)/epsiloncm*int(sigma,epsilonc=0..epsiloncm):  eq[1]:=T[1]=C:  epsilonl[1]:=epsilony:  eq[2]:=T[2]=C:  epsilonl[2]:=epsilonsh:  eq[3]:=T[3]=C:  a:=1         for j from 1 by 1 to 3 do          if a=1 then  cc:=max(solve(eq[j],c)):  epsilons:=subs(c=cc,epsilons):          if epsilons<=epsilonl[j] then T:=subs(c=cc,T[j]): a:=0:          end if:          end if:          end do:  M[i]:=(b*cc^2*fcm)/(epsiloncm^2)*int(sigma*epsilonc,epsilonc=0..epsiloncm)*10^(-6)+T*(d-cc)*10^(-6):  phi[i]:=epsiloncm/cc:  T[i]:=T:  cd[i]:=cc/d:  print(M[i],epsiloncm,phi[i]);  i:=i+1:  end do:                                        "

 

``

``

``

``

``

``

NULL

NULL

``

``

``

``

``

``

``

``

``

``

``

``

``

NULL

``

``

``

``

``

``

``

``

``


 

Download mathias.mw

 

Thanks in advance

 

Hi all, 

Please help me to solve this system of equations.


 

Download System.mwSystem.mw

Hi everyone,
I have to make a project that need to use openmaple in java. The problems is when i get the output data of the textCallBack method, i try to convert the output to utf-8 format using the UTF8Encoder provided in java Maple JAR file in order to display Vietnamese language but it doesn't work (the output type is String). Are there any solutions? Thank you.

Hi everybody,

I have two questions about DataSummarize when used with "summarize=embed"

  1. Is it possible to set the width of the table DataSummarize returns ?
    (I thought DataSummarize used DocumentTools[Tabulate] but the option "width=..." doesn't work)
  2. Is it possible to save this table programmatically (for instance in a jpeg file) ?
     

Thanks in advance.

 

 

 

Here are two inocent Vectors.

V1 := Vector[row]([`+`, 1, 2]);

_rtable[18446744074366295326]

(1)

V2 := Vector[row]([1, `+`, 2]);

_rtable[18446744074366288110]

(2)

lprint(V1);

Vector[row](3, {1 = `+`, 2 = 1, 3 = 2}, datatype = anything, storage = rectangular, order = Fortran_order, shape = [])

 

 

The questions:

 

1.  Try to guess the Maple's answer for:

 

V1, V2;

 

2. Try to explain.

 


Download v1v2.mw

Hi All;

I have a maple code made a text file which is readable by another software. My code works as a worksheet but when I run it as a .mpl file does not work well. The text file created in this situation has a limited in the word typing (split a long line into two line by "\" ) . So the text file is not readable by next software (shows a syntax error).

example:


f6 = -2*u*(-u^2+ux-1)-2*ux*(-2*u*ux+uxx)-2*uy*(-2*u*uy+uxy)-uy+10001/10000*u+.438112544881510e-4

 

in 


f6 = -2*u*(-u^2+ux-1)-2*ux*(-2*u*ux+uxx)-2*uy*(-2*u*uy+uxy)-uy+10001/10000*u+.\
438112544881510e-4

I would appreciate for your help.

Baran

HI!, 

first time post.

I am reletivly new to maple and am having some problems getting maple to solve my differential equasions

basically my equasion is -

m*dv/dt=mg-R 

v(0)=0

m is mass at 100kg, g is gravity at 9.81m/s and R is air resistence (k*v) where v is a constant yet to be determined

the constant k needs to be that which would make the terminal velocty 5m/s 

this is all i have so far and im scratching my head:

de := 100*diff(v(t),x)=100*9.81-R;
                      de := 0 = 981.00 - R
Y :=rhs(dsolve({de,v(0)=0}, v(t)));
Error, (in dsolve) not a system with respect to the unknowns [v(t)]
 

what am i doing wrong?. 

Below is a link to my worksheet.  I am attempting to evaluate a Fourier series for a particular number of terms & at a specific location xmax.  As far as I can tell xmax is assessed correctly.  Howver, when I go to evaluate the Fourier series for x = xmax using either the 'value' or 'evalf' command they do on seem to recognize that x = xmax.

So I am guessing I must have some syntax problem.  Can anyone show me what I have wrong?

command_syntax.mw

I'm having sone trouble with my matricies after running restart. When creating the matrice, everything is fine and they are defined as normal. But, if I run restart and try to erdefine them, they are defined as some form of 

The only way I've been able to fix this besides restarting Maple is to either add or remove the colon after the matrix, depending on wether there was one. If ther was, when I remove it the matrix is defined normally and vice versa. I'm relatively new to Maple coming from Matlab and would like to start using it more, but this matrix thing is very frustrating. I've attached a copy of the maple file. Thanks

First definition, everything is normal:

After restart, everything goes bad:

I want to animate a ball rolling on the surface cos(abs(x)+abs(y)).  The ball mass m is 1kg, radius r is 0.1meters starts at (0.5,0.5,cos(abs(5)+abs(5))) meters using g=9.8.  If we say the initial velocity of the ball is pushed in some random direction.  How do I show the path of the ball and animate?

Howto make all points start from [0,0] and goto [1,0] 
and make my loop not overwrite my initial conditions for the next steps?

M:=10; N:=1000;

for i to M do
X[i, 0] := 0;
Y[i, 0] := 0;
X[i, 1] := 1;
Y[i, 1] := 0;
for j to N do
Vinkel := c*rand()-Pi;
X[i, j] := X[i, j-1]+cos(Vinkel);
Y[i, j]:= Y[i, j-1]+sin(Vinkel)
end do:
end do :

 

 

 

Hi...  I use the "solve" command to solve an algebraic  system. Sometimes, the solution is given using some unknowns as parameters.

A toy example: x+y+z=1,x-y+3z=7 gives x=-2z+4,y=z-3 as a function of the parameter z (in the RHS).

Can I obtain directly (without going through the block-by-block solution) the set of parameters used in the given solution?

Thanks in advance.

First 872 873 874 875 876 877 878 Last Page 874 of 2215