clear all;
clc;
% Menuliskan variabel
syms w(t) theta a phi q s
assume([a q] > 0)
% Menuliskan persamaan diferensial
dw = diff(w,t);
eqn = dw == sin (theta) + a*cos(phi) - w/q
% Menuliskan kondisi awal
cond1 = w(0) == 0
cond2 = dw(0) == 0
% Transformasi Laplace
eqnLT = laplace(eqn,t,s)
syms w_LT 
eqnLT = subs(eqnLT,[laplace(w,t,s)], w_LT)
eqns = eqnLT;
vars = [w_LT];
[w_LT] = solve(eqns,vars)
% Transformasi Laplace balik
wsol = ilaplace(w_LT,s,t);
wsol = simplify(wsol);
% Substitusi Variabel dengan Nilai
vars = [a q theta phi w(0) dw(0)];
values = [1 2 -1*pi/2 pi/6 0 0];
wsol = subs(wsol,vars,values)
wsol = simplify(wsol);
pretty (wsol);
% Menggambar grafik
fplot(wsol, [0 10])                      
title('Kecepatan sudut')
ylabel('w(t)')
xlabel('t')