cantilever_beam_deflection_euler.sci
clear; clc; P = 1000.0 // N L = 2.0 //m E = 200e9 //GPa b = 0.1 //m h = 0.2 //m I = b*h^3/12.0 function dydx = f(x,y) M = P*(L-x); dydx(1,1) = y(2); dydx(2,1) = -M/(E*I); endfunction x = 0:h:L; n = length(x); y = zeros(2,n); y(1,1) = 0.0; // warunki poczatkowe dlaczego 1, 2 y(2,1) = 0.0; for i=1:n-1 k1 = f(x(i), y(:,i)); y(:,i+1) = y(:,i) + h*k1; end plot(x,y(1,:),'r') plot(x,y(2,:),'b'))