Originally Posted by

**anizatie** plz help me, im newbie in matlab.. ive tried so many ways to make sure it is correct but then still cannot work out..

i hv to do some iteration, and plot as well (save the answer in other place since the answer quite long)

function a = lukis(a0,s,n)

lamda = 1.0*(10e-6);

m = 2;

a(1) = a0;

for i = 1:n

a(i+1) = a(i) + lamda * a(i).^(m/2) * s.^m

if(a(i+1)>1)

break;

end

end

plot(a(i+1))

hold on

save('lukis', 'a(i+1)')

to run in c/w : lukis(1.0*10e-10,0.5,1000)

thanks.

Code:

function a = lukis(a0,s,n)
lamda = 1.0*(10e-6);
m = 2;
a=zeros(1,n);
a(1) = a0;
for i = 1:n-1
a(i+1) = a(i) + lamda * a(i).^(m/2) * s.^m
if(a(i+1)>1)
break;
end
end
plot(a(i+1))
hold on
save('lukis', 'a')

CB