egode45.m

来自「内燃机转子仿真」· M 代码 · 共 9 行

clear all;
clc;
close all;
format long;
%[T,y]=ode45(@liewenjingque,[0,40],[0,0]);
%he function liewenjingque yiqi shi yong 
[T,y]=ode45(@liewenjingque,[0:0.01:20],[0,0]);
plot(T,y(:,1));
save('egode45','y');