function [b] = f3(x) 

% Funktion aus Aufgabe 33 (c)

n=length(x);
b=zeros(n,1);

b(1)=-400*x(1)*(x(2)-x(1)^2)+2*x(1)-2;
b(2)=400*(x(2)-x(1)^2);