clc; clear all; close all;
format short g
t=(0:0.001:10)';

[Nob, ~ ]=size(t);

y1=8.7*sin(t)+0.001*cos(t);
y2=109.789*atan(t)-0.23*tanh(t);
y=[y1; y2];
fi=[sin(t), cos(t), zeros(Nob,1),  zeros(Nob,1); 
    zeros(Nob,1),  zeros(Nob,1), atan(t),  tanh(t)];
theta=mincuadm(y,fi,Nob,4,2);
theta