﻿函数功能:求空间直线与曲面的交点(geometric linear surface intersection)

函数格式:Geolsinter(x1,x2,s)

函数说明:x1,x2均是1*3的矩阵变量(其中三个元素依次代表坐标点相对于x,y,z轴的值).s是一个包含曲面方程的符号变量.s里面的表达式里只能以x、y、z名称分别来代表x、y、z轴的变量.如果曲面与直线有多个交点,只返回一个交点.如果曲面与直线没有交点,程序自动搜索一个比较近的点进行返回.

例子:
a =[ 0.00000000000000   0.00000000000000   0.00000000000000 ]
b =[ 1.00000000000000   1.00000000000000   1.00000000000000 ]
s="x^2+y^2+z^2-3"
Geolsinter(a,b,s)//执行命令
ans =
[ -1.0000000133083   -1.0000000133083   -1.0000000133083 ]//和我们的理想结果(-1,-1,-1)已经很接近了.

//完成 By2011/10/02 17:35:00