Question

怎么用MATLAB求解一维周期性势场中电子的能带并画图

只要在MATLAB中画图的编程就行啦，谁帮帮忙啊，谢谢了
这是一道量子力学的问题，求在MATLAB中的编码

Answer 1

%我想用matlab对函数dy/dx=（x-y）/（1-y-x）进行模拟，请问该怎样做，
%解微分方程！！
%归一化：
%令y=y（1）
%x=t=y(2)
%dy(1)/dt=(y(2)-y(1))/(1-y(1)-y(2))
%dy(2)/dt=1;
%函数文件
founction dy=fun1(t,y)
dy=zeros(2,1);
dy=[(y(2)-y(1))/(1-y(1)-y(2));1];
%以上保存为fun1.m文件
%以下是脚本程序
clear
ts=0:0.01:?%时间范围
y0=[?;?];%函数边界条件y(0),x(0)
[t,y]=ode45('fun1',ts,y0);
plot(y(:,1),y(:,2))