clf()

x=linspace(-5,5,51);

y=1/(1+x^2);

plot(x,y,'b-')

xgrid()

legend(['1/(1+x^2)'])

xtitle('','oś-x','oś-y')

scf(1)

clf(1)

x=linspace(-1,1,100)

y=linspace(-1,1,100)

x1=x'*ones(y)

y1=ones(x')*y

z=x1.^2+y1.^2

plot3d1(x,y,z, theta=70, alpha=70, flag=[7 2 4])

w1=poly([0:1:3],'x','c')

w2=poly([1,-2,3,2},'x','c')

w3=poly([0,-1,0,-1,-1],'x','c')

b) p1=derivat(w1) p2=derivat(p1)

c) horner(w1,2)

d) degree(w1)

function y=f3(x)

y=x^3+0.5*x^2-9.25*x+4.375

endfunction

x=(-5:0.01:4);

xgrid();

plot2d(x,f3(x))

x1=fsolve(-4,f3)

x2=fsolve(0,f3)

x3=fsolve(3,f3)

disp(x1)

disp(x2)

disp(x3)

function [xd]=f1(t, x)

xd(1)=x(2);

xd(2)=-a*x(2)-b*x(1);

endfunction

a=input('podaj wartość a: ');

b=input('podaj wartość b: ');

x01=input('podaj wartość x w chwili początkowej: ');

x02=input('podaj wartość x* w chwili początkowej: ');

tk=input('podaj czas symulacji: ');

t0=0;x0=[x01;x02]

T=linspace(t0,tk);

[x]=ode(x0,t0,T,f1);

clf();

plot2d(T',[x(1,:)',x(2,:)'],leg="x@x*")

xtitle("rozwiązanie graficzne")

xgrid()