Lab 5

% function [a,b] = lab5(p,L)
% if nargin<1
%     p=1;
%     L=1;
% elseif nargin<2
%     L=10;
% end
% a=100;
% x=linspace(-pi/2,pi/2,1000);
% y=L*(sqrt(cos(x*p)).*cos(a*x)+sqrt(abs(x))).*(4-x.^2).^(1/10);
% a=x;
% b=y;
% plot(x,y);

% function [g,h] = lab5(we)
% p=we.p;
% L=we.L;
% a=we.a;
% x=linspace(-pi/2,pi/2,1000);
% y=L*(sqrt(cos(x*p)).*cos(a*x)+sqrt(abs(x))).*(4-x.^2).^(1/10);
% g=x;
% h=y;
% plot(x,y);

% function f=lab5(x)
% p=1;
% a=100;
% f=(sqrt(cos(x*p)).*cos(a*x)+sqrt(abs(x))).*(4-x.^2).^(1/10);
% end
% xp=0;
% a=100;
% p=1;
% x1=(-pi/2,pi/2,1000);
% y=(sqrt(cos(x*p)).*cos(a*x)+sqrt(abs(x))).*(4-x.^2).^(1/10);
% plot(x1,y)
% hold on
% x0=fzero(@lab5,xp);
% plot(x0,0,'xm');

% ZLE