1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
|
Discontinuity in solution
0.00000000e+00 1.00000000e+00
2.34520788e-02 9.77085318e-01
4.69041576e-02 9.54695719e-01
1.26678715e-01 8.82381196e-01
2.06453273e-01 8.14935604e-01
2.86227830e-01 7.52436861e-01
4.50967660e-01 6.37579007e-01
6.15707490e-01 5.39988032e-01
7.80447320e-01 4.57240552e-01
9.45187150e-01 3.87139827e-01
1.00000000e+00 3.66442563e-01
1.00000000e+00 1.00000000e+00
1.02345208e+00 9.77085318e-01
1.04690416e+00 9.54695719e-01
1.12667872e+00 8.82381196e-01
1.20645327e+00 8.14935604e-01
1.28622783e+00 7.52436861e-01
1.45096766e+00 6.37579007e-01
1.61570749e+00 5.39988032e-01
1.78044732e+00 4.57240552e-01
1.94518715e+00 3.87139827e-01
2.00000000e+00 3.66442563e-01
Number of steps: 10 + 10 = 20
Discontinuity in RHS: Case 1 - explicit treatment
0.00000000e+00 1.00000000e+00
2.34520788e-02 9.77085318e-01
4.69041576e-02 9.54695719e-01
1.26678715e-01 8.82381196e-01
2.06453273e-01 8.14935604e-01
2.86227830e-01 7.52436861e-01
4.50967660e-01 6.37579007e-01
6.15707490e-01 5.39988032e-01
7.80447320e-01 4.57240552e-01
9.45187150e-01 3.87139827e-01
1.00000000e+00 3.66442563e-01
1.00000000e+00 3.66442563e-01
1.00504558e+00 3.57425469e-01
1.01009116e+00 3.48630259e-01
1.02686958e+00 3.20932037e-01
1.04364801e+00 2.95189893e-01
1.06042643e+00 2.71428693e-01
1.09875325e+00 2.23778534e-01
1.13708008e+00 1.84344072e-01
1.15795394e+00 1.65976090e-01
1.17882780e+00 1.49454458e-01
1.19970166e+00 1.34583024e-01
1.22057552e+00 1.21193305e-01
1.24144938e+00 1.09136406e-01
1.27686095e+00 9.13557908e-02
1.31227252e+00 7.65039734e-02
1.34768410e+00 6.40955034e-02
1.42252133e+00 4.42507080e-02
1.46951105e+00 3.50520256e-02
1.51650078e+00 2.77477961e-02
1.56349051e+00 2.19581210e-02
1.61048024e+00 1.73751345e-02
1.65746997e+00 1.37493280e-02
1.73029436e+00 9.57464470e-03
1.80311876e+00 6.66725545e-03
1.87594316e+00 4.63257604e-03
1.94876756e+00 3.21043602e-03
2.00000000e+00 2.47989685e-03
Number of steps: 10 + 26 = 36
Discontinuity in RHS: Case 2 - let CVODE deal with it
0.00000000e+00 1.00000000e+00
2.34520788e-02 9.77085318e-01
4.69041576e-02 9.54695719e-01
1.26678715e-01 8.82381196e-01
2.06453273e-01 8.14935604e-01
2.86227830e-01 7.52436861e-01
4.50967660e-01 6.37579007e-01
6.15707490e-01 5.39988032e-01
7.80447320e-01 4.57240552e-01
9.45187150e-01 3.87139827e-01
1.00000000e+00 3.66442563e-01
1.00000000e+00 3.66442563e-01
1.00220332e+00 3.63503469e-01
1.00440663e+00 3.59880618e-01
1.00660995e+00 3.56060814e-01
1.00881326e+00 3.52201766e-01
1.01101658e+00 3.48357263e-01
1.01810972e+00 3.36425359e-01
1.02520285e+00 3.24902144e-01
1.04291574e+00 2.97727832e-01
1.06062864e+00 2.72595694e-01
1.07834153e+00 2.49486977e-01
1.12193062e+00 2.00239352e-01
1.14404862e+00 1.79170304e-01
1.16616663e+00 1.60310191e-01
1.18828463e+00 1.43435282e-01
1.21040264e+00 1.28336161e-01
1.23252065e+00 1.14826356e-01
1.26887350e+00 9.56514614e-02
1.30522635e+00 7.97064617e-02
1.34157919e+00 6.64607835e-02
1.37793204e+00 5.54387340e-02
1.45323392e+00 3.81470272e-02
1.52853579e+00 2.62731987e-02
1.60383767e+00 1.81158684e-02
1.67913954e+00 1.24124410e-02
1.75444142e+00 8.47041328e-03
1.82974329e+00 5.73179862e-03
1.90504516e+00 3.89153977e-03
1.98034704e+00 2.63536234e-03
2.00000000e+00 2.37913555e-03
Number of steps: 10 + 29 = 39
|
