bank.sce~

source code which tells whether customer has to given loan or not depending on the details of the cu

function []=bank()
    printf("\n Enter File : ");
    fn=scanf("%s");
    
  
    fp=mopen(fn,'r');
      n=mfscanf(fp,'%d');
      incmin=mfscanf(fp,'%f');
      incmax=mfscanf(fp,'%f');
      savmin=mfscanf(fp,'%f');
      savmax=mfscanf(fp,'%f');

    for i=1:n,
      x(i)=mfscanf(fp,'%f');
      y(i)=mfscanf(fp,'%f');
      r(i)=mfscanf(fp,'%d');
    end
    mclose(fp);
    
    xsuml=0.0;xmeanl=0.0;xvarl=0.0
    xsumo=0.0;xmeano=0.0;xvaro=0.0;    
    ysuml=0.0;ymeanl=0.0;yvarl=0.0;
    ysumo=0.0;ymeano=0.0;yvaro=0.0;

    ////////////////////////////////////////////////////////////////////////////
     
     nl=0; no=0;     // Caliculations for mean
     
     for i=1:n,
      if r(i)==1 then
      nl=nl+1;
      xsuml=xsuml+x(i);
      ysuml=ysuml+y(i);
      end
      if r(i)==0 then
      no=no+1;
      xsumo=xsumo+x(i);
      ysumo=ysumo+y(i);
      end
     end
  

      xmeanl=xsuml/nl;
      xmeano=xsumo/no;

   
                         // Caliculations for Probabilities p(C=1) & P(C=0) 
      xpcl=nl/n;
      xpco=no/n;
      ypcl=nl/n;
      ypco=no/n;                      
                        // Caliculations for variance 
      for i=1:n,      
         if r(i)==1 then
            xvarl=xvarl+(abs(x(i)-xmeanl)^2);
            yvarl=yvarl+(abs(y(i)-ymeanl)^2);
          end,
          if r(i)==0 then
            xvaro=xvaro+(abs(x(i)-xmeano)^2);
            yvaro=yvaro+(abs(y(i)-ymeano)^2);
          end,
      end;
      
     xvarl=real(xvarl/nl);
     yvarl=real(yvarl/nl);
      xvaro=real(xvaro/no);
      yvaro=real(yvaro/no);
      
      
      
      for i=1:n,
        if r(i)==1 then
           gxl(i)=-log(sqrt(xvarl))-((x(i)-xmeanl)^2)/(2*xvarl)+log(xpcl);
           gyl(i)=-log(sqrt(yvarl))-((y(i)-ymeanl)^2)/(2*yvarl)+log(ypcl);
        end
        if r(i)==0 then
           gxo(i)=-log(sqrt(xvaro))-((x(i)-xmeano)^2)/(2*xvaro)+log(xpco);
           gyo(i)=-log(sqrt(yvaro))-((y(i)-ymeano)^2)/(2*yvaro)+log(ypco);
           
        end
     end
     
     minx1=9999;miny1=9999;minx0=9999;miny0=9999;
               
       for i=1:n,
        if r(i)==1 then
          if minx1>gxl(i) then  minx1=gxl(i); end
          if miny1>gyl(i) then  miny1=gyl(i); end
        end
         if r(i)==0 then
           if minx0>gxo(i) then minx0=gxo(i); end
           if miny0>gyo(i) then  miny0=gyo(i); end
         end
       end
       
       if minx1<minx0 then gx=minx1; else gx=minx0; end
       if miny1<miny0 then gy=miny1; else gy=miny0; end
       
                          //    we know the x, y values (threshhold values)

      printf("\n Enter Test input data ");
      printf("\n Enter Income and Savings: ");
      inc=scanf("%f");
      sav=scanf("%f");
      
      //////////////////////////////////////////////////////////////////////
      
      dgxl=-log(sqrt(xvarl))-((inc-xmeanl)^2)/(2*xvarl)+log(xpcl);
      dgyl=-log(sqrt(yvarl))-((sav-ymeanl)^2)/(2*yvarl)+log(ypcl);
      dgxo=-log(sqrt(xvaro))-((inc-xmeano)^2)/(2*xvaro)+log(xpco);
      dgyo=-log(sqrt(yvaro))-((sav-ymeano)^2)/(2*yvaro)+log(ypco);
      
//      printf("\n %f %f %f %f",dgxl,dgxo,dgyl,dgyo);
//      printf("\n %f  %f ",gx,gy);
      
      if dgxl<dgxo then
      dx=dgxl;
      else
      dx=dgxo;
      end
      if dgyl<dgyo then
      dy=dgy1;
      else
      dy=dgyo;
      end
       
      
      ///////////////////////////////////////////////////////////////////

          
      if(dx>gx & dy>gy) then
       printf("Accepted..");
      end
      if(dx<gx | dy <gy) then
        printf("Rejected..");
      end  
                              
endfunction;