DAY 19&20 8.6-8.9

8.6 常微分方程的符号解


8.7 平面曲线族的包络线
若有包络线，则必定满足 f(x,y,a)=0 df/da=0
例：8-38

syms a x v g
y=xtan(a)-gx^2/2/v2/cos(a)^2;
diff(y,‘a’)
ans1=ans;
solve(ans1,‘x’)
x=ans(2)
x=simplify(x)
y=v^{2/g*(1-1/2/sin(a)}2 );
g=9.8;v=100;
parabola_cover
未定义函数或变量 ‘parabola_cover’。

为双精度值的数值、日期时间、持续时间或数组。
注意 简化函数无simple 而是simplify

parabola_cover是自己定义的函数

g=9.8;v=100;
a=pi/50:pi/50:24pi/50;
x =(v^2cos(a))/(g*sin(a));
y=v^{2/g*(1-1/2./sin(a).}2 );
plot(x,y)
不知道为什么这个程序画不出包络线

8.8常微分方程的数值解

8.9 差分方程求解
例8-4