-
-
-
andysmile说:
說一個戀愛故事
這是一個充滿濃郁浪漫的戀愛故事。
哈利和芭芭拉是两个年轻而活力四射的大学生，哈利一直很喜欢芭芭拉，他每天都会在校园里跟着她，他也每天都会在社交媒体上打给她，但芭芭拉却不怎么理他。
一天，哈利决定要向芭芭拉表白，他去了她家，准备给她送上一束鲜花。他走到芭芭拉家门口时，发现芭芭拉也在等他，而且还带着一束鲜花！
哈利被惊讶到了，原来芭芭拉也喜欢他！他笑着说：“看来我们有缘分啊！”之后，两人就一起走在校园里，甜蜜地相依相伴。
从此，这对恋人就一直在一起了，他们经常一起出去旅行，每天也在社交媒体上分享彼此的点滴生活。一年后，哈利向芭芭拉求婚，芭芭拉爱不释手，欣然接受了哈利的求婚。
于是，这对恋人就完成了一段美丽的爱情，他们一起渡过了许多甜蜜的时光，一起走过了美丽的人生旅程。
-
-
peng111说:
求高阶微分方程y”-u(1-y2)y’ +y=0,y(0)= 1,y’(0)=0的数值解，并画出解的相应图形. (μ自己设)
解：
设μ=1
此高阶微分方程化为一阶微分方程组：
$$
\begin{cases}
\frac{dx}{dt}=y\
\frac{dy}{dt}=-1(1-y^2)y+1
\end{cases}
$$初值条件：$x(0)=0,y(0)=1$
用Euler法求解：
令步长h=0.02，则
$(x_0,y_0)=(0,1)$
$(x_1,y_1)=(h,1+hy_1)$
$(x_2,y_2)=(2h,y_1+h(-1(1-y_1^2)y_1+1))$
$(x_3,y_3)=(3h,y_2+h(-1(1-y_2^2)y_2+1))$
$\cdots$
$(x_n,y_n)=(nh,y_{n-1}+h(-1(1-y_{n-1}^2)y_{n-1}+1))$
用MATLAB编程求解得：
h=0.02; x(1)=0;y(1)=1; for n=1:50 x(n+1)=x(n)+h; y(n+1)=y(n)+h*(-1*(1-y(n)^2)*y(n)+1); end x,y结果如下：
把上述结果画在图像上：
