egode45.m

内燃机转子仿真

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');