matlab解符号方程输出数值解 限定变量范围

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

x=fsolve(fun,x0)求解fun(x)=0的解,x0是初值,fun是函数,x就是解比如解方你的方程可定义一个函数functiony=f(x)y=(ln(0.2-x)-ln(0.04-x))/

function[x,y,z]=sp(A,B,C)symsxyz[x,y,z]=solve(y*y-z*z-x*x,y+z-A,x*x-B*x-C);

=eval(x)

h=solve('1+a*x+b*x^2+c*x^3=0');a=1;b=2;c=3;eval(h)ans=-0.78390.0586+0.6495i0.0586-0.6495i或者>>subs(h)

symsa1a2c1c2%以下的值你自己改一下aa1=1;cc1=2;aa2=3;cc2=4;[x5,y5]=solve('y5=a1*x5+c1','y5=a2*x5+c2','x5','y5');

程序：symsfy=sqrt(2*(e/4+p/(2*x)-(m^2-1)/(8*x^2)-W/8*x));int(f,x,下限,上限)；OK了.注：一曝晒龙=W

symx;v=int(e^(x^2/(-2)),0,1)vpa(v)

syms t x0=1753000;y0=0;z0=0;vx0=0;vy0=1700;vz0=0;x1=1738000*cosd(10);y1=1738000*sind(10);z

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))

使用subs函数：subs(x1)subs(As)因为没有楼主全部的变量值,这里把相关变量全部赋值为1作为示例：[M,fcd,fy,h0,r0,v1,v2,v3,v4,v5]=deal(1);subs

不是所有符号方程能用solve()解,应减少符号量,最好将已知值代人,只含变量g,这样完全可以得到解.

编程观察了下F1(t)/F0(t)-t函数,图像如下：也就是说a/K需要是个正数才行.你给的数据a是正的,K是负的,无解.给定a,K,计算yita的程序是：%%示例a=1;K=2;Fn=@(n,yit

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

这个方程无解.如果|x|<1,积分项会出现复数,而另一项为实数,无解.而当|x|>=1时,积分项的值大于0.5,另一项最大为0.5,仍然无解. 可以很容易画出方程左端的函数图像：

这只有一个自变量啊,不是偏微分方程.用高数解决下就行了,y=-5x+1.不需要数值解法吧.即使要用数值解也很简单,步骤如下：变化为f(x)+xf'(x)+10x=0,f(x)用y来表示,用不大精确的差

设u=积分(0~x)ydx,那么原来方程就是u''=-0.04u'^2-sin(u)+0,44u''(0)=0,u'=3,u(0)=0(积分上下界都为0)再问：����������룬��ʤ�м����

subs()函数,我也是今天才看到的subs(a+b,a,2)就是把a+b里面的a用2替换subs(x/y,y,3)就是把x/y里面的y用3替换

symsmuLa=1.2825/(0.288*muL);b=muL-0.577215/a;L=muL;cdfL=vpa(exp(-exp(-a*(L-b))))pdfL=vpa(a*exp(-a*(L

symsxya1a2eq1=0.33*cos(a1)+0.3*cos(a1+a2)-x;eq2=0.33*sin(a1)+0.3*sin(a1+a2)-y;x=0.3;y=0.5;s=solve(ev