... 28f=++++∫Đểtăngđộchínhxáctadùngphươngphápchiađoạn[a,b]thànhmđoạn và trênmỗiđoạntadùngcôngthứcHardy.Taxâydựnghàminthardy()đểthựchiệncôngthứctrên:functiontp=inthardy(f,a,b,m) 311CHƯƠNG 6: ĐẠO HÀM VÀ TÍCH PHÂN SỐ ... (10)Chiahaivếcủa(10)cho‐3tanhậnđược:⋅⋅⋅+++′=−=6644ha165ha41)x(f4)1,1(D)1,2(D4)2,2(D (11)TrongkhiD(1,1) và D(2,1)saikhácf′(x)phụthuộcvàoh2thìD(2,2)saikhácf′(x)phụthuộcvàoh4.Bâygiờtalạichiađôibướch và nhậnđược: 330%TinhtichphanbangphuongphapHardyh=(b‐a)/(6*m);fork=1:6*mx(k)=a+k*h;endtp=0;j=1;fork=1:mtp=tp+(28*feval(f,a)+162*feval(f,x(j))+ ... cácđoạnhàmbiếnđổinhiều và saisốnhỏtrêncácđoạnhàmtươngđốibằngphẳng.Ngượclạiphươngphápcầuphươngthíchnghichiacácđoạnkhôngđều:ngắntrêncácđoạnhàmthayđổinhiều và dàitrêncácđoạnthayđổiít và sẽcó saisốnhỏkhisốđoạnchianhỏ. Thuậttoáncầuphươngthíchnghibắtđầubằngviệctínhtíchphânintđốiv ớitoànbộđoạn[a,b] và tổngtíchphânint12=int1+int2trên2đoạnbằngnhau.Dựatrênint và int12tatínhsaisố.Nếuchưađạtđộchínhxác,tachiađôimỗiđoạn và lặplạiquátrìnhtính.Tadùnghàmadaptivesimpson()đểthựchiệnthuậttoánnày:functionint=adaptivesimpson(f,a,b,tol)mid=(b+a)/2.0;int=simpsonapprox(f,a,b);int12=simpsonapprox(f,a,mid)+simpsonapprox(f,mid,b);if(abs(int‐int12)<15.0*tol)int=int12;elseleftint=adaptivesimpson(f,a,mid,tol/2);rightint=adaptivesimpson(f,mid,b,tol/2);...