Budownictwo 2006 Wydział Inżynierii Lądowej Politechniki Krakowskiej

Maciek_K

program działa poprawnie i nic (z tego co się orientuje Smile

) mu nie brakuje. ...nawet rysuje te głupie prostokąciki
pozdrawiam mq

function y=f(x)
y=x ;
end

function s=calka(a,b,n)
h=(abs(b-a))/n;
ps=0;
x=(2*a+h)/2;
for i=1:n;
y=f(x);
ps=ps+(h*y);

ps
x

plot([x-(h/2),x-(h/2),x+(h/2),x+(h/2),],[0,y,y,0],"b")

x=x+h;

end
ps
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

a=0
b=10
n=100
x=a:n:b
y=f(x)

plot(x,y,"c")
hold on

s=calka(a,b,n)
pause()