metoda_strzalow.sci
clear clc q = 1000.0 // Nm L = 2.0 //m E = 200e9 //GPa bb = 0.1 //m hh = 0.2 //m I = bb*hh^3/12.0 function dydx = f(x,y) M = -0.5*q*L*x + 0.5*q*x*x dydx(1,1) = y(2); dydx(2,1) = -M/(E*I); endfunction //pierwszy strzl s1 = 0.0 o1 = 0.0000449 //drugi strzal s2 = -0.00002 o2 = 0.0000049 theta0 = s1 theta0 = -o1/((o2-o1)/(s2-s1)) h=0.1 x = 0:h:L; n = length(x); y = zeros(2,n); y(1,1) = 0.0; y(2,1) = theta0; for i=1:n-1 k1 = f(x(i), y(:,i)); y(:,i+1) = y(:,i) + h*k1; end plot(x,y(1,:),'r') disp(y(1,$))