function sssss
clear all
clc
[t,y]=ode45(@odefun,[0,1],(1,1));
y1=y(:,1);
y3=y(:,2);
dist(t,y)
function [dydt]=odefun(t,y)
global dy1dt dy3dt
y1=y(1);
y3=y(2);
for i=1:length(t);
y2(i)=1+(t-y1(i));
end
dy1dt=t+y2;
dy3dt=t-y2;
dydt=[dy1dt dy3dt]';
为什么老是报错啊 |