区别于WhittakerM功能和替换

Question

我试图区分WhittakerM功能。 为了解决WhittakerM方程，我们有：

dsolve( 'D2y+(-1/4+Landa/r+(1/4-(L+1/2)^2)/r^2)*y=0' ,'r')

C1*WhittakerM(Landa,L+1/2,r)+C2*WhittakerW(Landa,L+1/2,r)

从边界条件我只需要WhittakerM(Landa,L+1/2,r)/r为问题的条件添加1/r。我尝试对其进行区分，然后在某些方面进行替换，但在subs和diff中存在一些错误。

Landa=1;L=0; % # for simplicity
R1=inline('WhittakerM(Landa,L+1/2,r)/r','r');
Rp1=diff(R1,r);
r=1:0.01:20;
R1sub=eval(R1,r);
Rp1sub=eval(Rp1,r);

你知道吗？

Answer 1

无论上面的错误如何（说实话，我并不完全理解你试图通过替换实现的目标），使用符号工具箱可能是一个好的开始：

syms r
%defines a symbolic function
R1(r)=whittakerM(landa,l+1/2,r)/r
%differentiate
Rp1=diff(R1,r);
%evaluate
Rp1_e=Rp1(1:0.01:20)