Вот Ф-я:
Код
function Go
global hEd1 hEd2 hEd3 hEd4 hPop hAxes
str1=get(hEd1,'String');str2=get(hEd3,'String');
str3=get(hEd2,'String');str4=get(hEd4,'String');
x0=str2num(str1);y0=str2num(str2);
xf=str2num(str3);N=str2num(str4);
index= get(hPop,'Value');
cellArr=get(hPop,'String');
funName =cellArr{index};
dx=(xf-x0)/N;
X=x0;Y=y0;
for k=1:N
Y=[Y,Y(end) + feval(funName,X(end))*dx];
X=[X,X(end)+dx];
end
axes(hAxes);
plot(X,Y);
grid on;
global hEd1 hEd2 hEd3 hEd4 hPop hAxes
str1=get(hEd1,'String');str2=get(hEd3,'String');
str3=get(hEd2,'String');str4=get(hEd4,'String');
x0=str2num(str1);y0=str2num(str2);
xf=str2num(str3);N=str2num(str4);
index= get(hPop,'Value');
cellArr=get(hPop,'String');
funName =cellArr{index};
dx=(xf-x0)/N;
X=x0;Y=y0;
for k=1:N
Y=[Y,Y(end) + feval(funName,X(end))*dx];
X=[X,X(end)+dx];
end
axes(hAxes);
plot(X,Y);
grid on;
Может я конечно чего-то недопонимаю, но если кто-нить понял в чем тут ошибка плз поделитесь.
P.S. Проверить себя можно элементарно, хотя бы даже на sin(x)

Bot: код лучше соотв. образом отображать, а тему я переношу в этот раздел - более к месту будет.