jacrot.m

来自「麻省理工ComputationalMethod」· M 代码 · 共 8 行

function J=jacrot(a,b,d)

tau=(b-a)/2/d;
t=sign(tau)/(abs(tau)+sqrt(1+tau^2));
c=1/sqrt(1+t^2);
s=c*t;
J=[c,s; -s,c];