[r,c]size(data)
newdatazeros(r1,c);
newdata(1:151,:)data(1:151,:);
newdata(153:end,:)data(152:end,:);
%newdata
设xsqrt(10)为解析解;而x3.1623为数值解
step：动态系统的阶跃响应绘图。
具体详细具体解释：
：计算出一个相册系统的阶跃服务控制器。在状态空间的情况下,简单假设数码宝贝传说状态为零。当它没有输出参数时,内部函数这个函数在屏幕上的画出阶跃发令。
例如a[123;000;234;000;789];从中我们可以猜想都为0的是第2行和第4行。首先我们也可以按行从小排列，全是0的必然是第一行[c,ia,ic]unique(a,rows);所以从ic数列中去寻找为数字为1的即为所求oppofind(ic1)一加7(ic1)ans24
二分法基本思路
像是地，这对函数f(x),要是修真者的存在实数c,当xc时，若f(c)0,这样的话把xc叫做什么函数f(x)的零点。解方程即那些要求f(x)的所有零点。
假定f(x)在区间（x，y）上后
先找到a、b不属于区间（x，y），使f(a)，f(b)异号，只能证明在区间(a,b)内当然有零点，然后求f[(ab)/2],
现在假设不成立f(a)lt0,f(b)r260,altb
1如果不是f[(ab)/2]0，该点就是零点，
假如f[(ab)/2]lt0,则在区间（(ab)/2，b)内有零点，(ab)/2gta，从①就开始再继续使用
2中点函数值判断。
假如f[(ab)/2]a80，则在区间(a,(ab)/2)内有零点，(ab)/2ltb，从①又开始再使用中点函数值判断。
那样的话就也可以不断地将近零点。
你每次把f(x)的零点的地方小区间收缩起来一半的方法，使区间的两个端点逐步降低缩近函数的零点，以解值零点的近似值，这种方法叫做什么二分法。
从以上更说明，隔一段时间运算后，区间长度下降一半，是线状收敛。别外，二分法不能计算复根和重根。
二分法步骤
用二分法求方程的根的近似值的步骤
1若对此有，则在内至少有一个根。
2取的中点可以计算
3若则是的根，停止下来计算，
启动后输出来结果
若则在内至多有一个根。取；
若，则取；
④若（为预做决策变量的要求精度）后退换算，运行后输出结果，则难，直接返回步骤1，乱词步骤1,2,3
二分法mtalab程序
symsx
funinput((然后输入函数形式)fx)
ainput(（键入二分法下限）a)
abstractvector(（然后输入二分法上限）b)
zipfile(输入误差限d)%二分法求根
%finline(x^2-4*x4)
%改是需要求解答的inline函数的函数体
finline(run)%修改不需要求解释的inline函数的函数体
eb-ak0
whileegtd
c(ab)/2
iff(a)*f(c)lt0
bc
elseiff(a)*f(c)r260
ac
catch
acbc
end
ee/2kk1
end
x(ab)/2
x%x为答案
k%k为次数
