自编近似求定积分的代码

纯萌新向

这几天粗略学了学matlab 突然想编个求定积分的程序

话不多说直接上代码

clear;clc;
syms x;
p(x)=input('请输入要求积分的函数为f(x)=');
a=input('请输入积分下限为='); 
b=input('请写出积分上限为=');
t=linspace(a,b,10000);
f=@(x)p(x);
q=sum(f(t))*(b-a)/10000;
sprintf('该函数在以%f为下限和%f为上限的区域内积分结果为%f',a,b,q)

理论上来说把步距调短 求得的值更精确 但还得看计算机性能

我试了两个简单的函数 基本都花了几秒才出结果 放图