< M A T L A B > Copyright 1984-2005 The MathWorks, Inc. Version 7.1.0.183 (R14) Service Pack 3 August 02, 2005 To get started, select MATLAB Help or Demos from the Help menu. >> % formule di integrazione dei trapezi e di Cavalieri Simpson >> % verifica con funzione lineare: devono dare risultati esatti entrambe >> f=@(x) 1-3*x f = @(x) 1-3*x >> integrale(-1,1,5,f) ans = 2.0000 >> integrale(-1,1,6,f) ans = 2 >> cavsimp(-1,1,3,f) ans = 2.0000 >> % verifica con funzione quadratica: Cavalieri Simpson deve essere esatta >> f=@(x) x*x f = @(x) x*x >> cavsimp(-1,1,3,f) ans = 0.6667 >> format long >> cavsimp(-1,1,3,f) ans = 0.66666666666667 >> integrale(-1,1,6,f) ans = 0.70370370370370 >> integrale(-1,1,100,f) ans = 0.66680000000000 >> integrale(-1,1,500,f) ans = 0.66667200000000 >> integrale(-1,1,1000,f) ans = 0.66666800000000 >> integrale(-1,1,10000,f) ans = 0.66666667999992 >> % sembrano funzionare >> % prova con la funzione 1/x >> f=@(x) 1/x f = @(x) 1/x >> log(2) ans = 0.69314718055995 >> integrale(1,2,100,f) ans = 0.69315343048182 >> integrale(1,2,1000,f) ans = 0.69314724305996 >> integrale(1,2,10000,f) ans = 0.69314718118497 >> integrale(1,2,100000,f) ans = 0.69314718056493 >> cavsimp(1,2,100,f) ans = 0.69314718057948 >> cavsimp(1,2,1000,f) ans = 0.69314718055997 >> log(2) ans = 0.69314718055995 >> % sembra soddisfacente