Notes

function La8Piorko_fun
%Autor Krystian Piorko
%nr Albumu 250612

% Funkcja frkw.m czytana dane z m-pliku,
% wyznacza miejsca zerowe za pomoca funkcji roots(),
% zapisuje wyniki w pliku tekstowych.
t0=clock; % poczatek obliczen
fd =fopen('frkwOUT.m','wt');
fprintf(fd, 'n Parametry rownania kwadratowego a,b,c nalezy przeczytac z pliku ');
fprintf(fd, 'n Date: %5d-%2d-%2d o''clock %2dh, %2dmin, %2.0fs',clock);
fprintf( 'n Date: %5d-%2d-%2d o''clock %2dh, %2dmin, %2.0fs',clock);
% wybor pliku z danymi
wdold=cd;
[fname,sciezka]=uigetfile('*dat.m','Wybierz dane do programu rkwewy, np. rkwewydat.m');
eval(['cd(''',sciezka,''')']); datafile=strtok(fname,'.');
fprintf(fd,'n Przeczytano dane z pliku: %s%s.mn',sciezka,datafile);
fprintf( 'n Przeczytano dane z pliku: %s%s.mn',sciezka,datafile);
% W(x) = a*x^2 + b*x + c - wielomian
[a,b,c,d,e]=feval(datafile); % czytanie danych z m-pliku za pomoca funkcji feval()
eval(['cd(''',wdold,''')']); % powrot do katalogu z programem
save 'La8Piorko_datBIN.nat' a b c d e -ascii

fprintf(fd,'n PARAMETRY ROWNANIA:');
fprintf(fd,'n a = %7.4f, b = %7.4f, c = %7.4f, d = %7.4f, e = %7.4f',a,b,c,d,e);
fprintf( 'n PARAMETRY ROWNANIA:');
fprintf( 'n a = %7.4f, b = %7.4f, c = %7.4f, d = %7.4f, e = %7.4f',a,b,c,d,e);
C=[a b c d e];

x=roots(C); % rozwiazanie wielomianu
rex=real(x); imx=imag(x);
fprintf(fd,'n ROZWIAZANIE ROWNANIA:');
fprintf(fd,'n Lp real(x) imag(x)');
fprintf( 'n ROZWIAZANIE ROWNANIA:');
fprintf( 'n Lp real(x) imag(x)');
nx=length(x);
t=-1:0.1:1;
y =a*t.^4 + b*t.^3 +c*t.^2 +d*t.^1 + e;
plot(t,y,'-b');
title('Krystian Piorko 250612 wykres');
grid on;
pause(4);
for i=1:nx
compass(rex(i), imx(i), '-b');
title('Krystian Piorko 250612 wykres typu compass');
hold on;
pause(4);
end


for i=1:2
fprintf(fd,'n%2d %15.6f %15.6f',i,rex(i),imx(i) );
fprintf( 'n%2d %15.6f %15.6f',i,rex(i),imx(i) );
end
et=etime(clock,t0); % czas obliczen
fprintf(fd, 'n Czas dzialania funkcji frkw() tobl = %.2f sekund !', et);
fprintf( 'n Czas dzialania funkcji frkw() tobl = %.2f sekund !', et);
fclose('all');
fprintf('nn K O N I E C obliczen programem frkw()');
return