bisect

这个是函数，不是直接运行的，要给一些参数，然后在命令行或者其他代码里调用它。比如你可以在命令行运行类似这样的agui_bisect(@sin, 0, 1, 1e-6) 其中agui_bisect是你这个函数保存时的文件名，有时保存的时候跟里面写的不一致，就要以文件名为准。当然，你这个保存的是agui_bisect.m，跟里面function后面定义的是相同的，就没问题了。 @sin是函数的句柄，这里用的是matlab的内部函数，正弦函数。可以用你需要的函数替换。做参数传递时，m文件那种函数需要加@，另外还有种inline函数，那个不加@。 后面还有3个参数，看样子像区间的边界a,b，以及精度e。 这些参数都可以换成别的值试试。 总之这种function开头的m文件，一般都不是直接运行的，大多是需要参数的，要由别的程序来调用。

中垂线

right bisector正确的平分线.很高兴为你解答！如有不懂，请追问。 谢谢！

你的code是找n以下所有的质数Q2: 你的循环条件就是sieve不能是空，所以不能是>1，一旦sieve空了，也就是sieve里的元素都能被p整除，循环停止。Q3:这是R语言，没有0的indexQ4:return是函数返回值，primes就是储存n以下所有质数的向量，不返回这个返回啥。Q5:stop就是当不满足if里的条件n>2时，错误信息提示：你的input value 必须大于等于2Q6: 感觉这是道English 的题目，你是不是直接翻译过来的。直接把原题贴上来。Q7: break加在sieve <- sieve [ (sieve %% p )!=0]的后面，整个code 的思想是2到100（我假设n=100），把2拿出来放到primes里，然后去掉2的倍数；接着把3push到primes里，之后去掉3的倍数；为什么没有4直接到5？因为在第一步4就被抹掉了(2的倍数)，接着5拿出来，然后去掉5的倍数。。。一直到最后一个数97 其实可以写个break条件，当p>=sqrt(100)就可以停止了，因为当时sieve里的2到100已经只剩下质数，为何你得自己想想。所以在第8行后面加：if(p>=sqrt(n)) break #大大节省运算时间 原来是从2试到100，现在就2到10，快了好几倍当然你输出的结果也得改一改 return(c(primes,sieve))两分法求根你google一下不就有了，现在都2013年，还有维基百科