% HBaa - 12.06.2015
%
%
clear all

A = [sin(pi/4)   16*sin(pi/2)   81*sin(3*pi/4);
     sin(pi/2)   16*sin(pi)     81*sin(3*pi/2);
     sin(3*pi/4) 16*sin(3*pi/2) 81*sin(9*pi/4) ];

% WOLFRAM ALPHA // analog //
Ainv=[648*sqrt(2) 1296 648*sqrt(2); 81 0 -81; 8*sqrt(2) -10 8*sqrt(2)]/2592

k = Ainv*[1 2 3]' % * q0/(4*EI_y) * l^4/pi^4

% Balkenachse: l = 1000 mm
x=0:10:1000;

% Näherunglsg.
w_sch = ( k(1)*sin(pi*x./1000) +k(2)*sin(2*pi*x./1000) +k(3)*sin(3*pi*x./1000) )/4/pi^4;



% exakte Lsg. aus TM2-Buch, S. 115 (4)
w_ex = 1*(7*x./1000-10*(x/1000).^3+3*(x/1000).^5)/360;


plot(x,-w_sch,'r.',x,-w_ex,'b-')

%
% *** ENDE ***