📄 untitled2.m
字号:
clc;close all;
tic;
a=0:0.02:0.3;
ts=linspace(0,12000,200);
p0=0.05*0.95.^[0:28];
for k=1:29;eval(['P',num2str(k),'=[];']);end;
for k=1:length(a);
[t,p]=ode23('haha',ts,p0,[],a(k));
for h=1:29;
E=p(:,h);
eval(['P',num2str(h),'=[P',num2str(h),',E];']);
end
end
figure;mesh(a,ts,P2);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P3);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P4);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P5);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
t=toc
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -