📄 lurabe1.m
字号:
function [p,dp,q,r,alp,dom,x,wint,wipr]=urabe1(n)
dom=[0,2*pi];alpha=(dom(2)-dom(1))/2;beta=(dom(2)+dom(1))/2;
e1=ones(n,1)/2;e2=e1;e2(2)=1;X=spdiags([e1,e2],[-1,1],n,n);
X=alpha*X+beta*eye(n);
r=1;
alp='c';%the known eigenvalues and the eigenfunctions will be given in p,dp respectively
[x,wint]=pd(n,dom);
wipr=ones(n,1);
p=zeros(10,1);dp=zeros(n,10);
for k=1:10 p(k)=2*k-1;end;
for k=1:10 dp(:,k)=sin(p(k)*x);end;
for i=1:length(p)
p(i)=p(i)^2;
I(i)=wip(dp(:,i),dp(:,i),wipr,wint,dom);
dp(:,i)=dp(:,i)/sqrt(I(i));
end
q=0*dp;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -