Wednesday, November 12, 2014

Matlab Program for finding positive root of a Non-Linear Function using Regula-Falsi Method

k=input('Enter the function =  ', 's');
f= inline(k);
a=0; b=1;
for i=1:inf
    if f(a)*f(b)<0
        break
    else
        a=a+1;
        b=b+1;
    end
end
for i=2:inf
   x(i)= (a*f(b)-b*f(a))/(f(b)-f(a));
   if f(x(i))*f(a)<0
       b=x(i);
   else
       a=x(i);
   end
   if abs(x(i)-x(i-1))<=0
      Iteration=i-1
      Root=x(i)
      break
   end
end

Tuesday, November 4, 2014

Matlab Program for finding positive root of a Non-Linear Function using Newton Rapson Method

a=input('Enter the function =  ', 's');
f= inline(a);
h= diff(a);
d= inline(h);
x(1)=0;
for k=1:inf
    if f(x(1))>0
      for i=2:inf
      x(i)= x(i-1)-((f(x(i-1)))/(d(x(i-1))));
       if abs(x(i)-x(i-1))<=0.001 && x(i)>0  
       Root=x(i)
       break 
       end
       end
       break
     else
      x(1)=x(1)+1;
    end
end

Monday, November 3, 2014

Matlab Program for finding positive root of a Non-Linear Funtion using Secant Method

a=input('Enter the function =  ', 's');
f= inline(a);
x(1)=0; x(2)=1;iteration=0;
for k=1:inf
 if f(x(1))*f(x(2))<0
  for i=3:inf
  x(i)=(x(i-2)*(f(x(i-1)))-x(i-1)*(f(x(i-2))))/(f(x(i-1))-f(x(i-2)));
  iteration=iteration+1;
    if abs(x(i)-x(i-1))<=0.001 && x(i)>0
    root=x(i)
    iteration=iteration
    break
    end
  end
  break
 else
 x(1)=x(1)+1;
 x(2)=x(2)+1;
 end
end