注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

徐小湛的博客

高等数学课件、函数图形动画、数学软件应用

 
 
 

日志

 
 

微分方程作图   

2012-05-17 11:39:52|  分类: 高等数学图形动画 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
微分方程作图 - calculus - 徐小湛的博客

 

微分方程作图 - calculus - 徐小湛的博客
 
 
with(DEtools): 调用微分方程工具
with(plots):  调用绘图工具
de:=diff(y(x),x)=2*x*y(x);  定义微分方程
fxc:=DEplot(wffc,y(x),x=-2..2,y=-2..2): 画斜率场
jfq:=contourplot(y/exp(x^2),x=-2..2,y=-2..2): 画积分曲线
display(fxc,jfqx);
 微分方程作图 - calculus - 徐小湛的博客

 

微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=diff(y(x),x)=2*x*y(x):dsolve(wffc);
fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(y/exp(x^2),x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);

 

微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=diff(y(x),x)=(cos(y(x))-y(x)*cos(x))/(x*sin(y(x))+sin(x)-1):
dsolve(wffc);fangxiangcang:=DEplot(wffc,y(x),x=0..4*Pi,y=0..4*Pi,thickness=2):
jifenquxian:=contourplot(y*sin(x)-x*cos(y)-y,x=0..4*Pi,y=0..4*Pi,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);

 

微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=diff(y(x),x)=(cos(y(x))-y(x)*cos(x))/(x*sin(y(x))+sin(x)-1):
dsolve(wffc);fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(y*sin(x)-x*cos(y)-y,x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);

 

微分方程作图 - calculus - 徐小湛的博客

 

微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=x*diff(y(x),x)+y(x)=sin(x):dsolve(wffc);
fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(x*y+cos(x),x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=(x-y(x)^3)*diff(y(x),x)+y(x)=0:dsolve(wffc);
fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(x*y-y^4/4,x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);
 
微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):
DEplot((x-y(x)^3)*diff(y(x),x)+y(x)=0,y(x),x=-2..2,
y=-2..2,[[y(0)=1],[y(0)=0.3],[y(0)=1.5],[y(0)=-0.5],[y(0)=-1],[y(0)=-1.5]],linecolor=[blue,black,gold,navy,green,maroon],
color=violet,stepsize=0.01,scaling=constrained);
 
微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):with(plots):
wffc:=diff(y(x),x)=2*x*y(x):dsolve(wffc);
fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(y/exp(x^2),x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);
 
微分方程作图 - calculus - 徐小湛的博客
 
with(DEtools):
DEplot(diff(y(x),x)=2*x*y(x),y(x),x=-2..2,
[[y(0)=-1],[y(0)=-0.5],[y(0)=0],[y(0)=0.5],[y(0)=1],[y(0)=1.5]],
y=-4..4,linecolor=[gold,black,blue,red,brown,green],
color=grey,stepsize=0.01,scaling=constrained);
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=3*x*y(x)^2*diff(y(x),x)=x^3+y(x)^3:dsolve(wffc,implicit);
fangxiangcang:=DEplot(wffc,y(x),x=-2..2,y=-2..2,thickness=2):
jifenquxian:=contourplot(x^2-2*y^3/x,x=-2..2,y=-2..2,contours=20,color=blue,thickness=2):
display(fangxiangcang,jifenquxian);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)-diff(y(x),x)-2*y(x)=0:
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-2..2),_C2=-2..2)]:
plot(toplot,x=-1..1,y=-10..10,thickness=3,color=red);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)+2*diff(y(x),x)+y(x)=0:
tongjie:=dsolve(wffc,y(x)):tejie:=dsolve({wffc,y(0)=4,D(y)(0)=-2},y(x)):
tongjiequxian:=[seq(seq(rhs(tongjie),_C1=3..5),_C2=0..4)]:
p1:=plot(tongjiequxian,x=-2..2,y=0..8,thickness=1,color=red):
p2:=plot(rhs(tejie),x=-2..2,y=0..8,thickness=5,color=blue):
display(p1,p2,scaling=constrained);
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)+2*diff(y(x),x)+y(x)=0:
tongjie:=dsolve(wffc,y(x)):tejie:=dsolve({wffc,y(0)=4,D(y)(0)=-2},y(x)):
tongjiequxian:=[seq(seq(rhs(tongjie),_C1=2..6),_C2=0..8)]:
p1:=plot(tongjiequxian,x=-2..2,y=0..8,color=blue):
p2:=plot(rhs(tejie),x=-2..2,y=0..8,thickness=5,color=red):
display(p1,p2);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)-4*diff(y(x),x)+13*y(x)=0:
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-2..2),_C2=-2..2)]:
plot(toplot,x=-2..1.5,y=-20..20,thickness=3,color=red);
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)-4*diff(y(x),x)+13*y(x)=0:
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-3..3),_C2=-3..3)]:
plot(toplot,x=-2..1,y=-20..20,thickness=2,color=blue);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)+y(x)=2*x^2-3:
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-3..3),_C2=-3..3)]:
plot(toplot,x=-2..2,y=-12..6,thickness=2,color=blue);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)-2*diff(y(x),x)-3*y(x)=exp(-x):
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-3..3),_C2=-3..3)]:
plot(toplot,x=-2..2,y=-12..12,color=blue,thickness=2);
 
微分方程作图 - calculus - 徐小湛的博客
 
微分方程作图 - calculus - 徐小湛的博客
 
wffc:=diff(y(x),x$2)-2*diff(y(x),x)+y(x)=(1+x)*exp(x):
tongjie:=dsolve(wffc,y(x)):
toplot:=[seq(seq(rhs(tongjie),_C1=-3..3),_C2=-3..3)]:
plot(toplot,x=-4..3,y=-10..10,color=blue);
 
 
  评论这张
 
阅读(1327)| 评论(3)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016