Laboratorinio darbo Nr.2 ataskaita

x0=input('iveskite x0 ===>');

xn=input('iveskite xn ===>');

xh=input('iveskite xh ===>');

a=input('iveskite a ===>');

nx=0;

for xv =x0: xh: xn

nx=nx+1;

X(nx)=xv;

if( xv<=10)

Y(nx)= [a]/[xv]+sqrt(xv^2);

elseif (xv <=20)

Y(nx)=2*a^2/(xv^2-2);

else

Y(nx)=(3/xv^2)-4;

end

end

disp('+----- ----- ----------+----- ----- --------+');

disp('| X | Y |');

disp('+----- ----- ----------+----- ----- --------+');

for ix=1:nx

strx=num2str(X(ix),'%10.4f');

stry=num2str(Y(ix),'%10.4f');

disp( ['|',blanks(19-length(strx)), strx, ' |' , blanks(17-length(stry)), stry, ' |']);

end

disp('+----- ----- ----------+----- ----- --------+');

str=num2str(x0,'%10.4f');

disp(['| x0=',blanks(35-length(str)), str, '|']);

str=num2str(xn,'%10.4f');

disp(['| xn=',blanks(35-length(str)), str, '|']);

str=num2str(xh,'%10.4f');

disp(['| xh=',blanks(35-length(str)), str, '|']);

str=num2str(a,'%10.4f');

disp(['| a=',blanks(36-length(str)), str, '|']);

disp('+-------- ----- ------ ------+');

plot(X, Y);