abdulganiy

140 Reputation

6 Badges

8 years, 220 days

MaplePrimes Activity


These are replies submitted by abdulganiy

@acer 

Thank you for your reply. I wanted to write an independent code different from the Maple inbuilt package for numerical. However, met with difficulty. The exact and numerical in my code do not agree.

Kindly see the attached for your perusal and possible modifications and enlightenment
 

restart

``

restart:

``

``

e1:=y[n+1]=y[n]+h*f(n):

``

``

inx:=2:

    h               Num.y                Ex.y             Error y
 0.10        5.1000000000        6.8858944423                1.79
 0.20        5.2100000000        6.4496474644                1.24
 0.30        5.3290000000        6.0739473917               0.745
 0.40        5.4561000000        5.7530324244               0.297
 0.50        5.5904900000        5.4816890703               0.109
 0.60        5.7314410000        5.2551999668               0.476
 0.70        5.8782969000        5.0692966676               0.809
 0.80        6.0304672100        4.9201169227                1.11
 0.90        6.1874204890        4.8041660239                1.38
 1.00        6.3486784401        4.7182818285                1.63

``


 

Download EULER_Correction.mw

.

@Kitonum 

Your suggestion is appreciated.

@acer 

Thank you. The different suggestions are appreciated.

@acer

Your suggestion is very useful. Thank you. It is appreciated.

However, I will appreciate it if I can learn some other ways to get the same effect from you.

Thank you and kind regards.

@Rouben Rostamian  

Thank you so much for your time. It is appreciated. Can you kindly include in the code the Matrix representation of the coefficients of the system of ODEs just immediately after " sys := { seq(eq__i, i=2..n-2) } union bc_left union bc_right: "? It would be appreciated. Thank you in anticipation of your positive response and best regards.

@acer Thank you so much. It is appreciated. 

Thank you it is appreciated

@Rouben Rostamian  

Your Suggestion is  highly apprciated. Thank you and kind regards

@acer Thank you very much. I really appreciate your suggestion and value the time you spent on it.

@Kitonum I really appreciate your suggestion. However, I have applied the "allvalues" command before settling for the line b the results dont just make any sense to me (too long and some contain imaginary) which shouldnt be.

@acer Thank you and kind regards

@acer 

Your suggestion and correction is appreciated . I have view questions as I am new to some of the command used in the correction

1. What is the function of

 

kernelopts(version);

 

2. Is it possible to display the final resuts " all2 " one after the other as contained in the original code. i.e alpha[1],beta[0],beta[1],beta[2],gamma0,gamma1,gamma2 separately and on separate lines. Your suggestion will be highly appreciated.

Thank you and kind regards.

@Preben Alsholm 

Your suggestion and correction is appreciated

@vv 

Thank you. Your suggestion is appreciated.

@Kitonum 

Thank you. the correction you made is greatly appreciated.

1 2 3 Page 2 of 3