Tags are words are used to describe and categorize your content. Combine multiple words with dashes(-), and seperate tags with spaces.
here my problem i work on ;
i want maple to end the enternal loop if it can't solve the equation after some iteration maple couldnt solve the equation and get stuck
Error, cannot determine if this expression is true or false: 0 < K
Download Q3.mw
Please log-in to your MaplePrimes account.
Wrong Email/Password. Please try again.
Error occurred during PDF generation. Please refresh the page and try again
You must be logged into your Facebook account in order to share via Facebook.
Click the button below to share this on Google+. A new window will open.
You must be logged in to your Twitter account in order to share. Click the button below to login (a new window will open.)
![with(LinearAlgebra); with(VectorCalculus); with(Student[LinearAlgebra]); with(SignalProcessing); with(Statistics); with(stats)](/view.aspx?sf=225291_question/b5dc63b157e1038c25adc9f0112a4fa2.gif)
![f[1] := n/R+sum(x[i], i = 1 .. n)-(sum((2+a[i])*x[i]*exp(R*x[i])/(exp(R*x[i])-1+Q), i = 1 .. n))](/view.aspx?sf=225291_question/2822ce3b5cbd6b010eab7b10b5404a93.gif)
![f[2] := m/S+sum(y[j], j = 1 .. m)-(sum((2+b[j])*y[j]*exp(y[j]*S)/(exp(y[j]*S)-1+Q), j = 1 .. m))](/view.aspx?sf=225291_question/5816b27bb0d74839d5ace3ab33efcdf2.gif)
![f[3] := (n+m+sum(a[i], i = 1 .. n)+sum(b[j], j = 1 .. m))/Q-(sum((2+a[i])/(exp(R*x[i])-1+Q), i = 1 .. n))-(sum((2+b[j])/(exp(y[j]*S)-1+Q), j = 1 .. m))](/view.aspx?sf=225291_question/544050b883232ecd9c49e440791b7065.gif)
![E1[1] := 0.5e-1](/view.aspx?sf=225291_question/002723aed3882e6945ef492b764251c6.gif)
![E2[1] := 0.5e-1](/view.aspx?sf=225291_question/d07915f70c5a27f9094e9841de9ff6a8.gif)
![E3[1] := 0.5e-1](/view.aspx?sf=225291_question/abb20ac48b75aef3c9d83f6d6e12dfdd.gif)
![a := [seq(0, i = 1 .. 21), 2, 2, 1, seq(0, i = 1 .. 21)]; b := [seq(0, i = 1 .. 21), 2, 2, 1, seq(0, i = 1 .. 21)]](/view.aspx?sf=225291_question/f3d4afb4ea6256f00bb372faa7cf1526.gif)
![so := 1; while so < K+1 do W := GenerateUniform(n, 0, 1); for iii to n do vv[iii] := W[iii]^(1/(iii+sum(a[jjj], jjj = n-iii+1 .. n))) end do; for sss to n do uu[sss] := 1-product(vv[n-jjj+1], jjj = 1 .. sss); x[sss] := fsolve(1-`α_α`/(exp(lambda1*t)-1+`α_α`) = uu[sss], t = 0 .. infinity) end do; U := GenerateUniform(m, 0, 1); for ii to m do v[ii] := U[ii]^(1/(ii+sum(b[jj], jj = m-ii+1 .. m))) end do; for ss to m do u[ss] := 1-product(v[m-jj+1], jj = 1 .. ss); y[ss] := fsolve(1-`α_α`/(exp(lambda2*t)-1+`α_α`) = u[ss], t = 0 .. infinity) end do; c := describe[quartile[1]]([seq(x[i], i = 1 .. n)]); cc := describe[quartile[3]]([seq(x[i], i = 1 .. n)]); L := describe[quartile[1]]([seq(y[i], i = 1 .. m)]); LL := describe[quartile[3]]([seq(y[i], i = 1 .. m)]); R[1] := fsolve(9*exp(R*c)-exp(R*cc) = 8, R = 0 .. infinity); S[1] := fsolve(9*exp(S*L)-exp(S*LL) = 8, S = 0 .. infinity); Q[1] := 3*(exp(R[1]*c)-1+(exp(S[1]*L)-1))*(1/2); for h to 40 while `and`(`and`(`and`(`and`(`and`(abs(E1[h]) > 0.5e-3, abs(E2[h]) > 0.5e-3), abs(E3[h]) > 0.5e-3), Q[h] > 2), S[h] > 0), R[h] > 0) do printlevel := 2; Q[h+1] := fsolve(eval(f[3], {R = R[h], S = S[h]}) = 0, Q = 0 .. infinity); S[h+1] := fsolve(eval(f[2], Q = Q[h+1]) = 0, S = 0 .. infinity); R[h+1] := fsolve(eval(f[1], Q = Q[h+1]) = 0, R = 0 .. infinity); if fsolve(eval(f[3], {R = R[h], S = S[h]}) = 0, Q = 0 .. infinity) = NULL then break end if; if fsolve(eval(f[2], Q = Q[h+1]) = 0, S = 0 .. infinity) = NULL then break end if; if fsolve(eval(f[1], Q = Q[h+1]) = 0, R = 0 .. infinity) = NULL then break end if; E1[h+1] := abs(R[h+1]-R[h]); E2[h+1] := abs(S[h+1]-S[h]); E3[h+1] := abs(Q[h+1]-Q[h]); KK := Matrix([[R[h]], [S[h]], [Q[h]]]) end do; A[so] := Determinant(KK[1]); B[so] := Determinant(KK[2]); C[so] := Determinant(KK[3]); so := so+1 end do](/view.aspx?sf=225291_question/35f7739c59ea5596c0e69e88dcbfa39c.gif)
