Integrasi numerik metode trapesium menggunakan scilab

    a=0
    b=2
    n=256
    h=(b-a)/n
    x=a
    sigma=0
 in1= ((4*0)-(0^3))*exp(0^2)
 printf('%2d \t %11.6f \t %11.6f \n',1,a,in1);
 for i=1:n-1;
      x=x+h
      fr=((4*x)-(x^3))*exp(x^2)
        printf('%2d \t %11.6f \t %11.6f \n',i,x,fr);
        sigma=sigma+2*((4*x)-(x^3))*exp(x^2)
  end
  inakhir= ((4*2)-(2^3))*exp(2^2)
  printf('%2d \t %11.6f \t %11.6f \n',n,2,inakhir);
  jumlah=(in1+sigma+inakhir)/2
  disp("Hasil perhitungan trapesium=")
  printf('%11.6f',jumlah)