Question: How do I solve a differential equation in my question(given initial point)?

> u := 0.4e-1;
                              0.04
> k := 2.536427;
                            2.536427
> s := 9.466073;
                            9.466073
> R := 4;
                               4
>
> Eq1 := diff(v(t)*sin(theta(t)), t) = k-u*sin(theta(t))*(s*sin(`ϕ`(t))+(v(t)*cos(theta(t)))^2/R); Eq2 := diff(v(t)*cos(theta(t)), t) = s*cos(`ϕ`(t))-u*cos(theta(t))*(s*sin(`ϕ`(t))+(v(t)*cos(theta(t)))^2/R); Eq3 := diff(`ϕ`(t), t) = v(t)*cos(theta(t))/R;
/ d      \                                    / d          \  
|--- v(t)| sin(theta(t)) + v(t) cos(theta(t)) |--- theta(t)| =
\ dt     /                                    \ dt         /  

                                /                        
  2.536427 - 0.04 sin(theta(t)) |9.466073 sin(ϕ(t))
                                \                        

     1     2              2\
   + - v(t)  cos(theta(t)) |
     4                     /
/ d      \                                    / d          \  
|--- v(t)| cos(theta(t)) - v(t) sin(theta(t)) |--- theta(t)| =
\ dt     /                                    \ dt         /  

                                                 /            
  9.466073 cos(ϕ(t)) - 0.04 cos(theta(t)) |9.466073 sin(
                                                 \            

                 1     2              2\
  ϕ(t)) + - v(t)  cos(theta(t)) |
                 4                     /
              d                1                  
             --- ϕ(t) = - v(t) cos(theta(t))
              dt               4                  
> IC := v(0) = 5., theta(0) = .174533, `ϕ`(0) = 0.;
        v(0) = 5., theta(0) = 0.174533, ϕ(0) = 0.
> vthetafi := dsolve({Eq1, Eq2, Eq3, IC}, {v(t), theta(t), `ϕ`(t)}, type = numeric);
proc(x_rkf45)  ...  end;
> vthetafi(1.);
Error, (in vthetafi) cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

Please Wait...