Question: Polarplot warning for a semistable limit cycle

Hey, i'm trying do demonstrate that a nonlinear system has a semistable limit cycle but i get a warning at the plot command saying "Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct" and i dont understand it. So i wonder if someone here could help me? 

restart; with(PDEtools); with(plots);
eq1 := diff(x(t), t) = x(t)*(x(t)^2+y(t)^2-1)^2-y(t);
                                             2       
           d              /    2       2    \        
          --- x(t) = x(t) \x(t)  + y(t)  - 1/  - y(t)
           dt                                        
eq2 := diff(y(t), t) = y(t)*(x(t)^2+y(t)^2-1)^2+x(t);
                                             2       
           d              /    2       2    \        
          --- y(t) = y(t) \x(t)  + y(t)  - 1/  + x(t)
           dt                                        
tr := {x(t) = r(t)*cos(theta(t)), y(t) = r(t)*sin(theta(t))};
     {x(t) = r(t) cos(theta(t)), y(t) = r(t) sin(theta(t))}
eq1b := dchange(tr, x(t)*eq1+y(t)*eq2, [r(t), theta(t)], simplify);
              / d      \       2 /        4         2\
         r(t) |--- r(t)| = r(t)  \1 + r(t)  - 2 r(t) /
              \ dt     /                              
eq1b := expand(eq1b/r(t));
                d                    5         3
               --- r(t) = r(t) + r(t)  - 2 r(t) 
                dt                              
eq2b := dchange(tr, y(t)*eq1-x(t)*eq2, [r(t), theta(t)], simplify);
                      2 / d          \        2
                 -r(t)  |--- theta(t)| = -r(t) 
                        \ dt         /         
eq2b := simplify(eq2b/(-r(t)^2));
                         d              
                        --- theta(t) = 1
                         dt             
sol1 := dsolve({eq1b, r(0) = r[0]}, r(t));
          /      /  /     2  \                 
          |      |  | r[0]   |          2     2
r(t) = exp|RootOf|ln|--------| (exp(_Z))  r[0] 
          \      \  \r[0] - 1/                 

                           2     2
   - ln(r[0] + 1) (exp(_Z))  r[0] 

       /             2\                                       
       |(exp(_Z) - 1) |          2     2            2        2
   - ln|--------------| (exp(_Z))  r[0]  + (exp(_Z))  _Z r[0] 
       \ exp(_Z) - 2  /                                       

                                /     2  \              
                2     2         | r[0]   |             2
   + 2 (exp(_Z))  r[0]  t - 2 ln|--------| exp(_Z) r[0] 
                                \r[0] - 1/              

                                2
   + 2 ln(r[0] + 1) exp(_Z) r[0] 

         /             2\                                   
         |(exp(_Z) - 1) |             2                    2
   + 2 ln|--------------| exp(_Z) r[0]  - 2 exp(_Z) _Z r[0] 
         \ exp(_Z) - 2  /                                   

                           /     2  \           
                   2       | r[0]   |          2
   - 4 exp(_Z) r[0]  t - ln|--------| (exp(_Z)) 
                           \r[0] - 1/           

                                 /             2\           
                           2     |(exp(_Z) - 1) |          2
   + ln(r[0] + 1) (exp(_Z))  + ln|--------------| (exp(_Z)) 
                                 \ exp(_Z) - 2  /           

                                          /     2  \        
              2                   2       | r[0]   |        
   - (exp(_Z))  _Z - 2 t (exp(_Z))  + 2 ln|--------| exp(_Z)
                                          \r[0] - 1/        

                                  /             2\        
                                  |(exp(_Z) - 1) |        
   - 2 ln(r[0] + 1) exp(_Z) - 2 ln|--------------| exp(_Z)
                                  \ exp(_Z) - 2  /        

              2                                    2            
   - (exp(_Z))  + 2 _Z exp(_Z) + 4 t exp(_Z) + r[0]  + 2 exp(_Z)

      \\    
      ||    
   - 1|| - 1
      //    
sol1 := simplify(sol1);
          /      /   /     2  \                
          |      |   | r[0]   |               2
r(t) = exp|RootOf|-ln|--------| exp(2 _Z) r[0] 
          \      \   \r[0] - 1/                

                                2
   + ln(r[0] + 1) exp(2 _Z) r[0] 

       /             2\                                     
       |(exp(_Z) - 1) |               2                    2
   + ln|--------------| exp(2 _Z) r[0]  - exp(2 _Z) _Z r[0] 
       \ exp(_Z) - 2  /                                     

                               /     2  \              
                     2         | r[0]   |             2
   - 2 exp(2 _Z) r[0]  t + 2 ln|--------| exp(_Z) r[0] 
                               \r[0] - 1/              

                                2
   - 2 ln(r[0] + 1) exp(_Z) r[0] 

         /             2\                                   
         |(exp(_Z) - 1) |             2                    2
   - 2 ln|--------------| exp(_Z) r[0]  + 2 exp(_Z) _Z r[0] 
         \ exp(_Z) - 2  /                                   

                           /     2  \          
                   2       | r[0]   |          
   + 4 exp(_Z) r[0]  t + ln|--------| exp(2 _Z)
                           \r[0] - 1/          

                                /             2\          
                                |(exp(_Z) - 1) |          
   - ln(r[0] + 1) exp(2 _Z) - ln|--------------| exp(2 _Z)
                                \ exp(_Z) - 2  /          

                                        /     2  \        
                                        | r[0]   |        
   + exp(2 _Z) _Z + 2 t exp(2 _Z) - 2 ln|--------| exp(_Z)
                                        \r[0] - 1/        

                                  /             2\        
                                  |(exp(_Z) - 1) |        
   + 2 ln(r[0] + 1) exp(_Z) + 2 ln|--------------| exp(_Z)
                                  \ exp(_Z) - 2  /        

                                                  2            
   + exp(2 _Z) - 2 _Z exp(_Z) - 4 t exp(_Z) - r[0]  - 2 exp(_Z)

      \\    
      ||    
   + 1|| - 1
      //    
sol2 := dsolve({eq2b, theta(0) = theta[0]}, theta(t));
                    theta(t) = t + theta[0]
theta[0] := (1/4)*Pi;
                              1   
                              - Pi
                              4   
plot1 := polarplot([subs(r[0] = .1, rhs(sol1)), rhs(sol2), t = 0 .. 10], color = red);
Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct
plot2 := polarplot([subs(r[0] = 2, rhs(sol1)), rhs(sol2), t = 0 .. 10], color = blue);
Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct
display({plot1, plot2}, scaling = constrained, tickmarks = [4, 3], view = [-2 .. 2, -2 .. 2]);