shimaa sadk

45 Reputation

5 Badges

7 years, 211 days

MaplePrimes Activity


These are questions asked by shimaa sadk

i have a big code some times maple get stuck in a certain evaluation  
How can I issue a command to Maple to stop an command  if the time expands for a specific time?

download the file

Download ask.mw 

 

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


 

``

lambda := .3:

omega := lambda+mu+xi:

alpha := 2*sqrt(lambda*mu)

.9165151390

(1)

``

B[1] := BesselI(k-1, alpha*(u-y))

BesselI(k-1, .9165151390-.9165151390*y)

(2)

B[2] := BesselI(k+1, alpha*(u-y))

BesselI(k+1, .9165151390-.9165151390*y)

(3)

``

F := evalf(Int(sum((B[1]-B[2])*exp(-omega*(u-y)), k = 1 .. infinity), y = 0 .. u))

Int(sum((BesselI(k-1., .9165151390-.9165151390*y)-1.*BesselI(k+1., .9165151390-.9165151390*y))*exp(-1.200000000+1.200000000*y), k = 1 .. infinity), y = 0. .. 1.)

(4)

``

``

``


 

Download int.mw

1 2 3 4 5 6 7 Page 1 of 8