블로그 이미지
Terapi

calendar

      1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30    

'c로쓴수치해석적분'에 해당되는 글 1

  1. 2008.12.11 c로 쓴 수치해석 적분
2008.12.11 10:25 System/System 이론

f(x)는 적분구하고자하는 함수입니다.

 

#include < stdio.h >
#include < math.h >

double f(double x)
{
    return(pow(x,3));
}

int main()
{
    int i, n;
    double a,b;
    double h;
    double s0, s1, s2, s;
    double x;

    printf("Input a= ");
    scanf("%lf", &a);
   
    printf("Input b= ");
    scanf("%lf", &b);

    printf("Input n= ");
    scanf("%d", &n);

    h = (b-a)/(2.0*n);

    printf("h= %lf\n",h);

    s0 = f(a) + f(b);
    s1 = 0;
    s2 = 0;

    printf("\n");

    for(i = 1; i <= 2*n-1; i++)
    {
        x = a+h*i;
        if(i % 2 == 0)
        {
            s2 = s2 + f(x);
        }
        else
        {
            s1 = s1 + f(x);
        }
    }
    s = h*(s0+4.0*s1+2*s2)/3.0;

    printf("y0 = %6.4lf\n", f(a));
    for(i=1;i<2*n-1;i++)
    {
        x=a+h*i;
        if(i%2 == 0) printf("y%d = %6.4lf\n", i/2, f(x));
    }
    printf("y%d = %6.4f\n", n, f(b));

    printf("\n /%4.2lf", b);
    printf("\n l    f(x)dx = %lf", s);
    printf("\n /%4.2lf", a);
    printf("\n");

 return 0;
}

신고

'System > System 이론' 카테고리의 다른 글

usb 모듈...  (0) 2009.01.02
c로 쓴 수치해석 적분  (0) 2008.12.11
고수준 접근과 저수준 접근  (0) 2008.12.11
posted by Terapi
prev 1 next

티스토리 툴바