NURBS曲线曲率计算

function curvature = clcCurvature（nurbsut）
%clcCurvature - 计算NURBS曲率
%
% Syntax: curvature = clcCurvature（nurbsut）
%       input:
%               nurbs:nurbs曲线的表达式，来源于nrbmak.m
%               ut:需要计算的曲率的点
%
% Description: 
%       该算法来源于 The nurbs book 算法：A4.2

dnurbs = nrbderiv（nurbs）;
ddnurbs = nrbderiv（dnurbs）;

[p w] = nrbeval（nurbsut）;
[p1 w1] = nrbev