matlab二分法求解一次方程的解法

clc;clear;y=sym('a+b*x+c*x^2+d*x^3+e*x^4+f*x^5')%是不是c*x^2啊y1=diff(y)y2=diff(y,2)y3=diff(y,3)eq1=subs

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再问：粘贴了您的

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

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

x=5y=3z=4

a1x+b1y=c1a2x+b2y=c2当a1b2-a2b1≠0,b1a2-b2a1≠0时x=(c1b2-c2b1)/(a1b2-a2b1)y=(c1a2-c2a1)/(b1a2-b2a1)当a1b2

function[x,nitr]=bisecton(f,a,b,d)%fÊÇ´ýÇ󷽳̵Ä&#

兄弟,西理工的?

a=1;b=2;f=@(x)x^3-x-1;c=(a+b)/2;whileabs(b-a)>1e-5iff(c)*f(b)

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

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

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

data=[];fora=1:19forb=1:a-1forc=1:b-1ford=1:c-1fore=1:d-1ifa+2*b-5*c+8*d-4*e==15data=[data;a,b,c,d,e

设在5元的基础上涨价X元,Y为利润,则由已知可得销售量为：500－20XY＝（500－20X）*（5+X－4）＝（500－20X）*（X+1）＝500X－20X^2－20X+500＝1760＝〉-20

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

把每个未知数分别用引号括起来,不是放在一起括起来.'a11','a12','a13','a14','a21','a22','a23','a24

少了一个单引号.>>[x,y]=solve('0.99820050*(1+80.01*x+80.01^2*y)=(1+75.9*x+75.9^2*y)','0.998005982*(1+99.99*x