과제중 Personal

format long
syms x
f = inline('x^-1'); %적분하고자 하는 함수
g = x^-1;
a = 1; %적분 시작 위치
b = 2; %적분 끝 위치
n = 4; %적분 구간 내의 지점 수
h = (b-a)/n;
s = feval(f,a);
for (i=1:n-1)
    x(i)=a+h*i;
    s=s+2*feval(f,x(i));
end
s = s+feval(f,b);
i = h*s/2;
r = int(g,'x',a,b);
error = abs(i-r);
double(i)
double(r)
double(error)

덧글

  • YsHead 2013/04/29 21:53 # 삭제 답글

    trapezoidal rule 이네요. ㅎㅎ
    1년만에, 또 우연찮게 이번엔 다비치 검색하다가, 들어오게 됬네요.
    일본에서 찍은 사진은 블로그에만 모셔 놓지 마시고,
    카카오 스토리에 올리면 사람들이 형한테 호감 많이 가지실것 같은데요?
    여행 다닌 사진들이 너무 아까워요, 자랑좀 하시지 ㅎㅎ
    아주 멋진 옷을 집에서만 입는 느낌이랄까.

    그리고 다른 관심사도 염탐하다 갑니다. 조금 충격이긴 하지만,
    똘레랑스 ㅎㅎ
댓글 입력 영역


유용원의 군사세계

다음위젯2