Maple & MapleSim分享 http://blog.sciencenet.cn/u/maplesim

博文

Maple 2021中的常微分方程Student[ODEs]学习包,自动显示微分方程的分步求解过程细节

已有 6147 次阅读 2021-6-22 10:54 |个人分类:Maple基础学习|系统分类:科研笔记

Maple 2021新增加了Student[ODEs]微分方程学生包,该学生包可以自动显示计算步骤和详细结果,与以前版本中的微积分、线性代数、数值分析、优化学生包一起辅助老师教学和学生自学。


Student[ODEs]

  命令:ODESteps

  功能:显示ode、ivp或odes的分步求解过程


学生包命令ODESteps支持以下类型的ode和odes:

• 一阶常微分方程

• 一阶IVPs常微分方程初值问题

• 二阶常微分方程

• 二阶IVPs常微分方程初值问题

• Cauchy-Euler 方程

• 级数解

• 特殊函数解

• 常微分方程系统

初始值问题常微分方程系统


示例1

非齐次一阶常系数线性微分方程:

Diff(u(x), x) = c*u(x) + x^2


restart;

with(Student[ODEs]);


ode1 := Diff(u(x), x) = c*u(x) + x^2; 

image.png

ODESteps(ode1);   #显示计算细节

image.png

示例2

with(Student:-ODEs);

ode2 := diff(y(x), x) - y(x) - x*exp(x) = 0;

image.png

ODESteps(ode2);

image.png

示例3:一阶IVPs

restart

with(Student:-ODEs);

ivp1 := {t^2*(z(t) + 1) + z(t)^2*(t - 1)*diff(z(t), t) = 0, z(3) = 1};

image.png

ODESteps(ivp1);

image.png

示例4:二阶ODEs

restart;

with(Student:-ODEs);

ode4 := 2*x*diff(y(x), x) - 9*x^2 + (2*diff(y(x), x) + x^2 + 1)*diff(y(x), x, x) = 0

image.png

ODESteps(ode4)

image.png

示例5:二阶IVPs

with(Student:-ODEs);

ivp5 := {diff(y(x), x, x) - diff(y(x), x) - x*exp(x) = 0, eval(diff(y(x), x), x = 0) = 0, y(0) = 1};

image.png


ODESteps(ivp5)  #显示计算细节

image.png


示例6:Cauchy-Euler方程

with(Student:-ODEs);

ode6 := x^2*diff(y(x), x, x) - 4*x*diff(y(x), x) + 2*y(x) = 0;

image.png

ODESteps(ode6)

image.png

示例7:级数解

restart;

ode7 := x^2*diff(y(x), x, x) + x*diff(y(x), x) + 5*x*y(x) = 0

image.png

Student:-ODEs:-ODESteps(ode7)

image.png

示例8:特殊函数解

restart

with(Student:-ODEs);

ode8 := x^2*diff(y(x), x, x) + 4*x*diff(y(x), x) + (25*x^2 - 9)*y(x) = 0;

image.png

ODESteps(ode8)

image.png

示例9:常微分方程系统

restart

with(Student:-ODEs);

high_order_ode1 := diff(y(x), x, x, x) + 3*diff(y(x), x, x) + 4*diff(y(x), x) + 2*y(x) = 0;

image.png

ODESteps(high_order_ode1)

image.png

示例9:Systems of ODEs with IVP

restart;

with(Student:-ODEs);

high_order_ivp1 := {diff(y(x), x, x, x) + 3*diff(y(x), x, x) + 4*diff(y(x), x) + 2*y(x) = 0, eval(diff(y(x), x), x = 0) = -1, eval(diff(y(x), x, x), x = 0) = 2, y(0) = 1};

image.png

ODESteps(high_order_ivp1)

image.png


你也可以下载免费手机app“Maple计算器”扫描获取相同的功能。




https://blog.sciencenet.cn/blog-516836-1292193.html

上一篇:数学软件Maple 2021发布
下一篇:求偏微分方程的数值解(有限差分法、配置法、Galerkin方法)
收藏 IP: 210.13.118.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-23 04:47

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部