Maple 2015 Questions and Posts

These are Posts and Questions associated with the product, Maple 2015

Hellow dears!!!

Hope everyone is fine with everything. I want the animation of the numerical solution of ODE i,e., f4 for delta=[0,1.5]. Please see the attachment and fix my problem. I Shall be very thankful to you.

Graph.mw

Special request to  acer 13834@Carl Love ,Preben Alsholm 10271

Irrational numbers: numbers that cannot be represented as a ratio of integers. The decimal form of a rational number is non-repeating and non-terminating.

Change to:

Irrational numbers: numbers that cannot be represented as a ratio of integers. The decimal form of an irrational number is non-repeating and non-terminating. 

or change to

Irrational numbers can be represented by decimal fractions in which the digits go on forever without ever repeating a pattern.  See Downing, Douglas. Dictionary of Mathematics Terms. 2nd ed. Hauppauge, NY: Barron's Ed. Series, Inc., 1995, p. 176).

It is a truth universally acknowledged, that a single man in possession of an algorithm for calculating a Groebner basis, must be in want of an algorthim for calculating a reduced Groebner basis.

It seems odd that i can't find something in the Groebner package - if there isn't something there, I assume that there is a well known piece of code for doing this!

Maple 2015

Using with(combinat) the permutation of {a,b,c} is determined.

>restart:
>with(combinat):
>permute({a, b, c})
                  [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]

The tree diagram of this permutation is

    

In Maple, using with(combinat) and with(GraphTheory), when I attempt to draw the permutation I get the following error:

>L := permute({a, b, c});
       L := [[a, b, c], [a, c, b], [b, a, c], [b, c, a], [c, a, b], [c, b, a]]
>DrawGraph(L);
  Error, invalid input: GraphTheory:-DrawGraph expects its 1st argument, H, to be of type       {GRAPHLN, list(GRAPHLN), set(GRAPHLN)}, but received [[a, b, c], [a, c, b], [b, a, c], [b, c,      a], [c,   a, b], [c, b, a]]

On Maple, again using with(combinat) and with(GraphTheory) the command permute(3) is used.  The results are manually configured as node-connection lines.  A fair representation of the tree diagram is configured by Maple, although the diagram has numeric instead of alpha configurations, and the a,b,c structure shown above is not easily recognized.

Any suggestions on developing a procedure that will graph (draw) an alpha-labeled permutation welcomed.  Thanks!  WC44_Permutation_Graph.mw

Dearz

Hope you would be fine with everything. I try to solve the following linear system of equations via fsolve command but the solution doesn't satisfied the system please see and put your valueable comments. Waiting your positive response.

-5.7167551941125971285 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 4.7759348859301130832 d[1, 3]

   + 82.882747548740738074 d[1, 4]

   + 1.5473302855836067493 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.2926823766365742120 d[2, 3]

   - 22.433527600870893213 d[2, 4]

   - 11.906076336447024126 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 9.9466643924764099316 d[3, 3]

   + 172.61685795222431091 d[3, 4]

   + 153.42462622364681378 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 222.04914007834331471 d[4, 3]

   - 2162.1913920527683546 d[4, 4] = 0
-6.3505370802317673052 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 5.4097167720492832599 d[1, 3]

   + 54.802782951629695640 d[1, 4]

   + 1.7188733853924696263 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.4642254764454370890 d[2, 3]

   - 14.833240696164645293 d[2, 4]

   - 13.226030621801645811 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 11.266618677831031617 d[3, 3]

   + 114.13574573628827681 d[3, 4]

   + 107.19584752215150208 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 175.82036137684800302 d[4, 3]

   - 1136.3239123361047712 d[4, 4] = 0
-6.7642088272251297212 d[1, 1] - 0.23520507704562101132 d[1, 2]

   - 5.8233885190426456759 d[1, 3]

   + 34.632657184619275137 d[1, 4]

   + 1.8308401918550417305 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.5761922829080091932 d[2, 3] - 9.373876878125528749 d[2, 4]

   - 14.087569594645296643 d[3, 1]

   - 0.48985298599265354856 d[3, 2]

   - 12.128157650674682449 d[3, 3]

   + 72.128164697121390121 d[3, 4]

   + 77.022155175221117487 d[4, 1]

   - 17.156128463674125233 d[4, 2]

   + 145.64666902991761843 d[4, 3]

   - 601.11088029977885095 d[4, 4] = 0
-1.5473302855836067487 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.2926823766365742115 d[1, 3]

   + 22.433527600870893218 d[1, 4]

   + 1.5473302855836067493 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.2926823766365742120 d[2, 3]

   - 22.433527600870893213 d[2, 4]

   - 7.7366514279180337465 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 6.4634118831828710599 d[3, 3]

   + 112.16763800435446606 d[3, 4]

   + 104.66490008068725185 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 181.31392067374404557 d[4, 3]

   - 1455.2623850848598494 d[4, 4] = 0
-1.7188733853924696257 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.4642254764454370885 d[1, 3]

   + 14.833240696164645297 d[1, 4]

   + 1.7188733853924696263 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.4642254764454370890 d[2, 3]

   - 14.833240696164645293 d[2, 4]

   - 8.5943669269623481316 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 7.3211273822271854450 d[3, 3]

   + 74.166203480823226458 d[3, 4]

   + 53.030427038219525869 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 129.67944763127631958 d[4, 3]

   - 668.89639482661771723 d[4, 4] = 0
-1.8308401918550417299 d[1, 1] - 0.06366197723675813430 d[1, 2]

   - 1.5761922829080091926 d[1, 3] + 9.373876878125528754 d[1, 4]

   + 1.8308401918550417305 d[2, 1]

   + 0.063661977236758134308 d[2, 2]

   + 1.5761922829080091932 d[2, 3] - 9.373876878125528749 d[2, 4]

   - 9.1542009592752086523 d[3, 1]

   - 0.31830988618379067154 d[3, 2]

   - 7.8809614145400459657 d[3, 3]

   + 46.869384390627643742 d[3, 4]

   + 19.328418292985322519 d[4, 1]

   - 19.162255148264198426 d[4, 2]

   + 95.977438886042116228 d[4, 3]

   - 305.71973224709969080 d[4, 4] = 0
 7.0561523113686303394 d[1, 1] - 1.9098593171027440292 d[1, 2]

    + 14.695589579779606456 d[1, 3]

    - 96.471127562654332340 d[1, 4]

    - 2.3520507704562101132 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.8985298599265354856 d[2, 3]

    + 32.157042520884777447 d[2, 4]

    + 16.464355393193470792 d[3, 1]

    - 4.4563384065730694016 d[3, 2]

    + 34.289709019485748399 d[3, 3]

    - 225.09929764619344213 d[3, 4]

    - 96.434081588704614639 d[4, 1]

    + 26.101410667070835066 d[4, 2]

    - 200.83972425698795490 d[4, 3]

    + 1318.4387433562758754 d[4, 4] = 0
-2.3520507704562101132 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 4.898529859926535486 d[1, 3] + 32.157042520884777450 d[1, 4]

   - 2.3520507704562101132 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 4.8985298599265354856 d[2, 3]

   + 32.157042520884777447 d[2, 4]

   + 7.0561523113686303394 d[3, 1]

   - 1.9098593171027440293 d[3, 2]

   + 14.695589579779606457 d[3, 3] - 96.47112756265433234 d[3, 4]

   - 11.760253852281050559 d[4, 1] + 3.183098861837906715 d[4, 2]

   - 24.49264929963267742 d[4, 3] + 160.7852126044238874 d[4, 4] = 

  1
1.9098593171027440291 d[1, 1] - 1.9098593171027440292 d[1, 2]

   + 9.5492965855137201456 d[1, 3]

   - 36.287327024952136554 d[1, 4]

   - 0.6366197723675813430 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 3.1830988618379067154 d[2, 3]

   + 12.095775674984045518 d[2, 4]

   + 4.4563384065730694010 d[3, 1]

   - 4.4563384065730694016 d[3, 2]

   + 22.281692032865347008 d[3, 3] - 84.67042972488831863 d[3, 4]

   - 26.101410667070835067 d[4, 1]

   + 26.101410667070835066 d[4, 2]

   - 130.50705333535417533 d[4, 3]

   + 495.92680267434586630 d[4, 4] = 0
-0.6366197723675813431 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 3.1830988618379067164 d[1, 3]

   + 12.095775674984045516 d[1, 4]

   - 0.6366197723675813430 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 3.1830988618379067154 d[2, 3]

   + 12.095775674984045518 d[2, 4]

   + 1.9098593171027440288 d[3, 1]

   - 1.9098593171027440293 d[3, 2] + 9.549296585513720146 d[3, 3]

   - 36.287327024952136560 d[3, 4] - 3.183098861837906717 d[4, 1]

   + 3.183098861837906715 d[4, 2] - 15.91549430918953358 d[4, 3]

   + 60.47887837492022764 d[4, 4] = 1
-1.4491448767744190950 d[1, 1] - 1.9098593171027440292 d[1, 2]

   + 6.1902923916365570215 d[1, 3]

   - 11.964006709004497915 d[1, 4]

   + 0.4830482922581396984 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 2.0634307972121856740 d[2, 3] + 3.988002236334832639 d[2, 4]

   - 3.381338045806977889 d[3, 1] - 4.4563384065730694016 d[3, 2]

   + 14.444015580485299718 d[3, 3] - 27.91601565434382847 d[3, 4]

   + 19.804979982583727629 d[4, 1]

   + 26.101410667070835066 d[4, 2]

   - 84.600662685699612634 d[4, 3] + 163.5080916897281382 d[4, 4] = 

  0
0.4830482922581396984 d[1, 1] + 0.6366197723675813431 d[1, 2]

   - 2.0634307972121856744 d[1, 3] + 3.988002236334832645 d[1, 4]

   + 0.4830482922581396984 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 2.0634307972121856740 d[2, 3] + 3.988002236334832639 d[2, 4]

   - 1.4491448767744190956 d[3, 1]

   - 1.9098593171027440293 d[3, 2]

   + 6.1902923916365570221 d[3, 3] - 11.96400670900449791 d[3, 4]

   + 2.415241461290698491 d[4, 1] + 3.183098861837906715 d[4, 2]

   - 10.317153986060928369 d[4, 3] + 19.94001118167416332 d[4, 4] = 

  1
 11.581726419330485018 d[1, 1] - 3.8605754731101616728 d[1, 2]

    + 27.024028311771131709 d[1, 3]

    - 158.28359439751662858 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 19.221163687741461135 d[3, 1]

    - 6.4070545625804870452 d[3, 2]

    + 44.849381938063409316 d[3, 3]

    - 262.68923706579996884 d[3, 4]

    - 172.31418534244454203 d[4, 1]

    + 57.438061780814847345 d[4, 2]

    - 402.06643246570393142 d[4, 3]

    + 2354.9605330134087411 d[4, 4] = 0
 7.0561523113686303394 d[1, 1] - 2.3520507704562101132 d[1, 2]

    + 16.464355393193470792 d[1, 3]

    - 96.434081588704614639 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 14.695589579779606456 d[3, 1]

    - 4.8985298599265354856 d[3, 2]

    + 34.289709019485748399 d[3, 3]

    - 200.83972425698795490 d[3, 4]

    - 96.471127562654332340 d[4, 1]

    + 32.157042520884777447 d[4, 2]

    - 225.09929764619344213 d[4, 3]

    + 1318.4387433562758753 d[4, 4] = 0
 1.9098593171027440291 d[1, 1] - 0.6366197723675813430 d[1, 2]

    + 4.4563384065730694010 d[1, 3]

    - 26.101410667070835067 d[1, 4]

    - 1.9098593171027440292 d[2, 1]

    + 0.63661977236758134308 d[2, 2]

    - 4.4563384065730694016 d[2, 3]

    + 26.101410667070835066 d[2, 4]

    + 9.5492965855137201456 d[3, 1]

    - 3.1830988618379067154 d[3, 2]

    + 22.281692032865347008 d[3, 3]

    - 130.50705333535417533 d[3, 4]

    - 36.287327024952136554 d[4, 1]

    + 12.095775674984045518 d[4, 2]

    - 84.670429724888318626 d[4, 3]

    + 495.92680267434586626 d[4, 4] = 0
-1.4491448767744190950 d[1, 1] + 0.4830482922581396984 d[1, 2]

   - 3.381338045806977889 d[1, 3] + 19.804979982583727629 d[1, 4]

   - 1.9098593171027440292 d[2, 1]

   + 0.63661977236758134308 d[2, 2]

   - 4.4563384065730694016 d[2, 3]

   + 26.101410667070835066 d[2, 4]

   + 6.1902923916365570215 d[3, 1]

   - 2.0634307972121856740 d[3, 2]

   + 14.444015580485299718 d[3, 3]

   - 84.600662685699612634 d[3, 4]

   - 11.964006709004497917 d[4, 1] + 3.988002236334832638 d[4, 2]

   - 27.91601565434382847 d[4, 3] + 163.50809168972813819 d[4, 4] = 

  0
Sols := solve([seq(`$`(F1[l1, l2], l1 = 2 .. 2^K*M-1), l2 = 2 .. 2^K*M), seq(`$`(F2[l1], l1 = 2 .. 2^K*M)), seq(`$`(F3[l1], l1 = 2 .. 2^K*M)), seq(`$`(F4[l1], l1 = 1 .. 2^K*M))], {seq(`$`(d[l1, l2], l1 = 1 .. 2^K*M), l2 = 1 .. 2^K*M)});
map(evalf, subs(Sols, convert(F4, list)));

 

Dearz!

Hope everyone is fine with everything. I am facing problem to solve the system of PDEs in the attached file. Is there any built-in command to the solve the attached system of PDEs via FEM, FDM, SIMPLER algorithm or some other efficient method? Please try to fix my problem. I am waiting your positive response. Thanks in advance.

PDEs_sol.mw

Hello,

Could you please help me with the following problem? I'm new to Maple and i need some help.

Solve the equation x^3 - a*x + 1 = 0 , in x. Determine the particular solution for a=1,2,... .Graphically represent the polynom that appears in the equation, in a case where the equation has a real root and in a case where the equation has 3 real roots.

Thank you !

Hi everibody 

I work with Maple 2015 under OS-X El Capitan.

Using more than one matrix vector product (either M.V  or MatrixVectorMultiply(M,V)  ; M is a n by p matrix and V a column vector of size p) within the same block of commands generates an error.

Do other people have the same problem ?
Thanks for your feedback.

SomethingGoesWrong.mw


PS : I know I can do this   X . <<1, 1, -1> | <-1, 2, 0>> but this doesn't explain the error I get

 

what are the dynamical system which act on invariant manifold?

Hi,

This sequence of commands works perfectly well

     plotsetup(jpeg, plotoutput=SomeJpegFile);
     plot(x, x=0..1);
     plotsetup(default);


Why this one doesn't create the file SomeJpegFile ?

f := proc()
     plotsetup(jpeg, plotoutput=SomeJpegFile);
     plot(x, x=0..1);
     plotsetup(default);
end proc;

f();


Thanks in advance

Is it a complete set ? How to search matrix?

Dear 

Hope everyone is good. I am face to attaine the converges solution of the attached problem. Please have a look and fix my problem. I am waiting your response

diverges.mw


 

لا شيء

-------------------------------------------------- -------------------------------------------------- -------------------------------------------------- -------------------------

إعادة بدء

مع (LinearAlgebra)

مع (orthopoly)

مع (طالب)

لا شيء

لا شيء

لا شيء

لا شيء

سيل (ألفا): = 2؛  سيل (بيتا): = 1؛  ألفا: = 1.5؛  بيتا: = .5

2

 

1

 

1.5

 

0.5

(1)

n: = 8؛  m: = 8

8

 

8

(2)

 

لا شيء

x [3]: = .611423302089630؛  x [4]: ​​= 1.09446605083631؛  x [5]: = 1.99636816302962؛  x [6]: = 3.38757178455234؛  x [7]: = 5.41873370919121؛  x [8]: = 8.49143699030089

،611423302089630

 

+1.09446605083631

 

1.99636816302962

 

3.38757178455234

 

5.41873370919121

 

8.49143699030089

(3)

# 1 / حساب مصفوفة (A). (طريقة الجمع)

A := array(1 .. n, 1 .. m); for j to m do A[1, j] := evalf(subs(x = 0, L(j-1, 2*x-1))) end do; for j to m do A[2, j] := evalf(subs(x = 0, diff(L(j-1, 2*x-1), x))) end do; for i from 3 to n do for j to m do A[i, j] := evalf(subs(x = x[i], fracdiff(L(j-1, 2*x-1), x, alpha, method = direct))+subs(x = x[i], fracdiff(L(j-1, 2*x-1), x, beta, method = direct))+subs(x = x[i], diff(L(j-1, 2*x-1), x))+subs(x = x[i], L(j-1, 2*x-1))) end do end do

print(`A=`, A)

`A=`, A

(4)

A := convert(A, Matrix)

A := Matrix(8, 8, {(1, 1) = 1., (1, 2) = 2., (1, 3) = 3.500000000, (1, 4) = 5.666666667, (1, 5) = 8.708333333, (1, 6) = 12.88333333, (1, 7) = 18.50972222, (1, 8) = 25.97658730, (2, 1) = 0., (2, 2) = -2., (2, 3) = -6., (2, 4) = -13., (2, 5) = -24.33333333, (2, 6) = -41.75000000, (2, 7) = -67.51666667, (2, 8) = -104.5361111, (3, 1) = 1., (3, 2) = -2.987486314, (3, 3) = -3.301220288, (3, 4) = .5119939327, (3, 5) = 9.171314221, (3, 6) = 23.72035697, (3, 7) = 45.59773916, (3, 8) = 76.72165628, (4, 1) = 1., (4, 2) = -4.549878909, (4, 3) = -1.208865530, (4, 4) = 6.408882482, (4, 5) = 16.03540544, (4, 6) = 27.10075251, (4, 7) = 40.26736031, (4, 8) = 57.11215315, (5, 1) = 1., (5, 2) = -7.181375466, (5, 3) = 6.777193107, (5, 4) = 12.19170970, (5, 5) = 9.600555508, (5, 6) = 7.084730200, (5, 7) = 11.13249218, (5, 8) = 24.60731420, (6, 1) = 1., (6, 2) = -10.92878792, (6, 3) = 28.28352183, (6, 4) = -10.19173665, (6, 5) = -20.04576479, (6, 6) = 9.17677094, (6, 7) = 39.97816692, (6, 8) = 49.07345342, (7, 1) = 1., (7, 2) = -16.09078867, (7, 3) = 78.08969329, (7, 4) = -166.5158779, (7, 5) = 129.0586058, (7, 6) = 104.8307190, (7, 7) = -104.838425, (7, 8) = -111.0119440, (8, 1) = 1., (8, 2) = -23.55908364, (8, 3) = 192.6052140, (8, 4) = -856.8131732, (8, 5) = 2255.610395, (8, 6) = -3256.154493, (8, 7) = 1577.05254, (8, 8) = 2063.443568})

(5)

NULL

# ------------------------------------------------- --------------------------
# 2 / حساب مصفوفة (ب) من قبل أدومين بوليس لمصطلح غير الخطية.

"G(y):=(e)^(y)"

proc (y) options operator, arrow; exp(y) end proc

(6)

"g(x):=evalf(((4*sqrt(x))/(sqrt(Pi)))+(8/(3))*((x^(3/(2)))/(sqrt(Pi)))+2*x+x^(2)+(e)^(x^(2)))"

proc (x) options operator, arrow; evalf(4*sqrt(x)/sqrt(Pi)+(8/3)*x^(3/2)/sqrt(Pi)+2*x+x^2+exp(x^2)) end proc

(7)

#Find أدومين بولي:

for k from 0 to n-1 do AP[k] := evalf(subs(lambda = 0, (diff(G(sum(y[t]*lambda^t, t = 0 .. k)), [`$`(lambda, k)]))/factorial(k))) end do

exp(y[0])

 

y[1]*exp(y[0])

 

y[2]*exp(y[0])+.5000000000*y[1]^2*exp(y[0])

 

y[3]*exp(y[0])+y[2]*y[1]*exp(y[0])+.1666666667*y[1]^3*exp(y[0])

 

y[4]*exp(y[0])+y[3]*y[1]*exp(y[0])+.5000000000*y[2]^2*exp(y[0])+.5000000000*y[2]*y[1]^2*exp(y[0])+0.4166666667e-1*y[1]^4*exp(y[0])

 

y[5]*exp(y[0])+y[4]*y[1]*exp(y[0])+y[3]*y[2]*exp(y[0])+.5000000000*y[3]*y[1]^2*exp(y[0])+.5000000000*y[2]^2*y[1]*exp(y[0])+.1666666667*y[2]*y[1]^3*exp(y[0])+0.8333333333e-2*y[1]^5*exp(y[0])

 

y[6]*exp(y[0])+y[5]*y[1]*exp(y[0])+y[4]*y[2]*exp(y[0])+.5000000000*y[4]*y[1]^2*exp(y[0])+.5000000000*y[3]^2*exp(y[0])+y[3]*y[2]*y[1]*exp(y[0])+.1666666667*y[3]*y[1]^3*exp(y[0])+.1666666667*y[2]^3*exp(y[0])+.2500000000*y[2]^2*y[1]^2*exp(y[0])+0.4166666667e-1*y[2]*y[1]^4*exp(y[0])+0.1388888889e-2*y[1]^6*exp(y[0])

 

y[7]*exp(y[0])+.5000000000*y[3]*y[2]*y[1]^2*exp(y[0])+.5000000000*y[5]*y[1]^2*exp(y[0])+y[5]*y[2]*exp(y[0])+y[6]*y[1]*exp(y[0])+y[4]*y[3]*exp(y[0])+.5000000000*y[3]^2*y[1]*exp(y[0])+.1666666667*y[2]^3*y[1]*exp(y[0])+0.1984126984e-3*y[1]^7*exp(y[0])+y[4]*y[2]*y[1]*exp(y[0])+0.8333333333e-2*y[2]*y[1]^5*exp(y[0])+0.8333333333e-1*y[2]^2*y[1]^3*exp(y[0])+0.4166666667e-1*y[3]*y[1]^4*exp(y[0])+.5000000000*y[3]*y[2]^2*exp(y[0])+.1666666667*y[4]*y[1]^3*exp(y[0])

(8)

NULL

#Find a ماتريسز b ^ (k) و C ^ (k): = A ^ (- 1) * b ^ (k)، ثم ايجاد حل تقريبي Y [k]: = سوم (C ^ (k) [i ] * L [i]، i = 1 .. n ):

# 1) البحث ب (0)

b0 := array(1 .. n, 1 .. m-7); for i to 2 do b0[i, 1] := 0 end do; for i from 3 to n do b0[i, 1] := evalf(subs(x = x[i], g(x[i]))) end do

print(`b0=`, b0)

`b0=`, b0

(9)

b0 := convert(b0, Matrix)

b0 := Matrix(8, 1, {(1, 1) = 0, (2, 1) = 0, (3, 1) = 5.533921684, (4, 1) = 10.78339161, (5, 1) = 69.22208674, (6, 1) = 96372.14332, (7, 1) = 0.5649990671e13, (8, 1) = 0.2063418920e32})

(10)

# 2) البحث عن ج (0)

C0 := LinearSolve(A, b0)

C0 := Matrix(8, 1, {(1, 1) = -0.11474558283495975e27, (2, 1) = -0.6041534517526968e26, (3, 1) = 0.28431046341368933e27, (4, 1) = -0.1109483456679843e28, (5, 1) = 0.2601411410469915e28, (6, 1) = -0.34736953613415415e28, (7, 1) = 0.23829217145639085e28, (8, 1) = -0.634449734180237e27}, datatype = float[8])

(11)

for i to n do k0[i] := C0[i, 1] end do

HFloat(-1.1474558283495975e26)

 

HFloat(-6.041534517526968e25)

 

HFloat(2.8431046341368933e26)

 

HFloat(-1.109483456679843e27)

 

HFloat(2.601411410469915e27)

 

HFloat(-3.4736953613415415e27)

 

HFloat(2.3829217145639085e27)

 

HFloat(-6.34449734180237e26)

(12)

# 3) البحث عن y (0)

y[0] := sum(k0[s]*L(s-1, 2*x-1), s = 1 .. 8)

-HFloat (5.083969685801073e25) -HFloat (1.4661238981264424e26) * س + HFloat (1.2387812172594187e26) * (2 * س 1) ^ 2-HFloat (1.9836944590452831e24) * (2 * س 1) ^ 3 HFloat (5.120751558697758 E25) * (2 * س 1) ^ 4 + HFloat (2.0830079097858884e25) * (2 * س 1) ^ 5 HFloat (2.8586478120802086e24) * (2 * س 1) ^ 6 + HFloat (1.2588288376592004e23) * (2 * س 1) ^ 7

(13)

# -------------------------

#Find b (1)

لا شيء

لا شيء

لا شيء

b1: = أري (1 .. n، 1 .. m-7)؛  ل i تو 2 دو b1 [i، 1]: = 0 إند دو؛  من i إلى n n b1 [i، 1]: = سوبس (x = x [i]، أب [0]) إند دو

برينت (`b1 =`، b1)

`b1 =`، b1

(14)

b1: = كونفيرت (b1، ماتريكس)

b1: = مصفوفة (8، 1، {(1، 1) = 0، (2، 1) = 0، (3، 1) = إكس (هفلوات (-1.3446720400287247e26))، (4، 1) = إكس هفلوت (-1.000132892371102e26))، (5، 1) = إكس (هفلوت (-1.7743764624635952e26))، (6، 1) = إكس (هفلوت (9.701444095568667e26))، (7، 1) = إكس 1.9741498268709318e28))، (8، 1) = إكس (هفلوات (4.2920269682087554e30))})

(15)

لا شيء

# 2) البحث ج (1)

لينيرزولف (A، b1)

المصفوفة ([هفلوات (هفلوات (وندفيند))]، [هفلوت (هفلوات (وندفيند))]، [هفلوات (هفلوات (وندفيند))]، [هفلوات (هفلوات (وندفيند))]، [هفلوات (هفلوات (وندفيند) )، [هفلوات (هفلوات (وندفيند))]، [هفلوات (هفلوات (وندفيند))]، [هفلوات (هفلوات (وندفيند))]])

(16)

لا شيء


 

تحميل jam.mw

Hi everybody,

I want to solve numerically an ode and I get this error (undocumented on the maplesoft web site https://www.maplesoft.com/support/help/errors/....)

Error, (in sol) maximum number of event iterations reached (100) at t=2.6610663

I understand where this error can come from but the help pages don't say anything to fix this.
There is some stuff about round-off that could help but I don't understand how to use it.

I would be grateful if you provide me some help.
Thanks in advance


Download ErrorWithDsolve.mw

 

 

First 21 22 23 24 25 26 27 Last Page 23 of 71