function[x,h,f,df,d2f,dy,d2y]=tab1
a=input('Podaj wartość a: ')
b=input('Podaj wartość b: ')
n=input('Podaj wartość n: ')
g=input('Podaj wartość f(x): ','s')
j=input('Podaj wartość df: ','s')
k=input('Podaj wartość d2: ','s')
x=zeros(1,n+1);
f=zeros(1,n+1);
df=zeros(1,n+1);
d2f=zeros(1,n+1);
dy=zeros(1,n+1);
d2y=zeros(1,n+1);
h=(b-a)/n;
for i=1:n+1;
x(i)=a+(i-1)*h;
f(i)=eval(g);
df(i)=eval(j);
d2f(i)=eval(k);
end
h2=1/(2*h)
h22=1/(h*h)
dy(1)=h2*(-3*f(1)+4*f(2)-f(3));
d2y(1)=h22*(f(1)-2*f(2)+f(3));
for i=2:n
dy(i)=h2*(-f(i-1)+f(i+1));
d2y(i)=h22*(f(i-1)-2*f(i)+f(i+1));
end
dy(n+1)=h2*(f(n-1)-4*f(n)+3*f(n+1));
d2y(n+1)=h22*(f(n-1)-2*f(n)+f(n+1));
%hold on;
%plot(x,f)
%hold on;
%plot(x,df)
%hold on;
end