function [ms_i]  = myfun2(gp, kp, tp, theta, ms_s, tc, type)

[gc kc_i ti_i] = ituning(kp, tp, theta, tc, type)
 ms_i = ms_s -  ms(gp, kc_i, ti_i);
 
end
