Forum grupy 6 WIL

Gość

function [f,d2f]=fun(x)

f=x^4;
d2f=12*x^2;

endfunction

function [x,f,d2f,h]=tabela(a,b,n)

x=zeros(1,n+1);
f=zeros(1,n+1);
d2f=zeros(1,n+1);

h=(b-a)/n;

for i=1:n+1

x(i)=a+(i-1)*h;
[f(i),d2f(i)]=fun(x(i));

end

endfunction

function [d2y]=poch(h,f,n)

d2y=zeros(1,n+1);

h22=1/(h*h);

for i=3:n+1

d2y(i)=h22*(f(i-2)-2*f(i-1)+f(i));

end

endfunction

a=-3;
b=abs(a);
n=15;

gset grid

[x,f,d2f,h]=tabela(a,b,n);
[d2y]=poch(h,f,n);

gset grid;
xlabel('X');
ylabel('Y');
title('Druga pochodna funkcji');

printf("wykres funkcji f(x) purpurowy\n");
plot(x,f,'m');
hold on;
pause;

printf("druga pochodna scisla d2f zielony\n");
plot(x,d2f,'g');
hold on;
pause;

printf("druga pochodna numeryczna d2y niebieski\n");
plot(x(3:n+1),d2y(3:n+1),'b');
hold on;
pause;

[/b]

Gość

Huge cocks, wet pussies, hardcore action!
[link widoczny dla zalogowanych]