MaplePrimes Questions

To determine whether or not equilibrium points of a nonlinear ordinary differential equation is globally stable, a Lyapunov function is often employed. Since there are no general methods for constructing Lyapunov functions, May I know if there are methods in MAPLE that can be used to determine these Lyapunov functions? Thanks

I have a function inside a file called proc2.mpl

When I read the file from a worksheet using the command read("proc2.mpl"); I see warning messages

Warning, incomplete string;  use " to end the string
Warning, incomplete string;  use " to end the string

printed on the screen.  Now if I copy the same exact code and past it into the worksheet and invoke the function, I see no warning messages.

Why is that? Here is the content of "proc2.mpl"

process_file := proc()
  local str,fileName;


  fileName := "output.txt";
  
  str:="
   \\begin{align*}
     A =& B  \\\\ 
       =& 3
   \\end{align*}
  ";
	
  writebytes(fileName, str);
  close(fileName);

end proc:

Here is screen shot loading the function from worksheet

But now when I call it using process_file(); it works as expected.

Now here is screen shot with same code inside the worksheet itself. No warning messages show up when I call it.

And the above call also works as expected.

So I am now ignoring these warning messages since I also do not know what causes them.

Any ideas why they show up when the code is inside a .mpl file and not in the work sheet?

Maple 2017.3 on windows.

 

 

 

Dear Community,

 

I tried to run a very simple dumped oscillator in MapleSim. Upon running however I do not get any result plots. Nothing appears in the result window :-/ Could you pls. have a short look at it?

Tx in advance,

best regards

Andras

RLC_Passive_Low_Pass_1.msim

Just want some input if anyone thinks this is a bug or not

Hiding the contents of equation labels in one table (table -> properties -> uncheck show equation labels) removes all reference to the labels within that table.  Is that supposed to occur?

The table below has show equation labels checked.  If I uncheck the show equation labels in the first table I would expect the reference labels (1) and (2) to disappear and (3) and (4) references in the next table to remain unchanged.

However unchecking show equation labels in the first table relabels the two equations in the second table to (1) and (2) as shown below.  Is this a bug?

However this doesn't disrupt further content in the worksheet if references were made to equation label (1). After unchecking show equation labels in the first table, all original references to label (1) are replaced with the actual value (sin(x))

Could someone help me with the following.  The syntax produces an unfinished graph with a warning.

> with(plots);
> z := polar(1.05, (1/10)*Pi);
                             
> display(polarplot(1, color = grey, axis[radial] = [color = "Blue"]), complexplot(seq(evalc(z))^n,
n = 1 .. 21));

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

the graphic looks like this

The graphic should look like this

Thanks, any help appreciated.

Les

 

Pl. help me remove the equation labels or hide the equation labels in my doc enclosed. What should i do to prevent equation labels from forming in my document. Please refer last two rows in the table below from the document.

Thanks for helping.Ramakrishnan V
 

restart

Conversion from ppm to % vol is 1ppm = 0.0001 %

Conversion from % vol to ppm is % = 10000*ppm

Avagadro's law between mass and volume is

"22.4 nm^(3)=1 kg mol; 22.4 lit = 1 g mol; at  101.325 kPa*(normal pressure) and (273+C ) K;"

ppm = 0.8205e-1*T*mg/(M*nm^3); mg/nm^3 = M*ppm/(0.8205e-1*T)

where T is abs temperature and M is molecular weight;

M/(0.08205T)]*ppm*((A/F)*EA+1) /ρex; Emission is mg/kg fuel

If SFC is x g/kWhr; then Mass flow of exhaust is  (x/1000)*[M/(0.08205T)]*ppm*((A/F)*EA+1) /ρex; Emission is mg/kWhr

If the distance covered per liter is DCL km/liter, then (SFC/1000)*(DCL)/`ρ__f`  is distance covered in km/kWhrNULL

[kg*km/(kWh*liter)*(liter/kg) = km/kWh]

1 g/cc = 1 kg/lit; 1kg/m3 = 1 g/lit; Fuel density`ρ__f`  is 0.77 kg/litNULL

Emission in mg/km = 0.77*Emission is mg/kg fuel/(km/liter fuel)

[(kg/liter*(mg/kg))*km/liter = mg/km]

SFC is in g/kWhr  

Measured emission is  E in ppm  ``

Excess Air percentage is EA in per unit point

Molecular Weight of pollutant is M g/mol

Molecular weight of air = .21*32+.79*28 = 28.84 g/mol

 

SFC

       SFC := 200 g/kWh

`ρ__f`

rho__f

(1)

        `ρ__f` := .77   kg/liter

DCL

   DCL := 50    km/liter fuel

A/F

AFR := 14.5``

EA

EA := .1NULL

T__exh

T__exh

(2)

   T__exh := 450     K

`ρ__ex`

rho__ex

(3)

      `ρ__ex` := .457 kg/m3   

 

c := (AFR*EA+1)/`ρ__ex` = 5.361050328 

NULL

CO

CO2

HC

NO2

O2

N2

Air

M [g/mol]

M__co := 28

NULL

M__co2 := 44

NULL

M__hc := 17

M__no2 := 46

M__o2 := 32

M__n2 := 28

M__air := 28.84

Emission [% ] or[ppm]

E__co := 10

                %

E__co2 := 10

                %

E__hc := 110

                ppm

E__no2 := 10

                ppm

E__o2 := 10

              %

E__n2 := 10

                 %

E__air := 10

                    %

Emission [mg/Nm3]

0.8205e-1*T__exh/M__co

1.318660714

(4)

0.8205e-1*T__exh/M__co2

.8391477273

(5)

0.8205e-1*T__exh/M__hc

2.171911765

(6)

0.8205e-1*T__exh/M__no2

.8026630435

(7)

0.8205e-1*T__exh/M__o2

1.153828125

(8)

0.8205e-1*T__exh/M__n2

1.318660714

(9)

0.8205e-1*T__exh/M__air

1.280253121

(10)

Emission

mg/kWh

k__1 := E__co*c/(1.318660714)
= 40.65526690NULL

k__2 := E__co*c/(.8391477273) = 63.88684799NULL

k__3 := E__co*c/(2.171911765)
 = 24.68355489NULL

k__4 := E__co*c/(.8026630435) = 66.79079562NULL

k__5 := E__co*c/(1.153828125)
 = 46.46316216NULL

k__6 := E__co*c/(1.318660714) = 40.65526690 

NULL

k__7 := E__co*c/(1.280253121) = 41.87492489NULL

Emission

mg/kg fuel

k__1f := k__1*SFC
 = 8131.053380NULL

k__2f := k__2*SFC
 = 12777.36960  NULL

k__3f := k__3*SFC
 =

4936.710978

(11)

 

k__4f := k__4*SFC =

13358.15912

(12)

 

k__5f := k__5*SFC
 =

9292.632432

(13)

 

k__6f := k__6*SFC
 = 8131.053380 =

``

k__7f := k__7*SFC
 = 8374.984978 

Emission

[mg/km]

 

Why in the above row, col2 and col3 have no equn nos?

What should i do to remove the equn numbers 11, 12 and 13?

 

 

Temperature of exhaust is in T K

Density of exhaust is say `ρ__ex` kg/m3 

DCL is distance covered per kWh

SFC := 100*g/kWh; ppm := 4.3; EA := .2; M := 28; `ρ__ex` := .457; T = 400

100*g/kWh

(14)

Emission in g/kWhr is

Emission g/g of fuel = Emission in g/kWhr/SFC in g/kWhrper

Emission in mg/km = 0.77*Emission is mg/kg fuel*(km/liter fuel)


 

Download ppm_to_g_per_km_conversion.mw

hi . i want to solve a equation with newton method but i cant, 

restart;
f(x):=x^2-3;
df(x):=diff(f(x),x);
x:=1;
for i from 2 to 5 by 1 do
    x:=x-(f(x)/df(x));
    print(x);
    end do

I want to save multiline string to a file. I am not sure Maple likes my string. Also when I look at the textfile after writing the string to it, the string does not look like what is expected. I given small example

First, created a "proc.mpl" in the folder. Here is the content of "proc.mpl"

process_file := proc()
local str;
  
str:="
   \\begin{align*}
     A =& B \\\\
       =& 3
   \\end{align*}
";

FileTools[Text][WriteString]("output.txt", str);
FileTools[Text][Close]("output.txt");

end proc:

Then opened Maple, and from a worksheet, called the above proc() like this

restart;
currentdir("C:\\where_the_file_is");
read("proc.mpl");
Warning, incomplete string;  use " to end the string
Warning, incomplete string;  use " to end the string
Warning, incomplete string;  use " to end the string
process_file();

You see these warnings there. So how does one make a multiline string in Maple? Please tell me it is possible. Else this whole process will not work.

When I look at the "output.txt", where the string is written, this is what shows


   \begin{align*}
     A =& B \      =& 3
   \end{align*}

Which is wrong, I expected this


   \begin{align*}
     A =& B \\     
        =& 3
   \end{align*}

 

I use \ to escape each \, since this is Latex code. 

QUestions are:

1) Can one have multiline string in Maple? as in

"
line one
line two
"

With implied carriage return, just as it appears in the input?  I do this all the time in Mathematica.

2) Is something wrong with how I am saving the string to the file?
 

On a side note: I noticed I had to close the textfile to see the output in it. It looks like FileTools[Text][WriteString] does not flush the string autmatically to a otuput file after each write() which is a little annoying.

Using Maple 2017.3 on windows.

This is how the above is done in Mathematica

And this is "test.txt" which shows the output as expected

   \begin{align*}
     A =& B \\
       =& 3
   \end{align*}

I'd like to do the same in Maple.

 

Hi, 

How do I get underlined text -- axis labels and legend -- in plot command? Tried many ways, but could not succeed. I would appreciate any inputs.

Regards,

Omkar

 

Hi, There is a problem in solving ODE using dsolve/numeric code.

 

Ba.mw

 

Hi, I do not understand how to solve errors in Maple on my project. My project is to solve Vehicle Routing Problem with Time Windows, and then the error is "Error, (in Optimization: -LPSolve) no feasible integer point found; use feasibilitytolerance option to adjust tolerance". I do not understand about feasibilitytolerance. Can anyone help me? Thankyou.

NULL

HASIL*MAPLE*UNTUK*KECAMATAN*COBLONGNULL

restart

with(Optimization);

[ImportMPS, Interactive, LPSolve, LSSolve, Maximize, Minimize, NLPSolve, QPSolve]

(1)

with(linalg);

[BlockDiagonal, GramSchmidt, JordanBlock, LUdecomp, QRdecomp, Wronskian, addcol, addrow, adj, adjoint, angle, augment, backsub, band, basis, bezout, blockmatrix, charmat, charpoly, cholesky, col, coldim, colspace, colspan, companion, concat, cond, copyinto, crossprod, curl, definite, delcols, delrows, det, diag, diverge, dotprod, eigenvals, eigenvalues, eigenvectors, eigenvects, entermatrix, equal, exponential, extend, ffgausselim, fibonacci, forwardsub, frobenius, gausselim, gaussjord, geneqns, genmatrix, grad, hadamard, hermite, hessian, hilbert, htranspose, ihermite, indexfunc, innerprod, intbasis, inverse, ismith, issimilar, iszero, jacobian, jordan, kernel, laplacian, leastsqrs, linsolve, matadd, matrix, minor, minpoly, mulcol, mulrow, multiply, norm, normalize, nullspace, orthog, permanent, pivot, potential, randmatrix, randvector, rank, ratform, row, rowdim, rowspace, rowspan, rref, scalarmul, singularvals, smith, stackmatrix, submatrix, subvector, sumbasis, swapcol, swaprow, sylvester, toeplitz, trace, transpose, vandermonde, vecpotent, vectdim, vector, wronskian]

(2)

with(ExcelTools);

[Export, Import, WorkbookData]

(3)

with(CodeTools);

[CPUTime, DecodeName, EncodeName, Profiling, RealTime, Test, Usage]

(4)

``

c := convert(Import("C:\\Users\\VaniaMR\\Documents\\SEMANGAT SKRIPSI\\skripsi\\Bab-bab\\data Bandung Utara.xlsx", 6, "B2:I9"), matrix)

array( 1 .. 8, 1 .. 8, [( 5, 8 ) = (2.9), ( 4, 1 ) = (2.5), ( 2, 2 ) = (0.), ( 8, 3 ) = (3.1), ( 2, 4 ) = (1.9), ( 7, 5 ) = (3.7), ( 6, 6 ) = (0.), ( 3, 7 ) = (3.9), ( 6, 8 ) = (2.7), ( 5, 1 ) = (.28), ( 7, 3 ) = (3.9), ( 1, 2 ) = (3.2), ( 3, 4 ) = (3.7), ( 8, 5 ) = (2.9), ( 5, 6 ) = (.26), ( 2, 7 ) = (1.8), ( 3, 8 ) = (3.1), ( 6, 1 ) = (.25), ( 8, 2 ) = (1.1), ( 2, 3 ) = (2.8), ( 4, 4 ) = (0.), ( 5, 5 ) = (0.), ( 8, 6 ) = (2.7), ( 1, 7 ) = (3.5), ( 4, 8 ) = (2.3), ( 7, 1 ) = (3.5), ( 7, 2 ) = (1.8), ( 1, 3 ) = (4.1), ( 5, 4 ) = (2.7), ( 6, 5 ) = (.26), ( 7, 6 ) = (3.4), ( 8, 7 ) = (1.6), ( 1, 8 ) = (3.3), ( 8, 1 ) = (3.3), ( 6, 2 ) = (2.9), ( 4, 3 ) = (3.7), ( 6, 4 ) = (2.7), ( 3, 5 ) = (4.0), ( 2, 6 ) = (2.9), ( 7, 7 ) = (0.), ( 2, 8 ) = (1.1), ( 5, 2 ) = (3.1), ( 2, 1 ) = (3.2), ( 3, 3 ) = (0.), ( 7, 4 ) = (1.3), ( 4, 5 ) = (2.7), ( 1, 6 ) = (.25), ( 6, 7 ) = (3.4), ( 7, 8 ) = (1.6), ( 4, 2 ) = (1.9), ( 6, 3 ) = (3.8), ( 8, 4 ) = (2.3), ( 1, 1 ) = (0.), ( 1, 5 ) = (.28), ( 4, 6 ) = (2.7), ( 5, 7 ) = (3.7), ( 8, 8 ) = (0.), ( 3, 1 ) = (4.1), ( 3, 2 ) = (2.8), ( 5, 3 ) = (4.0), ( 1, 4 ) = (2.5), ( 2, 5 ) = (3.1), ( 3, 6 ) = (3.8), ( 4, 7 ) = (1.3)  ] )

(5)

t := convert(Import("C:\\Users\\VaniaMR\\Documents\\SEMANGAT SKRIPSI\\skripsi\\Bab-bab\\data Bandung Utara1.xlsx", 7, "B2:I9"), matrix)

array( 1 .. 8, 1 .. 8, [( 5, 8 ) = (9.0), ( 4, 1 ) = (8.0), ( 2, 2 ) = (0.), ( 8, 3 ) = (8.0), ( 2, 4 ) = (6.0), ( 7, 5 ) = (10.0), ( 6, 6 ) = (0.), ( 3, 7 ) = (9.0), ( 6, 8 ) = (7.0), ( 5, 1 ) = (2.0), ( 7, 3 ) = (9.0), ( 1, 2 ) = (9.0), ( 3, 4 ) = (12.0), ( 8, 5 ) = (9.0), ( 5, 6 ) = (2.0), ( 2, 7 ) = (4.0), ( 3, 8 ) = (8.0), ( 6, 1 ) = (2.0), ( 8, 2 ) = (3.0), ( 2, 3 ) = (7.0), ( 4, 4 ) = (0.), ( 5, 5 ) = (0.), ( 8, 6 ) = (7.0), ( 1, 7 ) = (10.0), ( 4, 8 ) = (6.0), ( 7, 1 ) = (10.0), ( 7, 2 ) = (4.0), ( 1, 3 ) = (11.0), ( 5, 4 ) = (9.0), ( 6, 5 ) = (2.0), ( 7, 6 ) = (8.0), ( 8, 7 ) = (3.0), ( 1, 8 ) = (9.0), ( 8, 1 ) = (9.0), ( 6, 2 ) = (8.0), ( 4, 3 ) = (12.0), ( 6, 4 ) = (9.0), ( 3, 5 ) = (12.0), ( 2, 6 ) = (8.0), ( 7, 7 ) = (0.), ( 2, 8 ) = (3.0), ( 5, 2 ) = (10.0), ( 2, 1 ) = (9.0), ( 3, 3 ) = (0.), ( 7, 4 ) = (4.0), ( 4, 5 ) = (9.0), ( 1, 6 ) = (2.0), ( 6, 7 ) = (8.0), ( 7, 8 ) = (3.0), ( 4, 2 ) = (6.0), ( 6, 3 ) = (10.0), ( 8, 4 ) = (6.0), ( 1, 1 ) = (0.), ( 1, 5 ) = (2.0), ( 4, 6 ) = (9.0), ( 5, 7 ) = (10.0), ( 8, 8 ) = (0.), ( 3, 1 ) = (11.0), ( 3, 2 ) = (7.0), ( 5, 3 ) = (12.0), ( 1, 4 ) = (8.0), ( 2, 5 ) = (10.0), ( 3, 6 ) = (10.0), ( 4, 7 ) = (4.0)  ] )

(6)

a := `<,>`(0, 0, 0, 0, 0, 0, 0, 0)

a := Vector(8, {(1) = 0, (2) = 0, (3) = 0, (4) = 0, (5) = 0, (6) = 0, (7) = 0, (8) = 0})

(7)

b := `<,>`(30, 30, 30, 30, 30, 30, 30, 30)

b := Vector(8, {(1) = 30, (2) = 30, (3) = 30, (4) = 30, (5) = 30, (6) = 30, (7) = 30, (8) = 30})

(8)

n := sqrt(numelems(c)):

{1, 2, 3, 4, 5, 6, 7, 8}

(9)

z := add(add(c[i, j]*x[i, j], j = N), i = N);

3.2*x[1, 2]+4.1*x[1, 3]+2.5*x[1, 4]+.28*x[1, 5]+.25*x[1, 6]+3.5*x[1, 7]+3.3*x[1, 8]+3.2*x[2, 1]+2.8*x[2, 3]+1.9*x[2, 4]+3.1*x[2, 5]+2.9*x[2, 6]+1.8*x[2, 7]+1.1*x[2, 8]+4.1*x[3, 1]+2.8*x[3, 2]+3.7*x[3, 4]+4.0*x[3, 5]+3.8*x[3, 6]+3.9*x[3, 7]+3.1*x[3, 8]+2.5*x[4, 1]+1.9*x[4, 2]+3.7*x[4, 3]+2.7*x[4, 5]+2.7*x[4, 6]+1.3*x[4, 7]+2.3*x[4, 8]+.28*x[5, 1]+3.1*x[5, 2]+4.0*x[5, 3]+2.7*x[5, 4]+.26*x[5, 6]+3.7*x[5, 7]+2.9*x[5, 8]+.25*x[6, 1]+2.9*x[6, 2]+3.8*x[6, 3]+2.7*x[6, 4]+.26*x[6, 5]+3.4*x[6, 7]+2.7*x[6, 8]+3.5*x[7, 1]+1.8*x[7, 2]+3.9*x[7, 3]+1.3*x[7, 4]+3.7*x[7, 5]+3.4*x[7, 6]+1.6*x[7, 8]+3.3*x[8, 1]+1.1*x[8, 2]+3.1*x[8, 3]+2.3*x[8, 4]+2.9*x[8, 5]+2.7*x[8, 6]+1.6*x[8, 7]

(10)

conx := seq(add(x[i, j], i = `minus`(N, {j})) = 1, j = N);

x[2, 1]+x[3, 1]+x[4, 1]+x[5, 1]+x[6, 1]+x[7, 1]+x[8, 1] = 1, x[1, 2]+x[3, 2]+x[4, 2]+x[5, 2]+x[6, 2]+x[7, 2]+x[8, 2] = 1, x[1, 3]+x[2, 3]+x[4, 3]+x[5, 3]+x[6, 3]+x[7, 3]+x[8, 3] = 1, x[1, 4]+x[2, 4]+x[3, 4]+x[5, 4]+x[6, 4]+x[7, 4]+x[8, 4] = 1, x[1, 5]+x[2, 5]+x[3, 5]+x[4, 5]+x[6, 5]+x[7, 5]+x[8, 5] = 1, x[1, 6]+x[2, 6]+x[3, 6]+x[4, 6]+x[5, 6]+x[7, 6]+x[8, 6] = 1, x[1, 7]+x[2, 7]+x[3, 7]+x[4, 7]+x[5, 7]+x[6, 7]+x[8, 7] = 1, x[1, 8]+x[2, 8]+x[3, 8]+x[4, 8]+x[5, 8]+x[6, 8]+x[7, 8] = 1

(11)

conV := seq(add(x[i, k], i = N)-add(x[k, j], j = N) = 0, k = N);

x[2, 1]+x[3, 1]+x[4, 1]+x[5, 1]+x[6, 1]+x[7, 1]+x[8, 1]-x[1, 2]-x[1, 3]-x[1, 4]-x[1, 5]-x[1, 6]-x[1, 7]-x[1, 8] = 0, x[1, 2]+x[3, 2]+x[4, 2]+x[5, 2]+x[6, 2]+x[7, 2]+x[8, 2]-x[2, 1]-x[2, 3]-x[2, 4]-x[2, 5]-x[2, 6]-x[2, 7]-x[2, 8] = 0, x[1, 3]+x[2, 3]+x[4, 3]+x[5, 3]+x[6, 3]+x[7, 3]+x[8, 3]-x[3, 1]-x[3, 2]-x[3, 4]-x[3, 5]-x[3, 6]-x[3, 7]-x[3, 8] = 0, x[1, 4]+x[2, 4]+x[3, 4]+x[5, 4]+x[6, 4]+x[7, 4]+x[8, 4]-x[4, 1]-x[4, 2]-x[4, 3]-x[4, 5]-x[4, 6]-x[4, 7]-x[4, 8] = 0, x[1, 5]+x[2, 5]+x[3, 5]+x[4, 5]+x[6, 5]+x[7, 5]+x[8, 5]-x[5, 1]-x[5, 2]-x[5, 3]-x[5, 4]-x[5, 6]-x[5, 7]-x[5, 8] = 0, x[1, 6]+x[2, 6]+x[3, 6]+x[4, 6]+x[5, 6]+x[7, 6]+x[8, 6]-x[6, 1]-x[6, 2]-x[6, 3]-x[6, 4]-x[6, 5]-x[6, 7]-x[6, 8] = 0, x[1, 7]+x[2, 7]+x[3, 7]+x[4, 7]+x[5, 7]+x[6, 7]+x[8, 7]-x[7, 1]-x[7, 2]-x[7, 3]-x[7, 4]-x[7, 5]-x[7, 6]-x[7, 8] = 0, x[1, 8]+x[2, 8]+x[3, 8]+x[4, 8]+x[5, 8]+x[6, 8]+x[7, 8]-x[8, 1]-x[8, 2]-x[8, 3]-x[8, 4]-x[8, 5]-x[8, 6]-x[8, 7] = 0

(12)

conz := add(x[i, 1], i = N) = 1;

x[1, 1]+x[2, 1]+x[3, 1]+x[4, 1]+x[5, 1]+x[6, 1]+x[7, 1]+x[8, 1] = 1

(13)

conD := add(x[1, j], j = N) = 1;

x[1, 1]+x[1, 2]+x[1, 3]+x[1, 4]+x[1, 5]+x[1, 6]+x[1, 7]+x[1, 8] = 1

(14)

conTW := seq(seq(y[i]-y[j]+max(b[i]+t[i, j]-a[j], 0)*x[i, j] <= b[i]-a[j], i = `minus`(N, {j})), j = N);

y[2]-y[1]+39.0*x[2, 1] <= 30, y[3]-y[1]+41.0*x[3, 1] <= 30, y[4]-y[1]+38.0*x[4, 1] <= 30, y[5]-y[1]+32.0*x[5, 1] <= 30, y[6]-y[1]+32.0*x[6, 1] <= 30, y[7]-y[1]+40.0*x[7, 1] <= 30, y[8]-y[1]+39.0*x[8, 1] <= 30, y[1]-y[2]+39.0*x[1, 2] <= 30, y[3]-y[2]+37.0*x[3, 2] <= 30, y[4]-y[2]+36.0*x[4, 2] <= 30, y[5]-y[2]+40.0*x[5, 2] <= 30, y[6]-y[2]+38.0*x[6, 2] <= 30, y[7]-y[2]+34.0*x[7, 2] <= 30, y[8]-y[2]+33.0*x[8, 2] <= 30, y[1]-y[3]+41.0*x[1, 3] <= 30, y[2]-y[3]+37.0*x[2, 3] <= 30, y[4]-y[3]+42.0*x[4, 3] <= 30, y[5]-y[3]+42.0*x[5, 3] <= 30, y[6]-y[3]+40.0*x[6, 3] <= 30, y[7]-y[3]+39.0*x[7, 3] <= 30, y[8]-y[3]+38.0*x[8, 3] <= 30, y[1]-y[4]+38.0*x[1, 4] <= 30, y[2]-y[4]+36.0*x[2, 4] <= 30, y[3]-y[4]+42.0*x[3, 4] <= 30, y[5]-y[4]+39.0*x[5, 4] <= 30, y[6]-y[4]+39.0*x[6, 4] <= 30, y[7]-y[4]+34.0*x[7, 4] <= 30, y[8]-y[4]+36.0*x[8, 4] <= 30, y[1]-y[5]+32.0*x[1, 5] <= 30, y[2]-y[5]+40.0*x[2, 5] <= 30, y[3]-y[5]+42.0*x[3, 5] <= 30, y[4]-y[5]+39.0*x[4, 5] <= 30, y[6]-y[5]+32.0*x[6, 5] <= 30, y[7]-y[5]+40.0*x[7, 5] <= 30, y[8]-y[5]+39.0*x[8, 5] <= 30, y[1]-y[6]+32.0*x[1, 6] <= 30, y[2]-y[6]+38.0*x[2, 6] <= 30, y[3]-y[6]+40.0*x[3, 6] <= 30, y[4]-y[6]+39.0*x[4, 6] <= 30, y[5]-y[6]+32.0*x[5, 6] <= 30, y[7]-y[6]+38.0*x[7, 6] <= 30, y[8]-y[6]+37.0*x[8, 6] <= 30, y[1]-y[7]+40.0*x[1, 7] <= 30, y[2]-y[7]+34.0*x[2, 7] <= 30, y[3]-y[7]+39.0*x[3, 7] <= 30, y[4]-y[7]+34.0*x[4, 7] <= 30, y[5]-y[7]+40.0*x[5, 7] <= 30, y[6]-y[7]+38.0*x[6, 7] <= 30, y[8]-y[7]+33.0*x[8, 7] <= 30, y[1]-y[8]+39.0*x[1, 8] <= 30, y[2]-y[8]+33.0*x[2, 8] <= 30, y[3]-y[8]+38.0*x[3, 8] <= 30, y[4]-y[8]+36.0*x[4, 8] <= 30, y[5]-y[8]+39.0*x[5, 8] <= 30, y[6]-y[8]+37.0*x[6, 8] <= 30, y[7]-y[8]+33.0*x[7, 8] <= 30

(15)

batasan1 := seq(a[i] <= y[i], i = N);

0 <= y[1], 0 <= y[2], 0 <= y[3], 0 <= y[4], 0 <= y[5], 0 <= y[6], 0 <= y[7], 0 <= y[8]

(16)

batasan2 := seq(y[i] <= b[i], i = N);

y[1] <= 30, y[2] <= 30, y[3] <= 30, y[4] <= 30, y[5] <= 30, y[6] <= 30, y[7] <= 30, y[8] <= 30

(17)

binaryvariables = {seq(seq(x[i, j], i = `minus`(N, {j})), j = N)};

binaryvariables = {x[1, 2], x[1, 3], x[1, 4], x[1, 5], x[1, 6], x[1, 7], x[1, 8], x[2, 1], x[2, 3], x[2, 4], x[2, 5], x[2, 6], x[2, 7], x[2, 8], x[3, 1], x[3, 2], x[3, 4], x[3, 5], x[3, 6], x[3, 7], x[3, 8], x[4, 1], x[4, 2], x[4, 3], x[4, 5], x[4, 6], x[4, 7], x[4, 8], x[5, 1], x[5, 2], x[5, 3], x[5, 4], x[5, 6], x[5, 7], x[5, 8], x[6, 1], x[6, 2], x[6, 3], x[6, 4], x[6, 5], x[6, 7], x[6, 8], x[7, 1], x[7, 2], x[7, 3], x[7, 4], x[7, 5], x[7, 6], x[7, 8], x[8, 1], x[8, 2], x[8, 3], x[8, 4], x[8, 5], x[8, 6], x[8, 7]}

(18)

conu := seq(seq(u[i]-u[j]+n*x[i, j] <= n-1, i = `minus`(N, {1, j})), j = `minus`(N, {1}));

u[3]-u[2]+8*x[3, 2] <= 7, u[4]-u[2]+8*x[4, 2] <= 7, u[5]-u[2]+8*x[5, 2] <= 7, u[6]-u[2]+8*x[6, 2] <= 7, u[7]-u[2]+8*x[7, 2] <= 7, u[8]-u[2]+8*x[8, 2] <= 7, u[2]-u[3]+8*x[2, 3] <= 7, u[4]-u[3]+8*x[4, 3] <= 7, u[5]-u[3]+8*x[5, 3] <= 7, u[6]-u[3]+8*x[6, 3] <= 7, u[7]-u[3]+8*x[7, 3] <= 7, u[8]-u[3]+8*x[8, 3] <= 7, u[2]-u[4]+8*x[2, 4] <= 7, u[3]-u[4]+8*x[3, 4] <= 7, u[5]-u[4]+8*x[5, 4] <= 7, u[6]-u[4]+8*x[6, 4] <= 7, u[7]-u[4]+8*x[7, 4] <= 7, u[8]-u[4]+8*x[8, 4] <= 7, u[2]-u[5]+8*x[2, 5] <= 7, u[3]-u[5]+8*x[3, 5] <= 7, u[4]-u[5]+8*x[4, 5] <= 7, u[6]-u[5]+8*x[6, 5] <= 7, u[7]-u[5]+8*x[7, 5] <= 7, u[8]-u[5]+8*x[8, 5] <= 7, u[2]-u[6]+8*x[2, 6] <= 7, u[3]-u[6]+8*x[3, 6] <= 7, u[4]-u[6]+8*x[4, 6] <= 7, u[5]-u[6]+8*x[5, 6] <= 7, u[7]-u[6]+8*x[7, 6] <= 7, u[8]-u[6]+8*x[8, 6] <= 7, u[2]-u[7]+8*x[2, 7] <= 7, u[3]-u[7]+8*x[3, 7] <= 7, u[4]-u[7]+8*x[4, 7] <= 7, u[5]-u[7]+8*x[5, 7] <= 7, u[6]-u[7]+8*x[6, 7] <= 7, u[8]-u[7]+8*x[8, 7] <= 7, u[2]-u[8]+8*x[2, 8] <= 7, u[3]-u[8]+8*x[3, 8] <= 7, u[4]-u[8]+8*x[4, 8] <= 7, u[5]-u[8]+8*x[5, 8] <= 7, u[6]-u[8]+8*x[6, 8] <= 7, u[7]-u[8]+8*x[7, 8] <= 7

(19)

Sol := Optimization[LPSolve](z, {conD, conTW, conV, conu, conx, conz, batasan1, batasan2}, binaryvariables = {seq(seq(x[i, j], i = `minus`(N, {j})), j = N)})

Error, (in Optimization:-LPSolve) no feasible integer point found; use feasibilitytolerance option to adjust tolerance

 

X := eval(Matrix(n, symbol = x), {Sol[2][], seq(x[i, i] = 0, i = 1 .. n)})

Error, invalid input: eval expects its 2nd argument, eqns, to be of type {integer, equation, set(equation)}, but received {Sol[2][], seq(x[i, i] = 0, i = 1 .. n)}

 

f := [1, 5, 6, 3, 2, 8, 7, 4, 1];

[1, 5, 6, 3, 2, 8, 7, 4, 1]

(20)

add(c[f[i], f[i+1]], i = 1 .. nops(f)-1);

13.64

(21)

``

Download dataayosemangat.mw

1. The general case (the first output) is fine, but the second output is wrong. Surprisingly, replacing series(ee, ...) with series(sqrt(1/x)*erf(sqrt(x)), ...) does something different:

ee := sqrt(1/x)*erf(sqrt(x));

series(ee, x = 0, 2);
% assuming x < 0;
             -2/sqrt(Pi)+2*x/(3*sqrt(Pi))+O(x^2)

series(ee, x = 0, 2) assuming x < 0; # wrong
              2/sqrt(Pi)-2*x/(3*sqrt(Pi))+O(x^2)

series(sqrt(1/x)*erf(sqrt(x)), x = 0, 2) assuming x < 0;
Error, (in assuming) when calling '`series/signum`'. Received: 'no series at 0'

2. Should Maple be able to handle nested discontinuities? The constant term here is incorrect, so the remainder is not O(z^2).

series(ln(I*ln(-1+z)), z = 0, 2);
         ln(Pi)-(1/2*I)*csgn(ln(-1+z))*(csgn(I*(-1+z))+1)*Pi-I*csgn(I*(-1+z))*z/Pi+O(z^2)

3. This is wrong (it would mean that I^z is unbounded for large positive z):

series(a^z, z = infinity, 2);
                exp((1/2*I)*z*(1-signum(a))*Pi)*abs(a)^z

 

I have a series with an integral inside the series.  I have worked the problem 2 different ways using sum vs Sum.  The integration variables are independent of the series variables so swapping the order of operation should not matter, but in the case, (S2), I do get a difference & I do not understand why.  The explanation fo INERT vs ACTIVE I do not think explains this.  The reason why I say this is because S1 the ACTIVE sum concurs with the INERT expressions S3 & S4.  S2 is swapping the order of operation for the ACTIVE sum does not yield the same result as the other 3 cases.  Why is this?  I am at a loss so including examples would be helpful to me.

swapping_orders_of_operation.mw

I can't undersand, how I can create linear actuator model in MapleSim and to add it in 3D-visualized project. (Like this, for example: http://www.hiwin.com/images/lai2_series_linear_actuator_large.jpg).

Can anyone share an example of such task or explain to me?

It is possible without electric motor and transmission.
An example of translational motion with a certain force will suffice.

Thank you in advance for your cooperation!

First 837 838 839 840 841 842 843 Last Page 839 of 2385