sympson 1/3 yöntemi ile integral hesaplama

%sympson 1/3 kuralına göre integral heasaplama 
%sympson 1/3 kuralı n nin tek olduğu durumlarda çalışır
%f fonkisyounun kodunu ayrı bir m-file da açmak gerek


a=input('alt sını değerini girinz:');

b=input('üst sınır degerini giriniz:');

n=input('segment sayısını ggirniz:');//aralık sayısı

h=(b-a)/n;

if rem(n,2)==0;//mod görevi görüyor

else

disp('verilen n degeri gecersiz' )

n=input('yeni n degeri');

end



k=0;%0a atamayınca hata veriyor

p=0;



x=a:h:b;

for i=2:n

if rem(i,2)==0 %indisi çift olanlar(xin indisi 1den başladığı için)

k=k+4*f(x(i));

else %indisi tek olanalr

p=p+2*f(x(i));

end

end

I=(f(a)+f(b)+k+p)*(h/3)

syms x;

trueval=int(sin(x),a,b)

Et=trueval-I

Ea=-((((b-a)^5)/(180*n^4))*(int(sin(x),a,b)/(b-a)))
F FONKSİYONU
function y=f(x);
y=sin(x);
OUTPUT
alt sını değerini girinz:0

üst sınır degerini giriniz:pi
segment sayısını ggirniz:4
I =2.0046
trueval =2
Et =
-10267655424369/2251799813685248
Ea =
-7656611879816833/(576460752303423488*pi)