đây là code tính tích phân hàm e^(-x) * sin(x)
Code:
//tinh tich phan bang phuong phap hinh_thang;
#include <conio.h>
#include <stdio.h>
#include <math.h>
float f(float x)
{
float a=exp(-x)*sin(x);
return(a);
};
void main()
{
int i,n;
float a,b,x,y,h,s,tp;
printf("Tinh tich phan theo phuong phap hinh thang\n");
printf("Cho can duoi a = ");
scanf("%f",&a);
printf("Cho can tren b = ");
scanf("%f",&b);
printf("Cho so buoc n = ");
scanf("%d",&n);
h=(b-a)/n;
x=a;
s=(f(a)+f(b))/2;
for (i=1;i<=n;i++)
{
x=x+h;
s=s+f(x);
}
tp=s*h;
printf("Gia tri cua tich phan la : %10.6f\n",tp);
getch();
}
sử dụng phương pháp hình thang để tính
VD như 1 hàm e^(-x) * sin(x)
ta sẽ chia trên đồ thị thành các hình thang cong với độ dài hoành độ nhất định ( ta sẽ chọn ngẫu nhiên nhỏ hơn độ dài cận lớn - cận nhỏ )