222.html

Manual Maple 2.2.2.mws

Comandos Úteis e Simples

Section 2.2: Resolvendo Equações

Session 2.2.2: fsolve

> restart;

> Digits := 20;

> p := x^5-x+1;

> r := fsolve( p, x );

> subs( x=r, p );

> pseudozeros := fsolve( p, x, complex );

> map( evalf[1],[seq( abs(eval( p, x=z)), z=pseudozeros )]);

> Digits := trunc(evalhf(Digits));

> X := fsolve( x*tan(x) - 1, x );

> X*tan(X)-1;

> plot( [x*tan(x),1], x=-10..10, y=-5..5, discont=true, colour=black );

[Maple Plot]

> X := fsolve( x*tan(x)-1, x=3);

> X*tan(X);

> X := fsolve( x*tan(x)-1, x=6);

> X*tan(X);

> Digits := 20;

> r := fsolve( x*tan(x)-1, x=1000*Pi );

> r*tan(r);

> r - evalf(1000*Pi);

> 1/%;

> assume(n,integer);

> T := subs( x=(n*Pi+delta), x*tan(x) );

> T := series( T, delta );

> solve( T=t, delta );

series(1/(n*Pi)*t+(-1/(n^3*Pi^3))*t^2+(-1/3*(-6+n^2*Pi^2)/n^5/Pi^5)*t^3+1/3*(-15+4*n^2*Pi^2)/n^7/Pi^7*t^4+1/5*(70-25*n^2*Pi^2+n^4*Pi^4)/n^9/Pi^9*t^5+O(t^6),t,6)

> asympt( eval( convert( %, polynom ), t=1 ), n, 7 );

1/(n*Pi)-4/3*1/(n^3*Pi^3)+O(1/(n^5))

>

> restart;

> sys := {x^2+y^2-1,25*x*y-12};

$sys := {x^2+y^2-1, 25*x*y-12}$

> fsolve( sys, {x,y} );

> fsolve( sys, {x=-1..0,y=0..1} );

$fsolve({x^2+y^2-1, 25*x*y-12},{x, y},{x = -1 .. 0, y = 0 .. 1})$

> fsolve( sys, {x=0.2..1.2, y=0.2..1.2} );

> fsolve( sys, {x=-1.2..-0.2, y=-1.2..-0.2} );

>