Advanced Search

MaplePrimes Questions

hello guys, i really wondered if there any way to ...

Asked by: shimaa sadk 45
rootfinding fsolve
June 08 2020
0  4

it really confused, with generating different sample some times the command fsolve is working and some times don't . is there any other way to solve this system of equations. because I am going to generate 1000 samples and for each sample, I supposed to solve this system.
 

``

> 

with(LinearAlgebra):
NULL

 

 
> 

``

> 

f[1] := VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`+`(n_x, n_y), 1/theta1), VectorCalculus:-`-`(sum(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`+`(theta2, -1), x[t]), 1/VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), x[t]), 1)), t = 1 .. n_x))), VectorCalculus:-`-`(sum(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`+`(2, R_x[t]), x[t]), 1/VectorCalculus:-`*`(VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), x[t]), 1)^VectorCalculus:-`+`(1, VectorCalculus:-`-`(1/theta2)), VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), x[t]), 1)^(1/theta2), -1), theta3))), t = 1 .. n_x))), VectorCalculus:-`-`(sum(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`+`(theta2, -1), y[T]), 1/VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), y[T]), 1)), T = 1 .. n_y))), VectorCalculus:-`-`(sum(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`+`(2, R_y[T]), y[T]), 1/VectorCalculus:-`*`(VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), y[T]), 1)^VectorCalculus:-`+`(1, VectorCalculus:-`-`(1/theta2)), VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(theta1, theta2), y[T]), 1)^(1/theta2), -1), theta4))), T = 1 .. n_y))):
> 

N_x := 30:
> 

n_x := 30:
> 

NULL
> 

R_x := [seq(0, nl = 1 .. 30)];

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

  

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

 (1)
> 

`λλ` := 5; -1; QQ := 1; -1; `αα1` := 2; -1; `αα2` := 4

4

 (2)
> 

W := GenerateUniform(n_x, 0, 1):
``
> 

fsolve({f[1], f[2], f[3], f[4]}, {theta1 = 0 .. infinity, theta2 = 0 .. infinity, theta3 = .5 .. infinity, theta4 = .5 .. infinity})

 

 (3)
> 

``

``


 

Download ask_maple.mw

why this application don't work ?...

Asked by: panke 140
plot ode proc
June 08 2020
0  8

plotproc.mw

apparently non-existential procedure called by ran...

Asked by: Adam Ledger 360
module showstat
June 08 2020
0  5

I just need someone to explain to me why the following occurs:

 

 

 

> 

> 

showstat(rand)


rand := proc(r::{posint, numeric .. numeric}, $)
   1   if nargs = 0 then
   2     RandomTools:-MersenneTwister:-GenerateInteger()
       else
   3     try
   4       if type(r,('numeric') .. ('numeric')) then
   5         if rhs(r) < lhs(r) then
   6           error "invalid range, end points are out of order"
             elif not type(r,('integer') .. ('integer')) then
   7           return () -> RandomTools:-Generate(float(('range') = r,('method') = ('uniform')))
             end if
           end if;
   8       return RandomTools:-MersenneTwister:-NewGenerator(('range') = r)
         catch :
   9       error
         end try
       end if
end proc

  
> 

NULL
> 

showstat(RandomTools:-MersenneTwister:-GenerateInteger)


RandomTools:-MersenneTwister:-GenerateInteger := proc({range::{integer, integer .. integer} := 1000000000000}, $)
local bits, divisor, offset;
   1   if nargs = 0 then
   2     return MTKernelInterface(4)
       end if;
   3   if type(range,'integer') then
   4     divisor := range;
   5     offset := 0
       else
   6     if rhs(range) < lhs(range) then
   7       error "illegal range, end points out of order"
         end if;
   8     if lhs(range) = rhs(range) then
   9       return lhs(range)
         end if;
  10     offset := lhs(range);
  11     divisor := rhs(range)-offset+1
       end if;
  12   if divisor = 0 then
  13     return offset
       end if;
  14   bits := MTKernelInterface(5,divisor);
  15   MTKernelInterface(6,divisor,bits)+offset
end proc

  
> 

NULL
> 

showstat(MTKernelInterface)

Error, (in showstat) procedure name expected

  
> 

``
> 

``

 

 

Download rand_showstat.mw

Saddle point in maple...

Asked by: radaar 202
calculus saddle-point
June 08 2020
0  8

Does maple provide any method for finding the saddle point of a function? Thank you

How do I make nested for loop in Maple? ...

Asked by: abdgafartunde 20
programming linear-algebra homework
June 08 2020
0  8

I am trying to write maple program for Gauss Siedel method. I needed two for loop. One for running the iterative method and the other for checking where (the number of iteration required) the results converge.

Please, I need help. 

 

Thanks. 

int() not evaluating (Maple 2020)...

Asked by: stonsjen 10
integration syntax 2dmath
June 08 2020
2  3

Hi guys! I'm having trouble using int with a few expressions, for example:

is returning

however simpler expressions such as

int(x^2 + x, x)

or int(x/2 + 1/2, x)

are working correctly and returning 

respectively. I'm using Maple 2020 in worksheet mode using Math 2D Input/Output.

Here's my screen, any help would be much appreciated thanks!

Solving a Double Inequality?...

Asked by: paulmcquad 60
solve homework inequality
June 07 2020
0  7

Hi,

I'm trying to solve for Double Inequality. (1.1) and (1.3) combined.

Thanks in Advance
 

Example: Solving a Double Inequality

 

-3 <= 6*x-1

0 <= 6*x+2

(1.1)

"(->)"

[[-1/3 <= x]]

(1.2)

6*x-1 < 3

6*x < 4

(1.3)

"(->)"[[x < 2/3]]

"ex1(x):=6 x-1"

proc (x) options operator, arrow, function_assign; 6*x-1 end proc

(1.4)

"->"

 

"->"

"ex2(y):=3"

proc (y) options operator, arrow, function_assign; 3 end proc

(1.5)

"ex3(y):=-3"

proc (y) options operator, arrow, function_assign; -3 end proc

(1.6)

NULL


 

Download Double_Inequality.mw

Easier way to solving an Equation Involving an Abs...

Asked by: paulmcquad 60
solve abs
June 07 2020
1  3

Hi,

I was wondering is there an easier way for me to solve an Equation Involving an Absolute Value?


 

Example: Solving an Equation Involving an Absolute Value

 

x^2-3*x = -4*x+6

x^2-3*x = -4*x+6

(1.1)

"(->)"

[[x = 2], [x = -3]]

(1.2)

-x^2+3*x = -4*x+6

-x^2+3*x = -4*x+6

(1.3)

"(->)"

[[x = 1], [x = 6]]

(1.4)

"ex7(x):=|x^(2)-3 x|+4 x-6"

proc (x) options operator, arrow, function_assign; abs(x^2-3*x)+4*x-6 end proc

(1.5)

"->"

 

``

NULL


 

Download Absolute_Value.mw

Maple doesn' work at document mode...

Asked by: mesutc 0
document windows10
June 07 2020
0  0

Mapple doesn' work at document mode.How can i solve the problem.(My Os is Win 10 )

Realization of BVP[midrich]...

Asked by: Dwightun 5
June 07 2020
1  1

Hi!
I am trying to solve ODE with bvp[midrich], and I need to prove the validity of this method. 
Where can i read about realization of bvp[midrich]. What scheme does it use? Аre there any articles on this subject?
 

 

Truncation order of "series" function...

Asked by: greatpet 130
series order truncation
June 07 2020
1  1

In Maple 2019.2, I'm trying to compute series expansions. For the input
series(1/z^4/(1-z), z=0, 1);
Maple gives me a series of the expression around z=0, truncated to 1 term, i.e.
z^(-4) + O(z^(-3))

But with the input
series(1/z^4/(1-z), z=0, 2);
Maple gives me the result
z^(-4) + z^(-3) + z^(-2) + z^(-1) + 1 + z + O(z^2)

Apparently Maple has decided not to give a series with 2 terms, but a series up to O(z^2). I'm surprised that the third argument of "series" is not treated in a consistent manner. I've restarted the Maple server between the two computations, so it's not a memory effect.

In the above example, this is just a minor annoyance, but I'm actually dealing with huge expressions, and for performance reasons it's important to truncate the series expansion at the desired order and stop further computations. Has anyone encountered this problem?

calculate integral based on data set...

Asked by: John Eric 80
integration homework numeric
June 07 2020
0  6

Hello everyone,

Please how can I caluclate an integrale based on data set (Q) usiing any integral method

thanks


 

> 

restart:
> 

NULL
> 

NULL
> 

NULL
> 

Q := [21521.7, 11966.6, 7859.78, 7947.54, 9667.45, 10206.8, 10934.0, 10925.0, 11151.5, 10964.6, 11009.0, 11087.0, 11164.3, 11215.7, 11272.7, 11319.0, 11354.6, 11386.2, 11413.5, 11436.3, 11456.0, 11473.2, 11488.2, 11501.1, 11512.7, 11523.0, 11532.2, 11540.3, 11547.7, 11554.3, 11560.3, 11565.8, 11570.7, 11575.3, 11579.4, 11583.2, 11586.7, 11590.0, 11592.9, 11595.7, 11598.2, 11600.6]

[21521.7, 11966.6, 7859.78, 7947.54, 9667.45, 10206.8, 10934.0, 10925.0, 11151.5, 10964.6, 11009.0, 11087.0, 11164.3, 11215.7, 11272.7, 11319.0, 11354.6, 11386.2, 11413.5, 11436.3, 11456.0, 11473.2, 11488.2, 11501.1, 11512.7, 11523.0, 11532.2, 11540.3, 11547.7, 11554.3, 11560.3, 11565.8, 11570.7, 11575.3, 11579.4, 11583.2, 11586.7, 11590.0, 11592.9, 11595.7, 11598.2, 11600.6]

 (1)
> 

``

NULL


 

Download simpson.mw

How to make a vector- valued function in Maple ?...

Asked by: janhardo 335
vector-calculus
June 07 2020
0  17

I think there is a special needed package for this : student vector-calculus ?

 

How to get the name a colour in this case?...

Asked by: mmcdara 6149
colortools
June 07 2020
0  2

Hi, 

I ran this command
c := NeutralSpread(Color("RGB", "Yellow"), 10);
and I wonder how I could get the  NearestNameColor of all elements in list  c ?


Thanks in advance

how to calculate a permutation raised by a power i...

Asked by: miasene 20
homework permutation do-my-homework-for-me
June 06 2020
2  6 
example let P be a permutation: P: = Matrix ([[1,2,3], [4,5,6]])
calculate P ^ 4567 knowing that the order of P = 20
First 425 426 427 428 429 430 431 Last Page 427 of 2308