coeflagu.m

来自「various matlab programs to slove various」· M 代码 · 共 5 行

function [a, b] = coeflagu(n)
if (n <= 1), disp(' n must be > 1 '); return; end
a=zeros(n,1); b=zeros(n,1); a(1)=1; b(1)=1;
for k=2:n, a(k)=2*(k-1)+1; b(k)=(k-1)^2; end