lab 6

#a = sym('a');
#b = sym('b')
syms a b;
x = a+3*a+2*b^2-4*b+b^2;
#y = a^2-b^2;
#disp(factor(y));
f = @ (x) sin(x)+2*x;
disp(f(5));
disp(f(0.25));
disp(feval(f,pi/3));
disp(f(4:2:24));
f1 = (diff(f(x)));
f2 =(diff(diff(f(x))));
l1 =(int(f(x)));
l2 =(int(f(x),0,1));
l3 = (quad(f,0,1));
ezplot(f,[0,10]);

syms x;
c = @(x) [2*cos(x),2*sin(x)];
c1 = diff(c(x));
disp(c1);
c2 = diff(diff(c(x)));
disp(c2);
norma = norm(c1);
disp(norma);