Aire sous une courbe (2)

Soient \(a\) et \(b\) deux nombres réels (\(a < b\)) et \(f\) une fonction continue telle que l'on ait \(f(x)\ge0\) sur \(\left[a~;~b\right]\). On appelle \(C_f\) la représentation graphique de \(f\) dans un repère orthonormé.

On cherche à calculer une valeur approchée de l'aire du domaine situé entre :

  • l'axe des abscisses ;

  • la courbe \(C_f\) ;

  • les droites d'équations \(x=a\) et \(x=b\).

La figure ci-dessous (à gauche) illustre le problème avec :

  • \(f(x) = \frac12x^3 - \frac12 x^2 - \frac12 x + 1\) ;

  • \(a=-1\) ;

  • \(b = 2\).

Il est possible de calculer cette valeur approchée en utilisant la méthode des rectangles 1 suivante :

  • on découpe l'intervalle \(\left[a~;~b\right]\) en \(n>0\) subdivisions régulières. Chaque subdivision a pour largeur \(\text{d}x=\dfrac{b-a}n\) ;
  • il est possible de numéroter chaque subdivision avec un entier \(0 \leqslant k < n\) ;
  • dans la \(k\)-ième subdivision, on dessine le rectangle de largeur \(\text{d}x\) et de hauteur \(f\left(a+k\text{d}x\right)\) ;
  • la somme des aires de tous les rectangles donne une valeur approchée de l'aire du domaine étudié.

La figure de droite illustre cette méthode dans le cas \(n=16\). La valeur approchée obtenue vaut alors :

\[f(-1)\times\text{d}x+f(-1+\text{d}x)\times\text{d}x+f(-1+2\text{d}x)\times\text{d}x+\dots+f(-1+15\text{d}x)\times\text{d}x\simeq2,4888\]

Il est clair que plus le nombre de subdivisions est important, plus la valeur calculée est proche de la valeur exacte de l'aire cherchée.

Écrire la fonction aire qui prend en paramètres les entiers a et b ainsi que l'entier n strictement positif et renvoie la valeur approchée renvoyée par la méthode des rectangles décrite utilisant \(n\) rectangles.

assert ?

Le mot clé assert est utilisé en Python afin de vérifier que des propositions sont vraies.

Ainsi, l'expression assert 3 + 5*7 == 38 permet de vérifier que l'expression 3 + 5*7 est bien évaluée à 38.

Si c'est le cas, le programme continue de se dérouler normalement. Dans le cas contraire, le programme est interrompu et une erreur est signalée.

Comparaison de nombres flottants

Lorsqu'on écrit a = x ou x est un nombre réel, la valeur de a enregistrée en machine est une valeur approchée de x (quelques fois la valeur exacte). Cette valeur approchée a la forme d'un nombre flottant (le type float en Python). En conséquence, alors que des calculs et des comparaisons peuvent être effectués de manière exacte sur des réels, ils ne le sont que de manière approchée sur leur représentation en machine. On peut donc obtenir par exemple, avec a = x et b = y, l'expression a == b évaluée à True alors que x et y sont différents.

C'est pourquoi les tests ne vérifient pas l'égalité des résultats et des valeurs attendues mais leur proximité.

Ainsi, on peut vérifier que \(\sqrt{2} \approx 1,414214\) en faisantassert abs(1.414214 - sqrt(2)) < 1e-6. Ce test vérifie que les deux valeurs sont proches à \(10^{-6}\) près.

Exemples 
>>> a = -1
>>> b = 2
>>> aire(a, b, 1)
1.5
>>> aire(a, b, 2)
1.78125
>>> aire(a, b, 4)
2.1328125

###(Dés-)Active le code après la ligne # Tests (insensible à la casse)
(Ctrl+I)
Tronquer ou non le feedback dans les terminaux (sortie standard & stacktrace / relancer le code pour appliquer)
Si activé, le texte copié dans le terminal est joint sur une seule ligne avant d'être copié dans le presse-papier
Évaluations restantes : 5/5

.128013lS^]et-dA5f18umag,_/R=in{ 6)yàqPhcL\[(bsx.p;r4'90"ov+w73Ok:é *}2030c090a0k0r050I0%0C050k0I0I0q0S0a0r0L0S020u030I0i0j0j0k0N0x02060T050i0}0T0s0%000k0j0L0M0%0p09170N0z0i090I030o1416181a120L02031F1y1I0o1F120c0r0U0=0@0_0{0B0r0l0B051W0B0a10030-0H05091R0^0`0S1V1X1Z1X0a1(1*1%0a0N1G0a0B0=1d0I0L0k0s0{0)0S1,1T0S0f0/090s1l091%2224291.2c1*2f0j2h02040%0A0N0T0L0T0I0r1g1i0+200N0N090C2C1y2j0s1G0o1~2O1{1}1|1'0c2l0{1Z0s2e2z1%1O1Q0?1-2Y0r2!0s0T2'1%0L2H1G2M2O2^13231i2)2a2.0N1705100%0g2L2|112{2k2~1.3032340)3724392M2X0S3e0k33020%0Y3i2N123l3c0{3o3q0%0O3u3k2|3m3A340e3E3w3G3y3n0T313p340v3L3a2}1S3d3Q3f3r0X3V3x3Y3z3!3S3r0h3'3N3)3P3R3B0Q3/3b3;3I020g0R3_3X2*3=3#0g361z381J2?1y2'2R0c1}2W3O0C2/2r0*1P1G2=092@48473j034i0+4q3`420!100+0f3E0%3W3H0f104D4s2N4G3O0 020G3E4N3{100J4S3(424P0w0#3L0%4'4F4Y2a4A022H0a0i0N0s4E4T420j4I020R0K3K4L3r4@2a0T100'4?4*3d4V4X3:4253020'554 511.0j0r103t4 4)5b52100b565q5j4`4|4~2^5p4y5r5e5u5C58024W5h570{5d5f5a5G0{5k103h5o5i5M5s5F412a4_105y5Z3m5N5(3O0s595V5L0S5d0V5+3;5S3}3L0u3M5Q0S4,4C5@420s4`0.2H5P5!1.4P4R5K5v3z100k683m4P0m622 100H6i4O106l5/6e3n104=6d5~4!4$4 0u4(5B696f020I0T16096m1.5d0q6O5R5l4{4%4(5W6w020c5J5A6Y6Q6S0S6b6q4U026p6u5~5d5t6:6H6Z6h6z6^4!6)5d0o6)5-026y2^6E6X5:4,0f3Q71100!6~0W102,710H100N240l6N6{6j106c2`5:7274485:6B6W6F4'6Y726$386G5)106R6@3H6g6~105?7M5,7d7P5E7S6-6#7B7C7E106K6M7V7L6%7v7%6L2q7V7R7+6v724K7u6v6+7q7T5I6,4Z100w7V5g7=5~727Z6D6F6Y4,4.4:7x3j7I7}7'7/6D1y4v4p498o0o4c1y0a4e8t2U2P0k1)8q4c1E4x6^2H0j0n0f0k0!090n0B0Y101q1s1u1w0%6C2`1L391F0D3Q0I8T8h3;182B0B170a090J100E0G0s0E825o8.0s2W0U090N0;6Y8*1~8-8/8;0G0V0E2,0f0a0x8^4E0@0%8j1+0C3p0C0i050$1+8{8}8 0;050P661+0c0i0%0c6L0.2!0%0$4/1O0$0K0%0Z1i2.8.0=0,9R9y0%0y0%9w0.2E0c0P0i9G131{1h0l022,0a0$0l7m1*2L9*0s9,8V1J392'3m938,0k8.8:020E0F0d0q9a0s0a0n0k070H0%0f0G0J0w0Ea30a0t0c0(0J0E084b4m8D0D0k0%0U3p094:0%8/0k2J9z9k9T1+9.9:9=1+1w0a8(42a095a48=0)0m0v0)0e9f4 0K9|120B0N0%0o8Z02ay0%0N9:9o930$9C1+9h140H2'0r2A1h0%0L8 2qaS2E2A1m1Z2c8 9Y9k9m9oaF1i8J2J3Qb224aS23a-1*9v187o4:0;6C1M4l2(8)0N8+aW8;0F0E1*9c0Gai0k0w0VbO0V0gam8/aoaqakbR0G0k0V0)bV0JbXarbQ0E9D0a0Ib!b$8?0b0galanapb+0EbF1?b_0r2q0Ib(b*asbb090zaZ0O0hccatav0U9}1P3m1:1Y1!1$8E3m2n2e2g102t060CbF0LaS0A0x1~1h4S4o8E2_488nco3O60097^7H7E650r677|3;7{7_867OcU80026t856^726/cX6|6s7c737 2a7A89777?7-7(7X5c7K6)5_3 c?7D7,6!7G8g6'c}c{c;7sc:5Hc*cP5:6=c.6`c+7r028_c'7J0270d95H8f3v8a784J7bds6I7edA5;7g9-du507v7k4-7n7pdk6r4Qdc6IdH6Yc=757Cc@cY7~dD6(dD72djdf6v5=c.dCdo3O5*d%4Bd5dvdXaT6n6J7.dNd)6;d8d.6-9j7:c.cO4t7zdbc!d`d?cKcV814Ed_1.0C0g1001ahaj0w9iaD05aFbr0^enak0%0'9Cec76d26v8c0,8ec.e3d#54c.885Ad^eh0{ejel0%b40j9o0}0L1Za}00010Y0Q0M2e0I1r8Aa|9YazeJ758m4j2OcFcf398r0,0.0:02.

###(Dés-)Active le code après la ligne # Tests (insensible à la casse)
(Ctrl+I)
Tronquer ou non le feedback dans les terminaux (sortie standard & stacktrace / relancer le code pour appliquer)
Si activé, le texte copié dans le terminal est joint sur une seule ligne avant d'être copié dans le presse-papier
Évaluations restantes : 5/5

.128013lS^]et-dA5f18umag,_/R=in{ 6)yàqPhcL\[(bsx.p;r4'90"ov+w73Ok:é *}2030c090a0k0r050I0%0C050k0I0I0q0S0a0r0L0S020u030I0i0j0j0k0N0x02060T050i0}0T0s0%000k0j0L0M0%0p09170N0z0i090I030o1416181a120L02031F1y1I0o1F120c0r0U0=0@0_0{0B0r0l0B051W0B0a10030-0H05091R0^0`0S1V1X1Z1X0a1(1*1%0a0N1G0a0B0=1d0I0L0k0s0{0)0S1,1T0S0f0/090s1l091%2224291.2c1*2f0j2h02040%0A0N0T0L0T0I0r1g1i0+200N0N090C2C1y2j0s1G0o1~2O1{1}1|1'0c2l0{1Z0s2e2z1%1O1Q0?1-2Y0r2!0s0T2'1%0L2H1G2M2O2^13231i2)2a2.0N1705100%0g2L2|112{2k2~1.3032340)3724392M2X0S3e0k33020%0Y3i2N123l3c0{3o3q0%0O3u3k2|3m3A340e3E3w3G3y3n0T313p340v3L3a2}1S3d3Q3f3r0X3V3x3Y3z3!3S3r0h3'3N3)3P3R3B0Q3/3b3;3I020g0R3_3X2*3=3#0g361z381J2?1y2'2R0c1}2W3O0C2/2r0*1P1G2=092@48473j034i0+4q3`420!100+0f3E0%3W3H0f104D4s2N4G3O0 020G3E4N3{100J4S3(424P0w0#3L0%4'4F4Y2a4A022H0a0i0N0s4E4T420j4I020R0K3K4L3r4@2a0T100'4?4*3d4V4X3:4253020'554 511.0j0r103t4 4)5b52100b565q5j4`4|4~2^5p4y5r5e5u5C58024W5h570{5d5f5a5G0{5k103h5o5i5M5s5F412a4_105y5Z3m5N5(3O0s595V5L0S5d0V5+3;5S3}3L0u3M5Q0S4,4C5@420s4`0.2H5P5!1.4P4R5K5v3z100k683m4P0m622 100H6i4O106l5/6e3n104=6d5~4!4$4 0u4(5B696f020I0T16096m1.5d0q6O5R5l4{4%4(5W6w020c5J5A6Y6Q6S0S6b6q4U026p6u5~5d5t6:6H6Z6h6z6^4!6)5d0o6)5-026y2^6E6X5:4,0f3Q71100!6~0W102,710H100N240l6N6{6j106c2`5:7274485:6B6W6F4'6Y726$386G5)106R6@3H6g6~105?7M5,7d7P5E7S6-6#7B7C7E106K6M7V7L6%7v7%6L2q7V7R7+6v724K7u6v6+7q7T5I6,4Z100w7V5g7=5~727Z6D6F6Y4,4.4:7x3j7I7}7'7/6D1y4v4p498o0o4c1y0a4e8t2U2P0k1)8q4c1E4x6^2H0j0n0f0k0!090n0B0Y101q1s1u1w0%6C2`1L391F0D3Q0I8T8h3;182B0B170a090J100E0G0s0E825o8.0s2W0U090N0;6Y8*1~8-8/8;0G0V0E2,0f0a0x8^4E0@0%8j1+0C3p0C0i050$1+8{8}8 0;050P661+0c0i0%0c6L0.2!0%0$4/1O0$0K0%0Z1i2.8.0=0,9R9y0%0y0%9w0.2E0c0P0i9G131{1h0l022,0a0$0l7m1*2L9*0s9,8V1J392'3m938,0k8.8:020E0F0d0q9a0s0a0n0k070H0%0f0G0J0w0Ea30a0t0c0(0J0E084b4m8D0D0k0%0U3p094:0%8/0k2J9z9k9T1+9.9:9=1+1w0a8(42a095a48=0)0m0v0)0e9f4 0K9|120B0N0%0o8Z02ay0%0N9:9o930$9C1+9h140H2'0r2A1h0%0L8 2qaS2E2A1m1Z2c8 9Y9k9m9oaF1i8J2J3Qb224aS23a-1*9v187o4:0;6C1M4l2(8)0N8+aW8;0F0E1*9c0Gai0k0w0VbO0V0gam8/aoaqakbR0G0k0V0)bV0JbXarbQ0E9D0a0Ib!b$8?0b0galanapb+0EbF1?b_0r2q0Ib(b*asbb090zaZ0O0hccatav0U9}1P3m1:1Y1!1$8E3m2n2e2g102t060CbF0LaS0A0x1~1h4S4o8E2_488nco3O60097^7H7E650r677|3;7{7_867OcU80026t856^726/cX6|6s7c737 2a7A89777?7-7(7X5c7K6)5_3 c?7D7,6!7G8g6'c}c{c;7sc:5Hc*cP5:6=c.6`c+7r028_c'7J0270d95H8f3v8a784J7bds6I7edA5;7g9-du507v7k4-7n7pdk6r4Qdc6IdH6Yc=757Cc@cY7~dD6(dD72djdf6v5=c.dCdo3O5*d%4Bd5dvdXaT6n6J7.dNd)6;d8d.6-9j7:c.cO4t7zdbc!d`d?cKcV814Ed_1.0C0g1001ahaj0w9iaD05aFbr0^enak0%0'9Cec76d26v8c0,8ec.e3d#54c.885Ad^eh0{ejel0%b40j9o0}0L1Za}00010Y0Q0M2e0I1r8Aa|9YazeJ758m4j2OcFcf398r0,0.0:02.

  1. On propose ici une version simplifiée de cette méthode qui considère des rectangles ayant tous la même largeur. ↩