Pagina 1 van 1

Oefening 17.5: visCurt

Geplaatst: 23 feb 2008, 19:59
door Aushim

Code: Selecteer alles

function z = visCurt(a, guess)
% constanten
N = 50;
step = 0.1;
fout=0.0000001;
%functiegrafiek plotten
X = guess + step*(-N:N);
Y = X.^3-a;
plot(X,Y);

hold on
  %begingok
  z=guess;
  while abs(a-z^3)>fout
    %functiewaarde berekenen
    f = z^3 - a;
    %afgeleide berekenen
    Df = 3*z^2;
    %rico = afgeleide
    m = Df;
    %punt op de grafiek berekenen
    x0 = z;
    y0 = f;
    %punt op de raaklijn
    x1 = (m*x0-y0)/m;
    y1 = 0;
    %raaklijn tekenen
    plot([x0 x1],[y0 y1]);
    %volgend punt berekenen
    z=z-(z^3-a)/(3*z^2);
  end
hold off


Oplossing gehaald uit de website van Egon Geerardyn.

Directe link: http://homepages.vub.ac.be/~egeerard/fi ... MatLab.zip