Integrasi Numerik dengan menggunakan Scilab

• 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)

• Integrasi numerik metode simpson 1/3 menggunakan scilab

    a=0
    b=2
    n=256
    h=(b-a)/n
    x=a
    sigmaganjil=0
    sigmagenap=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);
        if modulo(i,2) ==0 then
        sigmagenap=sigmagenap+2*((4*x)-(x^3))*exp(x^2)
        else
           sigmaganjil=sigmaganjil+4*((4*x)-(x^3))*exp(x^2)  
        end
  end
  inakhir= ((4*2)-(2^3))*exp(2^2)
  printf('%2d \t %11.6f \t %11.6f \n',n,2,inakhir);
  jumlah=(in1+sigmaganjil+sigmagenap+inakhir)/3
  disp("Hasil perhitungan simpson=")
  printf('%11.6f',jumlah)

• Integrasi numerik metode titik tengah menggunakan scilab

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