Inicio

n=0;incx=(b-a)/30;SE=0;SEC=0;SEC2=0


Hacer

punx1=a+n*incx
punx3=a+(n+1)*incx
punx2=a+incx/2+n*incx
SE=SE+(f(punx1)+4*f(punx2)+f(punx3))*incx/6
SEC=SEC+(punx1*f(punx1)+4*punx1*f(punx2)+punx1*f(punx3))*incx/6
SEC2=SEC2+(f(punx1)^2+4*f(punx2)^2+f(punx3)^2)*incx/6
SE=(SE=0)?0.000001:SE
xc=SEC/SE
yc=SEC2/(2*SE)
n=n+1

Mientras

n<30