Quantcast
Channel: MaplePrimes - Questions and Posts tagged with 2dmath
Viewing all 296 articles
Browse latest View live

Multiplying Matrices--editing initial

$
0
0

Maple will multiply two matrices for me perfectly.  

However, when I change the value in one of the cells and try to re-evaluate the product, Maple starts giving me the sum of something in all of my cells.   

 

Are you simply not allowed to ever change the initial values you set for a matrix?  This doesn't seem to make sense to me.  

I have attached a screen shot to show an example with a random matrix.  


Error "unexpected ', ' in a for do loop (2 loops)

$
0
0

I attach a file in which i get an error though it performs.

I made a for do statement which performs perfectly for one loop, but gives an error "unexpected ',' , though performs (i cld not find the errors!).
What does this error state then? How do i remove when I use semicolon to end statement? Can any one help me. Thanks.
 

Data made for Data[i,j] [k] where i =1,2; j=1..8. Each Data has 6 values. k = 1..6

For each data i got an output list of 6 values using a module;

I made a for do statement which performs but gives an error "unexpected ','.

Can any onr help to find my mistake. Thanks.

 

restart

R1 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 2.36; BTE := .271; sfc := .308; Data[1, 1] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]
NULL

R2 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 2.75; BTE := .298; sfc := .282; Data[1, 2] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R3 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 3.14; BTE := .318; sfc := .264; Data[1, 3] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R4 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 3.53; BTE := .330; sfc := .253; Data[1, 4] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R5 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 3.88; BTE := .334; sfc := .249; Data[1, 5] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R6 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 3.93; BTE := .334; sfc := .249; Data[1, 6] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R7 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 4.32; BTE := .330; sfc := .253; Data[1, 7] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R8 := "SE2"; BP__ref := 3.93; BTE__optimum := .334; BP := 4.71; BTE := .318; sfc := .264; Data[1, 8] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

NULL

R1 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 2.36; BTE := .323; sfc := .258; Data[2, 1] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R2 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 2.75; BTE := .350; sfc := .239; Data[2, 2] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R3 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 3.14; BTE := .369; sfc := .229; Data[2, 3] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R4 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 3.53; BTE := .380; sfc := .219; Data[2, 4] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R5 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 3.88; BTE := .383; sfc := .217; Data[2, 5] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R6 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 3.88; BTE := .383; sfc := .217; Data[2, 6] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R7 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 4.32; BTE := .378; sfc := .221; Data[2, 7] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

R8 := "4S"; BP__ref := 3.88; BTE__optimum := .383; BP := 4.71; BTE := .365; sfc := .229; Data[2, 8] := [R1, BP__ref, BTE__optimum, BP, BTE, sfc]

 

``

Data[1, 8][4]

4.71

(1)

``

``

for i to 2 do for j to 8 do y[i, j] := action1:-Mproc(Data[i, j]); print(Data[i, j], y[i, j]) end do end do

["4S", 3.88, .383, 4.71, .365, .229], ["4S", .4311313490, 0.5411807293e-1, 0.3357509744e-1, .4235561571, 0.6183000001e-1]

(2)

NULL

N1 := NULL; for j to 8 do N1 := N1, Data[1, j][4] end do; print(N1)

2.36, 2.75, 3.14, 3.53, 3.88, 3.93, 4.32, 4.71

(3)

for i to 2 do N[i] := NULL; for j to 8 do N[i] := N[i], Data[i, j][4] end do end do; print(N[2])

Error, `,` unexpected

 

N[2]

2.36, 2.75, 3.14, 3.53, 3.88, 3.88, 4.32, 4.71

(4)

N[1]

2.36, 2.75, 3.14, 3.53, 3.88, 3.93, 4.32, 4.71

(5)

``


 

Download Doubt_A_ReactantValues.mw

Doubt_A_ReactantValues.mw

CustomGroup operations

$
0
0

This fails (the code computes the quotient group directly from the definition):

with(GroupTheory):

Qgr := (sgr, gr) -> CustomGroup(LeftCosets(sgr, gr),
  `.` = ((a, b) -> LeftCoset(Representative(a) . Representative(b), sgr)),
  `/` = ((a) -> LeftCoset(Representative(a)^(-1), sgr)),
  `=` = ((a, b) -> Elements(a) = Elements(b)));

gr := QuaternionGroup();
Generators(gr);
                [(1, 2, 3, 4)(5, 6, 8, 7), (1, 5, 3, 8)(2, 7, 4, 6)]

qgr := Qgr(Subgroup({Generators(gr)[1]^2}, gr), gr);
                  qgr := ` < a custom group with 4 generators > `

AreIsomorphic(qgr, DirectProduct(CyclicGroup(2), CyclicGroup(2)));
Error, (in =) invalid keyword expression

If I rename the parameters in the definition of Qgr, e.g., change the last line to

`=` = ((aa, bb) -> Elements(aa) = Elements(bb))

then everything works, AreIsomorphic gives true.

 

Bug in online Help page

How to display an unevaluated cross product like R x V with the LinearAlgebra package in Maple?

$
0
0

I am using the LinearAlgebra package to do dynamics between a rotating Cartesian coordinate system and a fixed Cartesian coordinate system. The VectorCalculus package is not what I need. 

Since I can't seem to get my test worksheet to paste into this post, I will manually enter an "approximation" to it. I assume that the notation [x, y, z] represents a column vector. I also assume that x represents the cross product operator from the operator pallete. 

I just want to get any one of the three ways of doing a vector cross product (see below) to simply display in math notation as R x V. What I get from the three methods below for an unevaluated cross product is "ugly".

Any help or advice will be greatly appreciated.

> restart

> with(LinearAlgebra):

> R := Vector(3, [x, y, z])

                                                 R := [x, y, z]

> V := Vector(3, [u, v, w])

                                                 V := [u, v, w]

>R x V

                                                [-vz + wy, uz - wx, -uy + vx]

>'R x V'

                                               Typesetting:-delayCrossProduct(R,V)

>CrossProduct(R, V)

                                               [-vz + wy, uz - wx, -uy + vx]

>'CrossProduct(R,V)'

                                               LinearAlgebra:-CrossProduct(R,V)

> R &x V

                                               [-vz + wy, uz - wx, -uy + vx]

'R &x V'

                                               LinearAlgebra:-&x(R,V)

2d input glitch

the elements of matrix was shown in a bad way in maple 2017

$
0
0

in maple 2017 I have enconterd a problem in showing of matrix elements. this error did not exist in maple 2016.

after some mathematical operations with(Linearalgebra)

I have obtained a matrix but maple shows it as below

 

`Vector(3, {(1) = \`+\`(\`*\`(0.19498673e10, \`*\`(\`^\`(x[1], 3))), \`-\`(\`*\`(0.36539123e10, \`*\`(\`^\`(x[1], 2), \`*\`(x[2])))), \`*\`(0.12514628e10, \`*\`(\`^\`(x[1], 2), \`*\`(x[3]))), \`-\`(\`*\`(7987.1506, \`*\`(\`^\`(x[1], 2)))), \`*\`(0.35170874e10, \`*\`(x[1], \`*\`(\`^\`(x[2], 2)))), \`-\`(\`*\`(0.41418339e10, \`*\`(x[1], \`*\`(x[2], \`*\`(x[3]))))), \`-\`(\`*\`(9187.8615, \`*\`(x[1], \`*\`(x[2])))), \`*\`(0.28788960e10, \`*\`(x[1], \`*\`(\`^\`(x[3], 2)))), \`*\`(2871.4325, \`*\`(x[1], \`*\`(x[3]))), \`-\`(\`*\`(0.95594712e9, \`*\`(\`^\`(x[2], 3)))), \`*\`(0.21616680e10, \`*\`(\`^\`(x[2], 2), \`*\`(x[3]))), \`-\`(\`*\`(16552.293, \`*\`(\`^\`(x[2], 2)))), \`-\`(\`*\`(0.20009067e10, \`*\`(x[2], \`*\`(\`^\`(x[3], 2))))), \`-\`(\`*\`(30447.582, \`*\`(x[2], \`*\`(x[3])))), \`*\`(0.42009564e9, \`*\`(\`^\`(x[3], 3))), \`-\`(\`*\`(42372.407, \`*\`(\`^\`(x[3], 2))))), (2) = \`+\`(\`-\`(\`*\`(0.14713991e11, \`*\`(\`^\`(x[1], 3)))), \`*\`(0.42489026e11, \`*\`(\`^\`(x[1], 2), \`*\`(x[2]))), \`-\`(\`*\`(0.25018214e11, \`*\`(\`^\`(x[1], 2), \`*\`(x[3])))), \`-\`(\`*\`(55498.093, \`*\`(\`^\`(x[1], 2)))), \`-\`(\`*\`(0.34645652e11, \`*\`(x[1], \`*\`(\`^\`(x[2], 2))))), \`*\`(0.52229109e11, \`*\`(x[1], \`*\`(x[2], \`*\`(x[3])))), \`-\`(\`*\`(399927.96, \`*\`(x[1], \`*\`(x[2])))), \`-\`(\`*\`(0.24172439e11, \`*\`(x[1], \`*\`(\`^\`(x[3], 2))))), \`-\`(\`*\`(367829.40, \`*\`(x[1], \`*\`(x[3])))), \`*\`(0.11603445e11, \`*\`(\`^\`(x[2], 3))), \`-\`(\`*\`(0.24618238e11, \`*\`(\`^\`(x[2], 2), \`*\`(x[3])))), \`*\`(163169.71, \`*\`(\`^\`(x[2], 2))), \`*\`(0.24814370e11, \`*\`(x[2], \`*\`(\`^\`(x[3], 2)))), \`-\`(\`*\`(711030.81, \`*\`(x[2], \`*\`(x[3])))), \`-\`(\`*\`(0.64758011e10, \`*\`(\`^\`(x[3], 3)))), \`*\`(427083.52, \`*\`(\`^\`(x[3], 2)))), (3) = \`+\`(\`*\`(0.13627334e11, \`*\`(\`^\`(x[1], 3))), \`-\`(\`*\`(0.67651416e11, \`*\`(\`^\`(x[1], 2), \`*\`(x[2])))), \`*\`(0.94045968e11, \`*\`(\`^\`(x[1], 2), \`*\`(x[3]))), \`*\`(46901.078, \`*\`(\`^\`(x[1], 2))), \`*\`(0.70616015e11, \`*\`(x[1], \`*\`(\`^\`(x[2], 2)))), \`-\`(\`*\`(0.13072873e12, \`*\`(x[1], \`*\`(x[2], \`*\`(x[3]))))), \`-\`(\`*\`(994642.52, \`*\`(x[1], \`*\`(x[2])))), \`*\`(0.41170263e11, \`*\`(x[1], \`*\`(\`^\`(x[3], 2)))), \`-\`(\`*\`(2768390.4, \`*\`(x[1], \`*\`(x[3])))), \`-\`(\`*\`(0.22189949e11, \`*\`(\`^\`(x[2], 3)))), \`*\`(0.67100205e11, \`*\`(\`^\`(x[2], 2), \`*\`(x[3]))), \`-\`(\`*\`(961344.42, \`*\`(\`^\`(x[2], 2)))), \`-\`(\`*\`(0.52533380e11, \`*\`(x[2], \`*\`(\`^\`(x[3], 2))))), \`*\`(2309741.6, \`*\`(x[2], \`*\`(x[3]))), \`*\`(0.23373385e11, \`*\`(\`^\`(x[3], 3))), \`-\`(\`*\`(2708629.4, \`*\`(\`^\`(x[3], 2)))))})`

 

 

however its elemets is shown correctly

for example the first element is shown as:

x(1) ^ 3 * 0.194986727121702600e10 + x(1) ^ 2 * x(2) * (-0.365391227687975264e10) + x(1) ^ 2 * x(3) * 0.125146284187383199e10 + x(1) ^ 2 * (-0.798715058230362229e4) + x(1) * x(2) ^ 2 * 0.351708738191481209e10 + x(1) * x(2) * x(3) * (-0.414183393126536083e10) + x(1) * x(2) * (-0.918786154503699254e4) + x(1) * x(3) ^ 2 * 0.287889598946906424e10 + x(1) * x(3) * 0.287143252997207219e4 + x(2) ^ 3 * (-0.955947121705965400e9) + x(2) ^ 2 * x(3) * 0.216166802966718006e10 + x(2) ^ 2 * (-0.165522927641286878e5) + x(2) * x(3) ^ 2 * 

2d input glitch


Parentheses around y'(t) in Maple?

$
0
0

Is there a way to remove the parentheses around y'(t)?

Thank you!

mapleatha
 

 

How to not execute part of the calculation for the display

$
0
0

In the context of homework, there are some calculations that I do through Maple, and others that I do mostly by hand, and just show my work with Maple's engine that is great at rendering special mathematical formats (like matrices, for example).

I would like to show something like that:

 

It's pretty easy to make the matrices, right now my line looks like this:

 

But I would like to refrain for executing the "+" operator and the multiplication by x6 and x7, to have the output like in the first picture. If I disable all execution and take the statement as plain text, then obviously I will not get the correctly rendered matrices.

Is there a way, like an escape character, to prevent specific operations from being executed in one line in particular?

ctrl+= doesnt display result inline

$
0
0

I'm learning Maple with the provided user manual, and it clearly says to use ctrl+= to display inline results. It simply doesn't work on  my system. Running Maple 2017.3. Is this a bug?

 

HELP !! Maple doesn t want to give me the integral expression

$
0
0


Hi

I have to create a function "I(epsilonc)" by integrating a function "sigma(y)"  from 0 to the variable "epsilonc"

but maple doesn t give me the actual formula of the integrated function, it only gives me equation 4 which is utterly useless to me.

 

How do I get Maple to actually calculate and give me the integrated function ??

Any help will be greatly appreciated

restart


#data

`&varepsilon;c1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`&varepsilon;cu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b = 250:

d := 450:


#coefficients

eta := epsilon/`&varepsilon;c1`:

k := 1.05*Ecm*`&varepsilon;c1`/fcm:

``


#formula

 

``

`&sigma;c` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc

(2)

plot(`&sigma;c`, 0 .. `&varepsilon;cu1`);

 

``

#Pressure as a function of y (y=`&varepsilon;c`*y/x) :

`&sigma;c` := proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc

(3)

Iy := proc (y) options operator, arrow; int(fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)), y = 0 .. `&varepsilon;c`/(`&varepsilon;c`+Fy/Es)) end proc;

proc (y) options operator, arrow; int(fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)), y = 0 .. `&varepsilon;c`/(`&varepsilon;c`+Fy/Es)) end proc

(4)

Iy(`&varepsilon;c`/(`&varepsilon;c`+Fy/Es))

Error, (in int) integration range or variable must be specified in the second argument, got `&varepsilon;c`/(`&varepsilon;c`+1/400) = 0 .. `&varepsilon;c`/(`&varepsilon;c`+1/400)

 

solve(Iyy*b = Fy*As, `&varepsilon;c`);

``

``


 

Download HW1_-_EC2_strain-pressure_graph.mw
 

restart


#data

`&varepsilon;c1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`&varepsilon;cu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b = 250:

d := 450:


#coefficients

eta := epsilon/`&varepsilon;c1`:

k := 1.05*Ecm*`&varepsilon;c1`/fcm:

``


#formula

 

``

`&sigma;c` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc

(2)

plot(`&sigma;c`, 0 .. `&varepsilon;cu1`);

 

``

#Pressure as a function of y (y=`&varepsilon;c`*y/x) :

`&sigma;c` := proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc

(3)

Iy := proc (y) options operator, arrow; int(fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)), y = 0 .. `&varepsilon;c`/(`&varepsilon;c`+Fy/Es)) end proc;

proc (y) options operator, arrow; int(fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)), y = 0 .. `&varepsilon;c`/(`&varepsilon;c`+Fy/Es)) end proc

(4)

Iy(`&varepsilon;c`/(`&varepsilon;c`+Fy/Es))

Error, (in int) integration range or variable must be specified in the second argument, got `&varepsilon;c`/(`&varepsilon;c`+1/400) = 0 .. `&varepsilon;c`/(`&varepsilon;c`+1/400)

 

solve(Iyy*b = Fy*As, `&varepsilon;c`);

``

``


 

Download HW1_-_EC2_strain-pressure_graph.mw

 

Help ...Maple will not plot

$
0
0

Hi

II ve managed to build an expression with one variable 
I me trying to plot this expression on a defined range but maple doesn t not let me.

How do I manage to plot this expression ?

If you have any advise on how to improve my code I m open to comments 

 

Thanks a lot in advance
 

restart


#data

`&varepsilon;c1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`&varepsilon;cu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b := 250:

d := 450:


#coefficients

eta := epsilon/`&varepsilon;c1`:

k := 1.05*Ecm*`&varepsilon;c1`/fcm:

NULL


#formula

 

NULL

`&sigma;c` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc

(2)

plot(`&sigma;c`, 0 .. `&varepsilon;cu1`);

 

NULL

#Pressure as a function of y (y=`&varepsilon;c`*y/x) :

`&sigma;c` := proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc

(3)

`assuming`([int(`&sigma;c`, 0 .. x)], [0 <= x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

int(25*(1302.000000*`&varepsilon;c`*x/x-226757.3696*`&varepsilon;c`^2*x^2/x^2)/(1+349.6190476*`&varepsilon;c`*x/x), x = 0 .. x)

(4)

C := `assuming`([int(fcm*(k*`&varepsilon;c`*s/(x*`&varepsilon;c1`)-(`&varepsilon;c`*s/(x*`&varepsilon;c1`))^2)*b/(1+(k-2)*`&varepsilon;c`*s/(x*`&varepsilon;c1`)), s = 0 .. x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

-0.7487980799e-19*x*(0.2706771684e26*`&varepsilon;c`^2+0.1331955800e22*ln(874047619.*`&varepsilon;c`+2500000.)-0.1962210817e23-0.4656771182e24*`&varepsilon;c`)/`&varepsilon;c`

(5)

T := `assuming`([(d-x)*`&varepsilon;c`*Es*As/x], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

377000000*(450-x)*`&varepsilon;c`/x

(6)

``

x := `assuming`([solve(C = T, x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

868.0577815*(0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.), -868.0577815*(-0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.)

(7)

``

``

``

l

(8)

plot(l, `&varepsilon;c` = 0 .. 0.34e-2)

Error, (in plot) expected a range but received `&varepsilon;c` = 0 .. 0.34e-2

 

``

NULL


 

Download HW1_-_EC2_strain-pressure_graph.mw
 

restart


#data

`&varepsilon;c1` := 2.1*10^(-3);

0.2100000000e-2

(1)

`&varepsilon;cu1` := 3.5*10^(-3):

Ecm := 31000:

fcm := 25:

Fy := 500:

Es := 200000:

#geometry

As := 1885:

b := 250:

d := 450:


#coefficients

eta := epsilon/`&varepsilon;c1`:

k := 1.05*Ecm*`&varepsilon;c1`/fcm:

NULL


#formula

 

NULL

`&sigma;c` := proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc;

proc (epsilon) options operator, arrow; fcm*(k*epsilon/`&varepsilon;c1`-epsilon^2/`&varepsilon;c1`^2)/(1+(k-2)*epsilon/`&varepsilon;c1`) end proc

(2)

plot(`&sigma;c`, 0 .. `&varepsilon;cu1`);

 

NULL

#Pressure as a function of y (y=`&varepsilon;c`*y/x) :

`&sigma;c` := proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc;

proc (y) options operator, arrow; fcm*(k*`&varepsilon;c`*y/(x*`&varepsilon;c1`)-`&varepsilon;c`^2*y^2/(x^2*`&varepsilon;c1`^2))/(1+(k-2)*`&varepsilon;c`*y/(x*`&varepsilon;c1`)) end proc

(3)

`assuming`([int(`&sigma;c`, 0 .. x)], [0 <= x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

int(25*(1302.000000*`&varepsilon;c`*x/x-226757.3696*`&varepsilon;c`^2*x^2/x^2)/(1+349.6190476*`&varepsilon;c`*x/x), x = 0 .. x)

(4)

C := `assuming`([int(fcm*(k*`&varepsilon;c`*s/(x*`&varepsilon;c1`)-(`&varepsilon;c`*s/(x*`&varepsilon;c1`))^2)*b/(1+(k-2)*`&varepsilon;c`*s/(x*`&varepsilon;c1`)), s = 0 .. x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

-0.7487980799e-19*x*(0.2706771684e26*`&varepsilon;c`^2+0.1331955800e22*ln(874047619.*`&varepsilon;c`+2500000.)-0.1962210817e23-0.4656771182e24*`&varepsilon;c`)/`&varepsilon;c`

(5)

T := `assuming`([(d-x)*`&varepsilon;c`*Es*As/x], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

377000000*(450-x)*`&varepsilon;c`/x

(6)

``

x := `assuming`([solve(C = T, x)], [0 < x and x <= (1/2)*d, 0 <= `&varepsilon;c` and `&varepsilon;c` <= 0.35e-2]);

868.0577815*(0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.), -868.0577815*(-0.2900000000e12*`&varepsilon;c`+(-0.7297483689e24*`&varepsilon;c`^2+0.1400157114e23*`&varepsilon;c`-0.4004807874e20*ln(874047619.*`&varepsilon;c`+2500000.)+0.5899803379e21)^(1/2))*`&varepsilon;c`/(0.2706771684e13*`&varepsilon;c`^2+133195580.*ln(874047619.*`&varepsilon;c`+2500000.)-0.4656771182e11*`&varepsilon;c`-1962210817.)

(7)

``

``

``

l

(8)

plot(l, `&varepsilon;c` = 0 .. 0.34e-2)

Error, (in plot) expected a range but received `&varepsilon;c` = 0 .. 0.34e-2

 

``

NULL


 

Download HW1_-_EC2_strain-pressure_graph.mw

 

 

Wrong output when solving an equation

$
0
0

Hi.

 

When I try to solve the equation 15-1/100M=5+1/600M. The result in maple is 7/6000. But the “real” result according to other calculating programs is 6000/7. Is there a setting that is causing this to happen?

A problem with matrix vector multiplication

$
0
0

Hi everibody 

I work with Maple 2015 under OS-X El Capitan.

Using more than one matrix vector product (either M.V  or MatrixVectorMultiply(M,V)  ; M is a n by p matrix and V a column vector of size p) within the same block of commands generates an error.

Do other people have the same problem ?
Thanks for your feedback.

SomethingGoesWrong.mw


PS : I know I can do this   X . <<1, 1, -1> | <-1, 2, 0>> but this doesn't explain the error I get

 


Error, (in fsolve) expecting an equation or set or list of equations or list of procedures, but received a combination

$
0
0

Hello everyone,

Compliment of the season. I received the error message while trying to solve some Pade approximant problem and try to find what could possibly be the error but could not.

Please, anyone with good informations on the possible solutions should help.

Thanks. The Link is below

 Error.mw

Typesetting:-delayDotProduct

$
0
0

hi every one

 im new in maple, i have this problem , the output result (nod) dont disply after the for loop how can i disply it ?( sorry for this question, im new in this prog.)

the 2nd question is the (Typesetting:-delayDotProduct) how can i resolve this???

plz help 

 

for i from 2 to 3 do for j from 2 to 3 do Typesetting[delayDotProduct](rod, [i, j], true) := s*imp*l*i*f*y*subs*({cx^2 = 1-cy^2-cz^2}, Typesetting[delayDotProduct](rod, [i, j], true)) end do end do;
Typesetting[delayDotProduct](rod, [1, 1], true) := collect(rod[1, 1], cos(theta));

Typesetting:-delayDotProduct

Can we overide Maple default dot deri...

$
0
0

Can we overide Maple default dot derivative with 'tau' instead of 't'?

In piecewise notation (real math symbol), how I could add more cases?

$
0
0

In piecewise notation (real math symbol), how I could add more cases?

print text and Equations

$
0
0

How do you way text below

Maple 2015 

if i use function cat() is not 

 

Viewing all 296 articles
Browse latest View live