Question: How can i get datas of implicit function

Equação26_12.mw

Download Equação26_12.mw

Hello!

How can i get datas and graphs of implicit function. I tried...

with(DEtools, odeadvisor);

ode := diff(y(x), x) = (-0.1*(y(x) - 0.7)^3 - 0.8*y(x)*(y(x) - 0.7)^2)/(0.7^2 - 2*0.7*y(x) + y(x)^2 - 0.7*0.2*0.8);
                                    3                        2
         d         -0.1 (y(x) - 0.7)  - 0.8 y(x) (y(x) - 0.7) 
 ode := --- y(x) = -------------------------------------------
         dx                                        2          
                            0.378 - 1.4 y(x) + y(x)           
odeadvisor(ode);
          /                            3                        2\
          | d         -0.1 (y(x) - 0.7)  - 0.8 y(x) (y(x) - 0.7) |
odeadvisor|--- y(x) = -------------------------------------------|
          | dx                                        2          |
          \                    0.378 - 1.4 y(x) + y(x)           /
ans := dsolve(ode, implicit);
                2        9                    199                
ans := x + ----------- + -- ln(10 y(x) - 7) + --- ln(90 y(x) - 7)
           10 y(x) - 7   28                   252                

   + _C1 = 0
ic1 := y(0) = 0;
                        ic1 := y(0) = 0
sol := dsolve([ode, ic1], implicit);
                2        9                    199                
sol := x + ----------- + -- ln(10 y(x) - 7) + --- ln(90 y(x) - 7)
           10 y(x) - 7   28                   252                

     2   10         10         
   + - - -- ln(7) - -- I Pi = 0
     7   9          9          
subs(x = 1, sol);
   9        2        9                    199                
   - + ----------- + -- ln(10 y(1) - 7) + --- ln(90 y(1) - 7)
   7   10 y(1) - 7   28                   252                

        10         10         
      - -- ln(7) - -- I Pi = 0
        9          9          
solve(9/7 + 2/(10*y(1) - 7) + (9*ln(10*y(1) - 7))/28 + (199*ln(90*y(1) - 7))/252 - (10*ln(7))/9 - 10*I*Pi/9 = 0, y(1));
1     /      /                                                  
-- exp|RootOf|-280 I Pi exp(_Z) - 280 ln(7) exp(_Z) + 15680 I Pi
90    \      \                                                  

          /1           56\                                       
   + 81 ln|- exp(_Z) - --| exp(_Z) + 199 _Z exp(_Z) + 15680 ln(7)
          \9           9 /                                       

                          /1           56\                   \\
   + 324 exp(_Z) - 4536 ln|- exp(_Z) - --| - 11144 _Z - 13608||
                          \9           9 /                   //

     7 
   + --
     90
implicitplot(sol, x = 0 .. 10, y = 0 .. 10);
                /         2        9                 
    implicitplot|x + ----------- + -- ln(10 y(x) - 7)
                \    10 y(x) - 7   28                

         199                   2   10         10           
       + --- ln(90 y(x) - 7) + - - -- ln(7) - -- I Pi = 0, 
         252                   7   9          9            

                              \
      x = 0 .. 10, y = 0 .. 10|
                              /
but I wasn't successful.

Thanks