< 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. >> % prova di calcolo del logaritmo mediante integrazione di 1/x >> % estremi di integrazione a=1 e b=2 >> a=1 a = 1 >> b=2 b = 2 >> % funzione da integrare f >> f=@(x) 1./x f = @(x) 1./x >> % precisione doppia >> format long >> % risultato aspettato >> log2=log(2) log2 = 0.69314718055995 >> % integrazione con la funzione interna di matlab >> quad(f,a,b) ans = 0.69314719986297 >> % integrazione col metodo dei trapezi implementato da noi >> 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 >> % esercizio: plottare l'errore relativo al variare di n