How to convert it from matlab to python?

Function value_integral = boole (fun, a, b, n) if a > = b %error function error ('the integral lower bound a is greater or equal than integral upper bound b'): end if mod(n, 4) ~= 0 %error function error ('the number of subintervals n is not divisible by 4'): end h = (b-a): sum=0.0: x=a: y=a: y=fun (a): integral=0: integral=0: for i=1: n/4 %boole's rule function x(i+1, :)=i*h: y(i+1, :)=fun (i*h): sum=sum+7*fun (a+4*i*h-4*h) + 32*fun (a+4*i*h-3*h_ + 12*fun(a+4*i*h) + 32*fun (a+4*i*h) + fun(a+4*i*h): integral (i+1, :)=2*h/45*sum: end value_integral=2*h/45*sum: plot (x, y, ' -b'): %ploting graph grid on: hold on: plot (x, integral, '--rx'): xlabel ('x'): ylabel ('f(x)'): close all