Voir le sujet précédent :: Voir le sujet suivant |
Auteur |
Message |
rosa
Inscrit le: 08 Juin 2004 Messages: 54
|
Posté le: 26 Juin 2004 3:10 Sujet du message: tjs des pbm en programmation |
|
|
salut
tjs avec mon pbm de syst, apres avoir resolu mon syst j ai bien sure plusieurs solutions sous forme de fct de x^4,x^6....ce qe je cherch maintenant c'est cmt declarer qe je dois deriver ts ces resultats par rapport à x^2 (derivee premiere: comme si j fais un chgt x^2=y) et resoudre apres ces derivees=0 pour trouver alors directement x et non pas x^2
j'espere q j ai pu bien expliqer ma qestion car c un peu compliq mais qoi faire!!! j essay tjs merciiiiiiiiiiii d'avance pr votre aide et j suis tjs à votre disposition pr mieux eclaircir ... |
|
Revenir en haut de page |
|
|
nono56 Invité
|
Posté le: 26 Juin 2004 6:52 Sujet du message: Re: |
|
|
Salut Rosa,
Code: |
> f:=x-> (65+73*x^2+105*x^4+215*x^6)/(56*x^2);
2 4 6
65 + 73 x + 105 x + 215 x
f := x -> 1/56 ----------------------------
2
x
> g:=y->f(sqrt(y));
g := y -> f(sqrt(y))
> g(y);
2 3
65 + 73 y + 105 y + 215 y
---------------------------
56 y
> dg:=y->normal(diff(g(y),y)):
> # récupération des valeurs réelles positives annulant la dérivée de f
> map(y->sqrt(y),fsolve(dg(y)));
> #on suppose x>0, ce qui n'est pas genant, f étant paire
dg := y -> normal(diff(g(y), y))
0.6800720243
> # vérification graphique
> plot(f(x),x=0..10,view=[0..10,0..20],numpoints=1000);
|
Je pense que ça répond à ton pb.
Bye. |
|
Revenir en haut de page |
|
|
rosa
Inscrit le: 08 Juin 2004 Messages: 54
|
Posté le: 26 Juin 2004 14:57 Sujet du message: tjs pbm d syst |
|
|
oui merci ça peut repondre à mon pbm mais j veux surtout savoir cmt declarer ceci pr ttes les solutions car à chaq fois j aurai plus d solutions dependant du nbre d'inconnus!!!! |
|
Revenir en haut de page |
|
|
nono56 Invité
|
Posté le: 26 Juin 2004 16:23 Sujet du message: Re: |
|
|
Tu fais pareil pour chaque solution
A+ |
|
Revenir en haut de page |
|
|
rosa
Inscrit le: 08 Juin 2004 Messages: 54
|
Posté le: 27 Juin 2004 0:48 Sujet du message: tjs le meme pbm |
|
|
le pbm c qe je ne sais combien de solution reel j' en aurai
je veux savoir cmt signaler dés le debut pr q ça s'appliq pr ttes solution reelle sans connaitre combien y en a
par exple pr un syst d'ordre 2 j ai trouvee 2 sol distinctes :
sol:= [{s = 3/28*(47040+6608*k^2+252*k^4+9*k^6)/k^2}, {s= 11/12*(181440+7632*k^2+108*k^4+k^6)/k^2}] mais pr + d'inconnues bien sure j aurai + d solutions j n peux pa y appliq pour chaq solut et meme si le cas cmt j peux indiqer ceci au debut d mon pgm pr avoir à la fin les solt directes j espere q vs allez m comprendre
merci pr l tps q vs m consacrez |
|
Revenir en haut de page |
|
|
nono56 Invité
|
Posté le: 27 Juin 2004 7:47 Sujet du message: Re |
|
|
Tu mets au départ toutes les solutions dans un ensemble ou une liste s,
exemple s:={solve(sys)} ou s:=[solve(sys)], ensuite nops(s) donne le nb total de solutions.
A+ |
|
Revenir en haut de page |
|
|
rosa
Inscrit le: 08 Juin 2004 Messages: 54
|
Posté le: 28 Juin 2004 0:56 Sujet du message: desolee tjs mme pbm |
|
|
voici la partie d mon pbm pr etre bien clair:
sys:=NULL:
> for j from 1 to N do sys:=sys,int(R1[1]*z*(1-z)*z^(j-1),z=0..1)end do:
> for i from 1 to N do
> sys:=sys,int(R1[2]*z*(1-z)^2*z^i,z=0..1)end do:
> {sys};
> with(linalg):
> A:=genmatrix([sys],[seq(a[i],i=1..N),seq(b[j],j=1..N)]);
> data:=det(A)=0:soln:=solve({data},{Ra});
> poly1:=[eval(soln,sigma=0)];
"sigma est un parametre q'on prend nul pr simplifier les calculs"
apres pr N=2 par exple, j aurai comme resultat :
poly1 := [{Ra = 3/28*(47040+6608*k^2+252*k^4+9*k^6)/k^2}, {Ra = 11/12*(181440+7632*k^2+108*k^4+k^6)/k^2}]
moi j cherch à n pas prendre à chaq fois les solut et poser f:=k-> 3/28*(47040+6608*k^2+252*k^4+9*k^6)/k^2
et puis
g:=y->f(sqrt(y)); ....
j veux q tt c calcul soit inclu directement ds le pgm
qel instruction dois-je utiliser?
car pr N=82 par exple j n'en peux plus!? cmt faire donc!
merci infiniment et d'avoir pris d votre tps pr cela je suis vraiment desolee!!! |
|
Revenir en haut de page |
|
|
|