数值积分

eq=@(y)(2.53e-7*y-5.6e-5)./(0.2872+5.85-0.019*y).*y/1.1.*exp(-0.5*y.^2/1.1);quad(eq,222,300)ans=0

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

你这个式子是没有办法求得解析解的!再问：为什么呢？可否说的清楚点啊，谢谢了。再答：对数函数是没有办法进行解析积分的，这个你应该是知道吧再问：所以我想用数值积分进行求解，譬如给定t的范围（0,100），

clear;clc;a=4;b=2;M=3;f=@(z)quadl(@(x)1./((1+(z-x)/a).*(1+x/b)).^M,0,z);quadl(@(z)arrayfun(@(z)log(1

f=@(x)2./(1+x.^2);%少点I=quad(f,-1,1)

Clear[f,t,x,y];f[t_]:=Piecewise[{{t,00](*求含有f[t]的变上限定积分y*)NIntegrate[y,{x,0,3}](*求y的数值积分*)

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

数值方法当然不能做不定积分啊,不定积分不会有数值结果啊,只有通解int()函数是可以算不定积分的,一般还是很好用的如果解不出来的话,很有可能是真解不出来

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

关键是写出正确被积函数,对其进行采样.x0=0;b=1;f=@(x)sqrt(1+(exp((x-x0).^2/(-2*b^2))).^2);%被积函数xs=(x0:0.01:x0+sqrt(log(

clearall;closeall;clcsigma1=1;mu1=0;mu2=0;sigma2=1;p=0.5;y=-4:0.1:4;symsxf=1/(2*pi*sigma1*sigma2*sqr

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

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

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

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

有20个v,那应该有20个结果吧?而且你的那个函数的定积分是没有解析解的.所以只能求数值解.建议你利用trapz函数,具体的看help.其实也可以用quad但是它需要函数的定义,所以有点麻烦,我不推荐

这是在网上找的,看你应该可以用

n=0.01;x=0:n:1;y=4./(1+x.^2);z=trapz(x,y)n取得越小,积分结果越接近精确值pi.再问：老师的意思是把n的取值看做一个函数直接出来不同n的不同结果不用自己一个一个

quadl(@(x)f(x),x1,x2),如y=sin(x)在（1,5）上的积分,quadl(@(x)sin(x),1,5)

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(