(* ::Package:: *)

(************************************************************************)
(* This file was generated automatically by the Mathematica front end.  *)
(* It contains Initialization cells from a Notebook file, which         *)
(* typically will have the same name as this file except ending in      *)
(* ".nb" instead of ".m".                                               *)
(*                                                                      *)
(* This file is intended to be loaded into the Mathematica kernel using *)
(* the package loading commands Get or Needs.  Doing so is equivalent   *)
(* to using the Evaluate Initialization Cells menu command in the front *)
(* end.                                                                 *)
(*                                                                      *)
(* DO NOT EDIT THIS FILE.  This entire file is regenerated              *)
(* automatically each time the parent Notebook file is saved in the     *)
(* Mathematica front end.  Any changes you make to this file will be    *)
(* overwritten.                                                         *)
(************************************************************************)



SetDirectory["C:\\drorbn\\AcademicPensieve\\2015-09\\K34"]


T=\[Pi](1-Cos[t/2]);
R=-Cos[t/2];
A=3.5;
b=-2.6;
a=1.95;
r=0.2*0.4*(9*R^4+1)*Sqrt[1-R^2];
circ= {0,Sin[\[Theta]]*r,Cos[\[Theta]]*r};
dx=-0.5545454545*Sin[T]^2+6.1*(10+Cos[T])/11*Cos[T]+1.3*Sin[T/2];
dy=4.875*Cos[1.95*T+0.05*\[Pi]];
n=Sqrt[dx^2+dy^2];
dx=dx/n;
dy=dy/n;
rotMat=({
 {dx, -dy, 0},
 {dy, dx, 0},
 {0, 0, 0.5*(1+2*(t/\[Pi]-1)^2)}
});
result=rotMat.circ;
coords={
0.2*((A-b)*((10.0+Cos[T])/11.0)*Sin[T]+b*Cos[T/2])+result[[1]],
0.5*Sin[a*T-(a-2)*\[Pi]]+result[[2]],
result[[3]],
If[t<7/8*\[Pi], (1-Cos[t*8/7])/2, 0.7*(Cos[8*t-7*\[Pi]])+0.3]
};



