求微分方程的符号解,并与数值解进行比较

dsolve('x^2*D2y+4*x*Dy+2*y=0','y(1)=2','Dy(1)=-3','x') ans 

结果：代码：clearallclcf=@(x,y)([y(2);   0.357*y(1)-0.1905*y(1)*y(2)]);[x,Y]=ode45(f,[0100]

请参考ode45函数的用法将所求值设为终值即可

%欧拉法解一阶常微分方程% y'=xy^(1/3)f = inline('x*y^(1/3)','x','y');

初值不够呀,应该有四个,题目能不能贴出来?

……一定要用matlab吗,只会mathematica的路过.而且啊,你一面说要求方程的数值解,一面又只给一个全是字母系数的方程,你是要怎样……再问：呃，abcdef都是常量，可以随便假定的，就赋1、

function[x,y]=Classical_RK4(odefun,xspan,y0,h,varargin)x=xspan(1):h:xspan(2);y(1)=y0;fork=1:length(x

一阶线性非齐次微分方程y'+p(x)y=q(x),通解为y=e^[-∫p(x)dx]{∫q(x)e^[∫p(x)dx]dx+C}用的方法是先解齐次方程,再用参数变易法求解非齐次.《高等数学》教科书上都

你这个属于边值问题.得用打靶算法.如果想简单,你令一个变量为x=y-2001,就变成正常的微分方程.

y=dsolve('D3y=D2y-Dy-y+t^2','y(0)=0','Dy(0)=1','D2y(0)=-1');tt=linspace(0,10,1001);%表示t从0到10fori=1:1

这个方程可以求符号解啊,虽然答案里用到了特殊函数,但是稍微带入下具体数值求个数值解也是可以的：eqn=r'[t]^2==cr[t]^2+2k^2r[t]-h^2/r[t];DSolve[eqn,r[t

y=dsolve('x*(D2y)+(1-2)*(Dy)+y=0','y(0)=0','Dy(0)=0','x')y=C6*x*besselj(2,2*x^(1/2))

首先建立vdp1.m文件functiondx=vdp1(t,x)dx=zeros(2,1);dx(1)=x(2);dx(2)=(10000*x+(96.321*exp(-t/378.47108)+28

trapz函数可以对只知道离散数值的函数做积分,示意：trapz(ts,ys);

f='sin(s+2*x)';>>int(f)%符号解ans=-cos(s+2*x)/2>>int(f,pi/2,pi)%数值解ans=-cos(s)

1、把G=1/3*((5*Pp-2*P1)/(P1-2*Pp)-P1*c1^2/Pp*cp^2);改成G=1/3*((5*Pp-2*P1)/(P1+2*Pp)-P1*c1^2/(Pp*cp^2));表

改动：functionS=history(t)S=ones(3,1);改为S=ones(4,1);OVER!PS.多翻翻help!

y'-xy/(1+x²)=0的解能用分离变量法求出来,是lny=1/2ln(1+x²)+C就是y=k√(1+x²)再设y'-xy/(1+x²)=x+1的通解是y

NDSolve里的语法错了,是{t,0,10}而不是{t,10}；另外对于常微分方程,所需的边界条件数目,为它的最高阶导数的阶数,并且边界的导数阶数要小于方程的阶数.也就是说你这里的条件多了一个,把那

x*dy/dx=ylnydy/(ylny)=dx/xlnlny=lnx+Alny=x*e^A=B*xy=e^(B*x)=(e^B)^x=C^x由x=1时y=2,C=2故特解是y=2^x