e00f97cd-9197-4de2-9ec5-d902d7556160 12 4 5 false false 0 0 decimal

&[DATE] &[TIME] - &[FILENAME]

One variable

x φ x 2^x 2^cos * 2 x ln * - 3 + :

a b 12 mat 1712 mat :

N 1000 :

0 a b N 1 - t u D t φ : Rkadapt N 2 el 1.1494 -

x φ x a b int maple 1.1494 -

Two variables - Rectangle domain

x y f 3 x 2^* 3 y 2^* + :

a b 12 mat 1412 mat :

Nx 200 :

Ny 200 :

c d 12 mat 1 - 612 mat :

x u Iy t v D x t f : 0 c d Ny 1 - D Rkadapt Ny 2 el 21 line :

0 a b Nx 1 - Iy Rkadapt Nx 2 el 72

Wrong

0 a b Nx 1 - Iy dn_ExplicitRK45 Nx 2 el 1092

x y f y c d int x a b int maple 1092

Two variables - Double Interal type I

x y f x y * y 2^+ :

a b 12 mat 2 - 312 mat :

Nx 200 :

Ny 200 :

x c 0.02 x 3^* 1 - :

x d 2 0.1 x 2^* - :

x u Iy t v D x t f : 0 x c x d Ny 1 - D Rkadapt Ny 2 el 21 line :

x c x a b lele * x d x a b lele * 2 1 sys

0 a b Nx 1 - Iy Rkadapt Nx 2 el 6.6667

Wrong

0 a b Nx 1 - Iy dn_ExplicitRK45 Nx 2 el 13.1114

x y f y x c x d int x a b int maple 13.1114

OK

x Iy x y f y x c x d int maple : 12520 x 2^- (3^* 50 - x 3^+ (3^- 75 x * 2520 x 2^- (2^* 50 - x 3^+ (2^- (* + 375000 /

x Iy x a b int 13.1114

The maple's Iy(x) and Rkadapt's Iy(x,0) are more or less the same

x Iy x a b lele * x 0 Iy 0.1 + ( x a b lele * 2 1 sys

Alvaro