shimaa sadk

45 Reputation

5 Badges

3 years, 118 days

MaplePrimes Activity


These are questions asked by shimaa sadk

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

 

Ask_maple.pdf

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/Ask_maple.mw .
 

Download Ask_maple.mw

 

hi, I am trying to evaluate the following integrations (Note that; the values of alpha, lambda_1 and lambda_2) are random and they change with every iteration 

`λ1_B`[so] := 10.94:

``

W[1] := evalf[4](int(`α_B`[so]^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[2] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[3] := evalf[4](int(1/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[4] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[5] := evalf[4](int(1/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[6] := evalf[4](int((z^(`λ2_B`[so]/`λ1_B`[so]))^2*ln(z)^2/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[7] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[8] := evalf[4](int(1/((z-1+`α_B`[so])^2*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^3), z = 1 .. infinity)); W[9] := evalf[4](int(z^(`λ2_B`[so]/`λ1_B`[so])*ln(z)/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity)); W[10] := evalf[4](int(1/((z-1+`α_B`[so])^4*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])), z = 1 .. infinity)); W[11] := evalf[4](int(1/((z-1+`α_B`[so])^3*(z^(`λ2_B`[so]/`λ1_B`[so])-1+`α_B`[so])^2), z = 1 .. infinity))

int(79.977249/((z+7.943)^2*(z^.1408592322+7.943)), z = 1 .. infinity)

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

0.1744e-1

 

Warning,  computation interrupted

 

0.3348e-2

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

Warning,  computation interrupted

 

``

``


 

Download ask_maple.mw

 

1 2 3 4 5 6 7 Page 1 of 7