matlab题用解析法和数值法分别计算sin(lnx)从1到e的积分

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

半椭圆y=3√（1-x^2/4）椭圆周长=2×∫√（1+y′²）dxmatlab>>symsx>>y=3*(1-x^2/4)^(1/2);>>y1=diff(y,x);>>f=(1+y1^2

Pi=[...],Pj=[...]>>Gi=[...],Gj=[...]>>Dij=[...]>>Rij=(sqrt(Pi.*Gi).*sqrt(Pj.*Gj)./Dij.^2)再问：��ʦ�����

这个如果点列很密集的话,数值微分就是差分那你就用循环结构咯forii=1:mforjj=1:mdy(ii)=(y(ii+1)-y(ii))/(x(ii+1)-x(ii));endenddy这样每相邻两

%给你参考一下我当年做的一个类似的程序%cleardt=0.004;t=0:dt:5;a=-pi:0.2:pi;l=100;r=8;d=10;%导弹尺度参数：长度、半径、火焰尺度L=l;D=l/3;c

使用quad有问题>>quad('0.0732*exp(-(x+0.02287).^2/64.90)',0,inf)Warning:InfiniteorNot-a-Numberfunctionvalu

你没有给出参数的具体值,我只好编了.functionhhfsolve(@myfun,[100,10])functionf=myfun(x)T=x(1);t=x(2);c1=1;c2=9;a=8;b=3

你可能误解了解析解与数值解的来源.solve给你的是解析解,你得到解析解以后,给参数赋值,然后再使用eval()函数,即可得到结果（但不是你说的数值解）.如果需要数值解,用积分来说吧,我们有辛普森(s

给出实例看看吧再问：给图了再答：看到图了，这个简单，先求出极值点，然后拟合就行了。你可以把数据给我，我给你做！再问：这个我查了一下，好像也能做了，还有一个，包络是过零点的，RT，帮忙看下怎么拟合，可以

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

a=sprintf('%.2f',10736865691775629/549755813888)b=sprintf('%.2f',-94400051914681359/549755813)y=-600

m=2;cnt=0;k=1;whileabs(m-k)>10e-10cnt=cnt+1;k=m;m=sqrt(m);s(cnt)=m;endformatlongcntS=s'运行结果：cnt=30S=

记得好像是solve('x^2+4*x-100=0')也可以symsx;y=x^2+4*x;solve(y-100,'x')

[v,x,y]=solve('x+a=5','y+b=6','v=x+y')%意思是把v,x,y作为未知量,解方程%结果v=-a+11-bx=-a+5z=-b+6通过helpsolve命令,可以了解更

symst>>ansa=int(exp(-t^2/2),0,3)/(2*pi)^0.5+0.5ansa=1125899906842624/5644425081792261*erf(3/2*2^(1/2

xi_A=[0.20.40.60.81.0];fxi_A=[0.980.920.810.640.38];xi_B=[0.2,0.2+0.08*1,0.2+0.08*11,0.2+0.08*10];%所

程序修改如下：clearall;clc;symsxy;H=[(1+x)*(1+y)/4;(1-x)*(1+y)/4;(1-x)*(1-y)/4;(1+x)*(1-y)/4];fori=1:4forj=

数值型数据是我们通常的数据,但在计算机中一般是有限的长度数字,在一些计算中不能保证完全正确.符号型数据就是把数据看成是符号,就如同我们用x,y代替数据一样,它是完全正确的.

这个很简单,一个命令就够了：ode45(inline('y-2*y^3*sin(t)','t','y'),3,1)先用inline创建一个内联函数描述微

matlab不会把q=acos((x.^2-1+x.^2.*i^2)./(2*x.^2.*i)); 代入f的改为下面的即可i=2;f=@(x)((x.^2-1).*x.^2.*i.*sin(