1340/CH7/EX7.10/7_10.sce


1
2
3
4
5
6
7
8
9

clc;
syms s a K;
G = K/(s^2+a*s+K);
disp(G,"Closed loop transfer function:");

DG = diff(G,K);
sta = (a/G)*diff(G,a);
sta = simple(sta);
disp(sta,"Sensitivity wrt parameter a = ")