怎么用ode23求解一阶方程?
看了一本很基础的书,就是不知道怎么输入表达式才行
书是这么说的:
“dy/dt=cos(t)
我们很容易求得方程的积分为y(t)=sin(t)+C,它将可以用来检验我们得到的数值解.
首先我们定义函数,我们创建一个.m文件,输入下面的内容.
functionydot=eq1(t,y)
ydot=cos(t);
我们的函数称为eq1.让我们求0≤t≤2π内的解并假设y(0)=2.调用的语句是:
>>[t,y]=ode23('eq1',[02*pi],2);”
我按说着的步骤File-New-M-File,然后输入
functionydot=eq1(t,y)
ydot=cos(t);
接着保存到安装目录下的work文件夹下,左侧的currentdirectory也切换到work文件夹了,然后就是在commandwindow下输入:
[t,y]=ode23('eq1',[02*pi],2);
提示没有定义eq1,直接输入eq1文件名,又提示Inputargument"t"isundefined.
请问如何输入命令?