LINUX.ORG.RU

процедура в maple.


0

2

А если серьёзно, то бери любую удобную систему компьютерной алгебры для такой задачи:

gramschmidt:=proc(f)
F:=[(t->f[1](t)/sqrt(scalprod(f[1],f[1])))];
for i from 2 to N do
       F:=[op(F),     t->(f[i](t)-add(scalprod(f[i],f[j])*f[j](t),j=1..i-1))*1/sqrt(scalprod(f[i](t)-add(scalprod(f[i],f[j])*f[j](t),j=1..i-1),f[i](t)-add(scalprod(f[i],f[j])*f[j](t),j=1..i-1)))];
end do:return F;end proc:
★★☆☆☆

Последнее исправление: dikiy (всего исправлений: 1)
Ответ на: комментарий от stevejobs

это мои потуги реализовать то, что я описал в прошлом треде на maple. При это f инициализируется так:

N:=7; unassign('f'); f:=[];
for j from 1 to N do
f:=[op(f), unapply(orthopoly[H](j,t),t)];
end do:
dikiy ★★☆☆☆
() автор топика
Последнее исправление: dikiy (всего исправлений: 1)
Вы не можете добавлять комментарии в эту тему. Тема перемещена в архив.