Решение обыкновенного диференциального уравнения.

[romchic]

Помогите, пожалуйста.У меня такая проблема.
мне дано уравнение xy'-y/(lnx)=0, y(e)=1.
Сначала я пишу три m-файла.
Например м-файл для метода Эйлера:

 
function [x,y]=eiler(x1,y1,n,b)
x(1)=x1
y(1)=y1
h=(b-x1)/n
for i=1:n
    x(i+1)=x(i)+h
    y(i+1)=y(i)+h*func(x(i),y(i))
end

Далее в Command Window решаю диф. уравнение:
n=10
x0=pi=3.1416
y0=1
b=x0+1
b=4.1416
и вызываю м-файл с помощью [x,y]=eiler(x0,y0,n,b)

Он начинает решать и выдаёт ошибку, что в строчке y(i+1)=y(i)+h*func(x(i),y(i)) есть какая-то ошибка или вообще не решает ничего и пишет "??? Undefined function or method 'func' for input
arguments of type 'double'."
Кто знает, в чём проблема?
Заранее спасибо!

[romchic]

PC Программа в MatLab

[p51x]

Переведите, програма ясно все написала