วันอังคารที่ 23 มกราคม พ.ศ. 2561

Sl 15



Sl 11



Sl 14



Sl 13




Sl 12



Sl 10


Sl 9


Sl 8




วันพุธที่ 17 มกราคม พ.ศ. 2561

วันอังคารที่ 16 มกราคม พ.ศ. 2561

Sl 2


Sl 1

http://www.scilab.org/en/download/latest
SciLab 6.0.0 for Windows

Ex 25

ตอบ 
มีกองสต็อกเท่ากันที่ 51661.13 ตัน
ที่เวลา 57.57ชั่วโมง

Ex 24

3*x^2+3*y^2+14, 5x+2*y^3

Google แก้สมการไม่ได้

Ex 23





Ex 22


วันจันทร์ที่ 15 มกราคม พ.ศ. 2561

Ex 21


Ex 20

จุดตัดใน Q1 (1.388, 0.778)
จุดตัดใน Q3 (-0.720, -3.446)

สิ่งที่ อ.สอนในคาบ

การค้นหาไฟล์ใน Google
พิมพ์
Blasting ext pdf (หาไฟล์ pdf แบบระบุชื่อ)
* inurl:www...... ext:pdf  (หาไฟล์ pdf แบบไม่ระบุชื่อจากเว็บ)

กูเกิ้ลวาดกราฟ
ไม่ต้องพิมพ์ Y=.....
3x^2 = 3*x^2

Scilab
Clear(a) ......เคลียเฉพาะ a
Clear(all) = clc .........เคลียทุกตัวแปร


exp a ...........exponantial a 
save a........คือเก็บ a ไว้ แม้เคลียแล้ว แต่ยังสามาใช้คำสั่ง load a เรียกค่ากลับมาได้
ทำแนวตั้ง แนวนอน



x=-10.0:0.5:10.0
xf=x'
y=sin(x*%pi/10)
plot(x,y)
help plot
u=[-1,2,3]
v=[6,-10,0]
u+v
u-v
u*v
u*v'
u'*v
u(3)
u(2)+v(1)
a=1;b=2;c=3;r=[a,b,c]
s={-1,2}
letters={'a','b','c','d'}
A=[1,2,3;4,5,6;1,-1,0]
B=[1,1,1,
2,3,-1,
5,7,-2]
u=[1,3,-5];v=[4,2,3];w=[-1,0,1];
C=[u;v;w]
D=[u'v'w']
r=[u,v,w]
A+B
C-D
A*B
B*A
C*u
D*v'
rank(A)
inv(A)
cond(B)
det(C)
A*inv(A)
inv(B)*B
spec(A)
trace(C)
-----------------------
เว้นวรรคแทน , ได้ แต่โอกาสพลาดง่ายไม่แนะนำ 
x={.5 .7 .9 1.3 1.7 1.8]
y=[.1 .2 .75 1.5 2.1 2.4]
วาดกราฟ 2D เป็นจุด
(-1 เป็น +)
(-2 เป็น X)
(-3 เป็น วงกลมมีบวก) 
-4 เป็นหน้าจั่วทึบ
-5 เป็นหน้าจั่วโปร่ง
-6 เป็นสามเหลี่ยมโปร่งชี้ขึ้น
-7 เป็นสามเหลี่ยมโปร่งชี้ลง
-8 เป็นหน้าจั่วโปร่งมีบวก
-9 เป็นวงกลมโปร่ง
-10 เป็นดอกจันมีบวก
-11 เป็นสี่เหลี่ยมโปร่ง
plot2d(x,y,style=-1)
ไม่ต้องปิดกราฟแรก คำสั่งต่อมาจะได้เส้นเพิ่ม ถ้าปิดไปแล้วจะได้แต่เส้น
plot2d(x,y)
พลอตกราฟรูปหัวใจ
x1=linspace(-1,0,50);
y1=-x1+sqrt(3-3*x1.^2);
y2=-x1-sqrt(3-3*x1.^2);
title('heart shape curve');
plot(x1,y1,'r');
plot(-x1,y1,'r');
plot(-x1,y2,'r');
plot(x1,y2,'r');
------------------------------
n=[3 5 9 17 33 65]; s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4]; plot2d(n,s,style=-1) plot2d(n,s,style=5)

พิมพ์ใน Scinote
n=[3 5 9 17 33 65];
s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=4,logflag='ll')
plot2d(n,s,style=4,logflag='nl')
plot2d(n,s,style=4,logflag='ln')
plot2d(n,s,style=-1) ใส่จุด
ln คือเซมิล็อก y เป็น log แกน x เป็น normal
nl คือเซมิล็อก x เป็น log แกน y เป็น normal
xgrid(10); ใส่เส้นกริดในกราฟ

clear
n=10
p=100
i=5
for j=1:n
    f(j)=p*(1+(i/100))^j;
    printf('\nFuture value for each period.%3f\n',f(j))
        j=j+1;
    end
 
clear
n=input('input number of period...');
p=input('input Present value...');
i=input('input interest rate in percent...');
for j=1:n
    f(j)=p*(1+(i/100))^j;
    printf('\nFuture value for each period.%3f\n',f(j))
        j=j+1;
    end
clear
n=input("จำนวนปี...");
p=input("เงินฝาก...");
i=input("อัตราดอกเบี้ยต่อปี...");
for j=1:n
    f(j)=p*(1+(i/100))^j;
    printf('\nจำนวนเงินทั้งหมดที่ได้รับในแต่ละปี..%3f\n',f(j))
        j=j+1;

    end
กดรันใน Scinote แล้วพิมพ์ n=10,P=100,i=5

clear
a=input("เงินเดือน...");
b=input("จำนวนปี...");
i=input("อัตราดอกเบี้ยต่อปี...");
for j=1:b
    p=a*12*b
    f(j)=p*(1+(i/100))^j;
    printf('\nจำนวนเงินทั้งหมดที่ได้รับในแต่ละปี.%d..%3f\n',j,f(j))
        j=j+1;
    end

การคำนวณแฟคทอเรียล
function fact=factorial(n)
fact=1
for k=1:n
    fact =k*fact
end
endfunction
Run ใน Scinote แล้วพิมพ์
factorial(5) ใน Scilab ตามใจใส่เลขไร
เปลี่ยนชื่อฟังก์ชั่นได้นะ 
function fact=Jinna(n)
fact=1
for k=1:n
    fact =k*fact
end

endfunction
Run ใน Scinote แล้วพิมพ์
Jinna(5)
-----------------------------
การใช้ฟังก์ชั่น IF ELSE
function s=signum(x)
if(x>0)
    s='positive'
elseif(x<0)
    s='negative'
else
    s=0
end
endfunction
Run ใน Scinote แล้วพิมพ์
signum(-123)
signum(123)
signum(0)
-------------------------------------------
clear
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัว=%.11f\n',Vt)
อันนี้คือเขียนสมการ

-----------------------------------------------------------------------------
x_vais=[.5 2 2.5 3 3.5 4 4.5 5] y_vais=[0 2.7 4 5 6 7 8 9] plot2d(x_vais,y_vais,leg='Line of Best Fit',style=10) plot2d(x_vais,y_vais,style=-5) xtitle('Data.Analysis')
------------------------------------
x_vais=[.5 2 2.5 3 3.5 4 4.5 5] y_vais=[0 2.7 4 5 6 7 8 9] plot2d(x_vais,y_vais,leg='Line of Best Fit',style=10') plot2d(x_vais,y_vais,style=-5) xtitle('Data Analysis') a=[.5 2 3 4 5 6] b=[0 2.5 3.5 4.5 5.5 6.5] plot2d(a,b,style=5) plot2d(a,b,style=-9) xtitle('Data Analysis')
-------------------------------------
n=[3 5 9 17 33 65];
s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=-1)
plot2d(n,s,style=5)
----------------------------------------------------------------------------------------------
clear
a=input("แร่ป้อน=","string");
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัวของ %s เท่ากับ%.11f\n',a,Vt)

การให้มันประกาศคำ

------------------------------------------------------------------------------------------
clear

a=input("แร่ป้อน=","string");

d=[1:10:1000]; ทำให้มันเป็นค่า 1 ถึง1000 โดย ห่างช่วงละ 10

u=input("ความหนืดของตัวกลาง(kg/m.s.)=");

g=input("แรงโน้มถ่วง(m/s^2)=");

s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");

s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");

S=(s1-s2);

Vt=(1/18)*(S/u)*g*((d*10^-6)^2);
xtitle("การตกตัวของอนุภาค"); ใส่ชื่อกราฟ
plot2d(d,Vt,style=9);
xgrid(10);

 ------------------------------------------------------------------------
D=input ("Blast Hole diameter in mm.=");
BH=D/15
printf('\nBench height >= %.2f metres\n',BH);
B=30*D
printf('\nBurden = %.2f\n',B);
S=1.5*B
printf('\nSpacing = %.2f\n',S);
Sd3=3*D
Sd15=15*D
printf('\nSubdrill Range= %.2f to %.2f\n',Sd3,Sd15);
SD=input ("Please select the Subdrill you want to use = ");
C=20*D
printf('\nRecommend Charge length >= %.2f\n',C);
SL=20*D*B
printf('\nStemming >= %.2f\n',SL);
BSR=BH/B
printf('\nBurden stiffness ratio = %.5f\n',BSR);
if(BSR>2&BSR<3.5)
BSR='good'
else BSR='very good'
end
SMS10=D/10
SMS20=D/20
printf('\nStemming material size= %.2f to %.2f\n',SMS10,SMS20);
L=BH+SD
printf('\nHole length = %.2f\n',L);
C2=L-SL
printf('\nCharge length = %.2f\n',C2);
N=input ("Number of Holes in a blast=");
DS=input ("Density of rock in t/m^3 =");
V=B*S*BH*N
printf('\nBlast Volume = %.2f\n',V);
T=V*DS
printf('\nBlasted tonnes = %.2f\n',T);
Vb=3.14*((D^2)/4000)*L
printf('\nVolume of blasthole = %.2f\n',Vb);
Ex=input ("Explosive density=");
Kg=Vb*Ex
printf('\nMass of explosive per hole = %.2f Kg\n',Kg);
PF=(Kg*N)/T
printf('\nPower Factors = %.2f Kg/m^3\n',PF);
-----------------------------------------------------------------------
ทำทิ้งไว้
clear;
clc
printf('Do you want shape 2D or 3D?');
T=input("If you want 2D Print 20 , 3D Print 30 =");


if T==20 then
printf('You want to find area of "Circle" or "Annulus"?');
T=input("If you want Circle Print 21 , Annulus Print 22 =");
if(T==21);
r=input("Radius =");
d=2*r
printf('\nDiameter =%.2f \n',d);
A=3.14*(r^2);
C=2*3.14*r;
printf('\nArea =%.2f \n',A);
printf('\nCircumference =%.2f \n',C);


elseif(T==22);
r=input("Inner radius =");
R=input("Outer radius =");
P=0.5*(r+R);
W=R-r;
A=3.14*((R^2)-(r^2));
printf('\nAverage Radius =%.2f \n',P);
printf('\nWidth =%.2f \n',W);
printf('\nArea =%.2f \n',A);

3D
if(T==30);
printf('You want to find area of "General cone or Pyramid" or "Frustum of cone"?');
T=input("If you want General cone or Pyramid Print 31 , Frustum of cone Print 32 =");
if(T==31);
A=input("Area of base =");
h=input("Height =");
V=(1/3)*A*h;
printf('\nVolume =%.2f \n',V);
if(T==32);
r=input("Top Radius =");
R=input("Base radius =");
h=input("Height =");
s=input("Slant height =");
V=(3.14/3)*((r^2)+(R*r)+(R^2))*h
S=3.14*s*(R+r)+(3.14*(r^2))+(3.14*(R^2))
printf('\nVolume =%.2f \n',V);
printf('\nSurface Area =%.2f \n',S);

!= คือเครื่องหมายไม่เท่ากับ
----------------------------------------------------------
clc
x=1;
while x==1

D=input("How many Dimentions =?");
       if  D==2 then
                 ch1=input("Sector of circle press 1\nPararellogran press 2   ");
             if ch1==1 then
                    r=input("radius= ?");
                     a=input("Angle in ra =?");
                     Ar=0.5*a*r*r;
                     s=a*r;
                     printf("Area = %.2f",Ar);
                     printf("\n Arc Length = %.2f",s);
               elseif ch1==2 then
                      b=input("base = ?");
                      h=input("height = ?");
                      a1=input("Side = ")
                     Ar1=b*h;
                     P= (2*a1)+(2*b);
                     printf("Area = %.2f",Ar1);
                      printf("\n Perimeter = %.2f ",p);
                    else
                    disp("Error404");     
            
             end
             
       
     elseif D==3  then
            ch1=input("Right circular cone press 1\nQuare Pyramid press 2   ");
             if ch1==1 then
    
                    r2=input("radius= ?");
                     h2=input("height = ?");
                    Vo1=0.33*3.14*r2*r2*h2;
                    s=(3.14*((r2^2)+(h2^2))^0.5)+(3.14*r2*r2);
                    printf("Volume = %.2f m3",Vo1);
                    printf("\n Surface Area  = %.2f  ",s);
             elseif ch1==2 then
                     s=input("side = ?");
                    h=input("height = ?");
                    vo=0.33*(s*s)*h;
                     sf=s*(s+((s^2)+(4*h*h))^0.5);
            
            printf("Vol = %.2f   \n Surface Area = %.2f",vo,sf);
            
            end
            
end
x=input("\n \n \n **** Repeat press 1   ****");
end












Ex 19


Ex 18

จุดตัดแกน X (-3.426, -5.660x10^-5)


จุดตัดแกน Y (5.000, -7.231x10^-6)

Ex 17


จุดตัดแกน X (-5.043, -0.043)


จุดตัดแกน Y (0.418, 5.042)