滷蛋與貢丸

clear
p0=0.78;
p1=0.7;
Tol=10^(-10);
N0=10;
stop=0;
i=2;
g0=cos(p0)-p0;
g1=cos(p1)-p1;
while(i<N0 & stop==0)
    temp=g1*(p1-p0)/(g1-g0);
    p=p1-temp;
    tempabs(i,1:3)=[i p abs(p-p1)];
    if abs(p-p1)<Tol
        stop=1
        p
    end
    i=i+1;
    p0=p1;
    g0=g1;
    p1=p;
    g1=cos(p)-p;
end

clear
p0=5.5;
Tol=10^(-10);
N0=10;
stop=0;

temp1=cos(p0)-p0;
temp2=-sin(p0)-1;
p=p0-temp1/temp2;

clear format long e
p0=1.5
tol=0.000000000001;
N0=100;
allp=p0;

clear format long e
p0=1.5
tol=0.000000000001;
N0=100;
allp=p0;

下面是函數那一部份的寫法： %p=p0-p0^3-4*p0^2+10;  p=0.5*(10-p0^3)^0.5; % temp1=p0^3+4*p0^2-10;
% temp2=3*p0^2+8*p0;
% p=p0-temp1/temp2;