rosa
Inscrit le: 08 Juin 2004 Messages: 54
|
Posté le: 06 Avr 2005 20:41 Sujet du message: mon pgm |
|
|
voilà mon programme et je ne veux pas de solution négative dans (soln)
Digits:=20:
> DE:=diff(f(y),y$4)-lambda^4*f(y)=0;
> N:=20;
> trial:=sum(a[i]*y^(4+i)*(1-y)^(4+i), i=1..N);
> residual:=lhs(subs(f(y)=trial,DE));
> R:=simplify(residual);
> sys:=NULL:
> for i from 1 to N do
> sys:=sys,int(R*y^(4+i)*(1-y)^(4+i),y=0..1)=0:
> od:
> sys;
> with(linalg):
> A:=genmatrix([sys],[seq(a[i],i=1..N)]);
> data:=det(A)=0;
> assume(lambda,positive);
> about(lambda);
> soln:=fsolve(data,lambda);
> |
|