/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 12.04.0 ] */

/* [wxMaxima: title   start ]
Oscilator sa faznim pomerajem,
izvodjenje
   [wxMaxima: title   end   ] */

/* [wxMaxima: comment start ]
sistem jednacina po potencijalima cvorova
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
e1: v1=vin;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e2: -Y*v1 + (1/R+2*Y)*v2 - Y*v3 = 0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e3: -v2*Y + v3*(2*Y+1/R) - Y*vout = 0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
e4: -v3*Y + vout*(1/R + Y) = 0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
eliminacija
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sol: eliminate([e1, e2, e3, e4], [v1, v2, v3]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sol: sol[1];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sol: subst(s*C, Y, sol);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
funkcija prenosa
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
tf: solve(sol, vout);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
tf: tf[1];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
FUNKCIJA PRENOSA

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
tf: rhs(tf)/vin;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
normalizacija,

NORMALIZOVANA FUNKCIJA PRENOSA

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
tf: subst(s/(R*C), s, tf);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
prelaz na jw osu
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
hjw: subst(%i*w, s, tf);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
declare(w, real);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
"fazna rezonansa"
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
ihjw: imagpart(hjw);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
wr: solve(ihjw, w);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
FREKVENCIJA OSCILOVANJA, frekvencija "fazne rezonanse"

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
wr: rhs(wr[2]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
hjwr: subst(wr, w, hjw);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
SLABLJENJE BETA MREZE NA FREKVENCIJI OSCILOVANJA

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
hjwr: ratsimp(hjwr);
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$