www.gusucode.com > mbcmodels 工具箱 matlab 源码程序 > mbcmodels/@localsurface/evalpev.m

    function p= evalpev(x,L,varargin);
%EVALPEV

%  Copyright 2000-2004 The MathWorks, Inc. and Ford Global Technologies, Inc.



FX= CalcJacob(L,x);

Ri= var(L);

% calculate inverse and use multiplication
pev= FX*Ri;

p= sum(pev.^2,2);