Доброго времени суток.
При решении в maple стационарного уравнения Шредингера возникла проблема нахождения уровней энергии. А именно не получается оперировать численно найденной пси-функцией и ее производной. Необходимо разделить
значение производной в конкретной точке на значение функции в этой же точке. Но насколько понял, машина пытается делить сами функции друг на друга, а потом считать вычислять значение получившегося в точке.
Код:
L := dsolve({conL, shred}, y(x), numeric);
R := dsolve({conR, shred}, y(x), numeric);
fl := proc (x) options operator, arrow; L(x)[2] end proc;
gl := proc (x) options operator, arrow; L(x)[3] end proc;
fr := proc (x) options operator, arrow; R(x)[2] end proc;
gr := proc (x) options operator, arrow; R(x)[3] end proc;
al := subs(x = 6, fl(x));
al := %;
subs(x = 6, gl(x));
bl := %;
subs(x = 6, fr(x));
ar := %;
subs(x = 6, gr(x));
br := %;
raz := br/ar-bl/al
Пробовал присваивать значения переменным al, bl, ar, br по разному, так же пробовал задать raz как функцию от точки сшивки, с тем же успехом. Как можно обойти эту проблему?