SG Fórum - Maple

Fórum / Felhasználói programok

xDJCx 2011. okt. 12. 19:29 | válasz | #527

Van benne több hiba. Javítottam egy-kettõt,amit megtaláltam. Már szintaktikailag helyes, de tartalmilag még nem.

Jeleztem a javított sorokat, és azt is ami kérdéses.

Pl. boundRoots eljárás még biztos hibás: pl. a bb4: kezdetû sorban a t változó nincs elõtte definiálva és itt mindenképpen kellene neki érték.

polyops:=module()

export leen, height, qnorm, boundRoots , boundRootsDense;

options package;

leen := proc (p::(polynom(constant)), x::name)

convert(map(abs, [coeffs(expand(p), x)]), `+`)

end proc;

height := proc (p::(polynom(constant)), x::name)

max(coeffs(collect(p, x, distributed, abs), x))

end proc;

qnorm := proc (p::(polynom(constant)), x::name)

sqrt(convert(map(z->z*conjugate(z) , [coeffs(expand(p), x)]), `+`))

end proc;

boundRoots := proc (p::(polynom(constant)), x::name)

local q, d, s, n, t, b, b1, b2, b3, bb4;

q := collect(p, x, distributed, evalf);

d := degree(q, x);

s := lcoeff(q, x);

n := nops(q);

if d = ldegree(q, x) then RETURN(0) end if;

if s <> 1. then q := expand(q/s) end if;

b1 := 1+polyops:-height(q-x^d, x);

b2 := max(1., leen(q-x^d, x));#javítva leen ?

b := min(b1, b2);

b3 := 0;

bb4 := (n*abs(lcoeff(t, x)))^(1./(d-degree(t, x)));# mi az a t?

for t in q-x^d while b3 < b do

b3 := max(b3, bb4);

end do;

userinfo(3, polyops, [b1, b2, b3]);

min(b, b3);

end proc;

boundRootsDense:=proc(p::polynom(constant), x::name)

local q, d, s, k, b, b4, b5;

q:=collect(p,x, distributed,evalf);

if ldegree(q)>0 then

q:=expand(q/(x^(ldegree(q))));

end if;

if type(q, constant) then

RETURN(0);

end if;

d:=degree(q,x);

s:=lcoeff(q,x);

if s<> 1. then

q:=expand(q/(s))

end if;

b:=polyops:-boundRoots(q,x);

b4:=convert([abs(coeff(q,x,0)),seq(abs(coeff(q,x,k-1)-coeff(q,x,k)),k=1..d)],`+`);

if nops(q)<d+1 then b5:=infinity;

elif type(q, polynom(positive,x)) then

b5:=max(seq(coeff(q, x, k-1)/coeff(q, x, k), k=2..d));# javítva;

else

b5:=max(abs(coeff(q, x, 0)/coeff(q, x, 1)), seq(2*abs(coeff(q, x, k-1)/coeff(q, x,k)), k=2..d));#javitva

end if;

#end if; javitva

userinfo(3, polyops, [b, b4,b5]);

min(b, b4, b5);# javítva

end proc;

end module;

Válasz 'Gimre' üzenetére (#526)