matlab求变步长积分

把f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);改为f=a^2*(1-e^2)*(cos(x))./((1-((e^2)*(sin(x)).^2)

积分区间[a,b]等分为n段,积分步长h＝(b-a)/n

s1='(x+1)^10*exp(-x)'s1=(x+1)^10*exp(-x)>>int(s1,0,1)ans=-26813184*exp(-1)+9864101>>numeric(int(s1,0

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

楼上的是对的,可是要注意大小写啊set(gca,'Xtick',[xmin:step:xmax])如：标记0pi2pi并将坐标轴改成0pi2pi,代码如下：set(gca,'Xtick',[0pi2*

积分结果是f(0)再问：是直接应用δ函数的性质吗？再答：嗯

最好给出具体的函数表达式.一般的做法差不多可以写成：fminunc(@(yz)quadl(@(x)f(x,yz(1),yz(2)),a,b),yz0)其中,yz0是y和z的初值,a和b是积分上下限.&

x=0.6时,y=68.581359956511627574834151414947closeallclear,clcx=0:0.05:1;y=zeros(1,length(x));symsyYex1

我有excel的程序,是用复化simpson求积的工具.要启用宏.前几天刚考完计算方法,复习完了就做了个小程序.你要的话回我一个邮件,我发给你.

symsaf=((2.72.^a-a-1).^(-1/2)-1);g=int(f,a,a,3)结果：g=int(1/((68/25)^a-a-1)^(1/2)-1,a=a..3)

这条语句太长,分成几个表达式写,是括号错误.我调了一下,有两个地方：1.mu=0.4l;%滑动摩擦系数,即μ应为0.41；你是0.4L.2.倒数第二句：f=dW*ds,这是symbol对象,改为inl

http://numericalmethods.eng.usf.edu/mws/gen/07int/mws_gen_int_txt_gaussquadrature.pdf

%bydynamic%2009.2.%积分Matlab中提供了数值和解析解法两种1.数值方法>>fun=@(x)738.5*x.^2.*exp(-0.00105*x.^2);>>quadl(fun,0

其实你可以考虑尝试用数组的算法,毕竟MATLAB并不擅长循环运算.如果不是特别复杂的运算,数组的算法可能并不难想出来.如果你非要用变步长的循环,我想到一个折中的办法,因为你提前知道了m应该取那些数值,

走的距离÷走的步数=平均步长

楼主是不是编译的时候显示sin没定义?编译的时候给最后面加个-lm就行了

symsxtf=abs(sin(cos(t)));%需要积分的公式I=int(f,x,0,2*pi)%int(f,x,a,b)表示定积分,f是方程式,x表示参数,a表示下限,b表示上限sim=simp

我的建议是把e,p的值代入int()中,可以计算,这样是有值的

#include"stdio.h"#include"math.h"doublef(doublex){returnsin(x);}doubleJifen(doublea,doubleb)

f=inline('sin(x)')>>g=@(x)x.*f(x)>>quad(g,1,2)