octave и lambda-функции

А можно чуть побольше код .... проверю на свежей версии.

И да всегда можно написать в рассылку

>> basis=@(x,k) quad(@(t) polyval(legendrepoly(k),t),-1,x);
>> quad(@(y) basis(y,k),-1,1);
error: 'k' undefined near line 1 column 19
error: evaluating argument list element number 2
error: called from:
error:    at line -1, column -1
error: quad: evaluation of user-supplied function failed
>>

GNU Octave, version 3.7.6+

удалось немного локализовать. Вот этот код вызывает уже ошибку.

function retval=basis(x,k)
       retval=quad(@(t) sin(k*t),-1,1);
end

quad(@(y) basis(y,2),-1,1);

quad(@(y) basis(y,2),-1,1);

//fixed.

так надо. я просто опечатался в топике.

пользуясь случаем, реквестую статей, чтобы научиться строить графики с наложением

А почему в толксы, а не в Development?

Скорее всего потому как тут народу больше шастает

А мне за такую фигну пару лет назад шкворец порезали

«О времена, о нравы!» (с)

hold on и hold off должны помочь, если я тебя правильно понял.

Почитай мануалы по gnuplot. Один день убьешь, зато наделаешь себе шаблонов для дальнейшего использования.

P.S. Пользуясь случаем, реквестирую статей, чтобы понять, что за непонятную штуку написал автор.

basis=@(x,k) quad(@(t) polyval(legendrepoly(k),t),-1,x);

определяю функцию

S_k(x)=\int_{-1}^x P_k(t) dt, где P_k(t) - это полином Лежандра степени k.

quad(@(y) basis(y,2),-1,1);

здесь я пытаюсь взять интеграл от функции S_2(x):

\int_{-1}^1 S_2(x) dx

и на этом моменте октава фейлится.

А зачем ты дважды интеграл вычисляешь? Я, конечно, с этими @ никогда не работал, но сдается мне, что у тебя получается, что basis — это интеграл от -1 до x по полиному Лежандра степени k. Т.е. тебе разве не basis(2, 1) надо во втором случае писать?

да. два раза вычисляю. Если я напишу

basis(2, 1)

то это будет просто

S_2(x)=\int_{-1}^1 P_k(t) dt

а мне нужен именно набор функций

S_k(x), который будет является базисом в пространстве W_2^1[-1,1]