Zadanie 1. Nie jestem pewien czy dobrze, bez tego co na niebiesko

function wczytywanie(a)

Zadanie. 2

function [x1,x2] = zad1rkw

a = input('podaj a = ');

b = 2*a; c=a;

delta = b*b-4*a*c; %taka delta zawsze bedzie = 0

if a~=0

if delta >0

x1 = (-b - sqrt(delta))/(2*a) ;

x2 = (-b + sqrt(delta))/(2*a);

fprintf ('x1 = %1.5f\n' , x1);

fprintf ('x2 = %1.5f\n' , x2);

elseif delta==0

x2 = -b/(2*a); %to zawsze bedzie rowne -1

x1 = x2;

fprintf ('x1 = x2 = %1.5f\n' , x1);

else

disp ('funkcja nie ma pierwiastkow rzeczywistych');

end

x = -2*a:0.01:2*a;

y = a*x.^2+b*x+c;

plot (x,y); grid on;

else

disp ('funkcja nie jest rownaniem kwadratowym')

end

return

Zadanie 3.

function [par] = zad2dat;

a = input ('Podaj a = ');

a3 = a; a2=-0.5*a; a1=-0.25*a; a0=-0.125*a;

par = [ a3 a2 a1 a0];

fp=fopen('par2dat.m','w');

fprintf(fp,'a3=%1.4f a2=%1.4f a1=%1.4f a0=%1.4f ',a3,a2,a1,a0);

fclose(fp);

return

function zad2obl

parametry= feval ('zad2dat');

mzer = roots (parametry);

fd = fopen ('zad2OUT.m' , 'w');

fprintf ( fd , ' x1 = %1.4f+(%1.4f)j x2 = %1.4f+%(1.4f)j x3 = %1.4f+%(1.4f)j',...

real(mzer(1)), imag(mzer(1)),real(mzer(2)), imag(mzer(2)), real(mzer(3)), imag(mzer(3)));

fclose(fd);

return

Zadanie 4

function [f]= fun4(x,a)
b=a/100;
f=([x(1)*x(2)+b*x(1)-b*x(2)-b*b;
x(1)*x(2)-2*b*x(2)-2*b*b]);
return

function [x]=zad4(a)
a=input('Podaj wartosc a=');
b=a/100;
if (b>=0.2 & b<=0.8)
x0=[2*b+0.1 b+0.1];
[x,fval]=fsolve(@fun4,x0)
else
disp('zle podane b’);
return
end
return

function [g]=grad4(x,a)
a=input('Podaj wartosc a=');
b=a/100;
if (b>=0.2 & b<=0.8)
x0=[2*b+0.1 b+0.1];
[x,fval,jacobian]=fsolve(@fun4,x0)

g=(jacobian)’
else
disp('zle podane b');
return
end
return