matlab 中用二分法求bessel 方程的根

先建立二分法的fun.m文件,代码如下：functionfun(a,b,e)%f是自定义的函数%a为隔根区间左端点,b为隔根区间右端点,e为绝对误差限ifnargin==2e=1.0e-6;elsei

i=12;s=1;while(i>0)s=s*i;i=i-1;end

f=inline('x^2-x-2');>>[c,err,yc]=bisect(f,0,3,0.01)c=2.0010err=0.0059yc=0.0029-----------%使用二分法求解上面超

functiontesta=1;b=2;tol=1e-4;N=10000;k=0;fa=f(a);fork=1:Np=(a+b)/2;fp=f(p);if(fp==0||(b-a)/2再问：粘贴了您的

输入a的数据；x=a(:,1);y=a(:,2);plot(x,y);

代码如下:functionrtn=newton1(fx,dfx,x0,tol,N)%NewtonMethod%Thefirstparameterfxisaexternalfunctionwithres

这是源代码:在matlab中保存为:bisection.mfunction rtn=bisection(fx,xa,xb,n,delta)% 二分法解方程% fx是由方程

可以,只要在二分法外面,加一个循环就可以了.比如,以i为变量从1到1000,只要f(x(i-1))*f(xi)

这个有例子区间二分法：与对分查找法相同1区间二分法求出的仅仅是方程的一个单根,如果方程有重根或者多个根时,在做区间二分法时就会出现分叉,这样方程有几个根,就会产生几个实数序列,每一个实数序列的极限便是

用二分法,你得先找出一个根所在的区间function[a,b]=findbracket(f,x0)%fisthefunctionevaluated%x0isthestartingpoint%aisth

程序如下：clear,clc;a=0;%a=input('inputa:');b=1;%b=input('inputb:');err=10^-5;y1=a*exp(a)

对于求平方根,变成方程模式为f(x)=x^2-a,即求此方程的实根；下面编写了两个function函数,可以直接调用.二分法：functionx=sqrt_bisect(a)f=@(x)x^2-a;i

那调用两次不就ok了.区间你是用a、b定义的嘛.

建议你在第三行插入一个end,以完成函数的定义而且函数的定义只有在m文件里面才能做到.（这很重要!）还有,好像没有与while相应的end对了,你的函数在[-pi/2,pi/2]里面只有一个解0,所以

function [x,fval]=myerfen()clc;f=@(x) sqrt(x^2+9)-tan(x);a=0;b=pi/2-eps;tol=1e-4;while&nbs

在matlab中使用fsolve函数就行了再问：不知道具体的该怎么使用呢？是直接fsolve（inline('H'),x0）就可以么？我初学，很多不懂，多谢再答：你可以例如x=fsolve(@(x)s

你是在说符号运算么,亲?再问：恩，其实我感觉应该不太算是符号运算。。。因为整个运算，在源程序中是符号运算，但是因为我是计算的物体有具体的尺寸限制，而且我是需要得到具体的数值的，所以我会另外进行一系列的

functiontestclearclc%实验方程：3*x.^2+x+2*exp(x)=0%原函数f=@(x)3*x.^2+x-2*exp(x);%导函数df=@(x)6*x+1-2*exp(x);%

A=[123-4-2;-3-4-12135;2100-310;4149-137];%A为系数矩阵fori=2:4forj=i:5fork=1:(i-1)A(i,j)=A(i,j)-A(i,k)*A(k

举个简单例子：问题1：求f(x)=1-x-sinx=0在【0,1】的根误差不超过0.5*10^(-4)clc;cleara=0;b=1;fa=1-a-sin(a);fb=1-b-sin(b);c=(a