سؤال حول matlab

[khawla MI]

السلام عليكم

سؤالي بسيط لكن لم أفلح في حله وهو أنه لدي برنامج matlab تنفيذه يستغرق مدة ساعة والسبب يعود الى أني استخدمت deux boucle for حاولت استبدالهما par faire ce que l'on appelle la vectorisation لكني لم أتحصل على النتيجة المرجوة

هل بالامكان مساعدتي في أيجاد طريقة لاستبدال les deux boucle

function Ing=integrale(f,a,b,c,d)
%l'approximation de l'integrale double par la methode de rectangle gauche
N = b*10-a;
M=d*10-c;
hx = (b - a)/N;
hy=(d - c)/M;

y = [c:hy:d];
x=[a:hy:b];
Ing = 0;
for i = 1:N
for j=1:M
Ing = Ing + hx*hy*f(x(i),y(j));
end
end
le probleme se pose lorsque le b et le d qui sont les bords de l'integrale tend vers un grand nombre
f dans ce programme est une fonction à deux variable definie comme suit
function r=f(x,y).
r=x+y;%par exemple
merci