Я думаю, что вы должны построить график, используя plot(x,int_f);
таким образом, чтобы вы построили два массива, а не один массив против одного единственного числа.
Также ваш интегральный шаг неверен, он должен быть int_f(n) = f(n)*dx+int_f(n-1)
, за исключениемдля первого запуска, где оно должно быть int_f(n) = f(n)*dx