%
% Ricardo Carretero 08/Dec/2003
% Plots the bridge
%

function [] = bridgeplot(x,u,v,theta,omega,N,opt)

global e1 e2 d1 d2 L l K m
global allx allt ally allz alltheta allomega;
 
Npts=30;
w=linspace(-l,l,Npts);
[W,X]=meshgrid(w,x);

H=u'*ones(1,Npts)+sin(theta)'*w;

fs=16;

if(opt)
 figure(2)
 subplot(2,2,1)
 set(gca,'FontSize',[fs]);
 plot(allx,u);
 xlabel('x'); ylabel('u(x)');
 subplot(2,2,2)
 set(gca,'FontSize',[fs]);
 plot(allx,theta);
 xlabel('x'); ylabel('\theta(x)');
 subplot(2,2,3)
 set(gca,'FontSize',[fs]);
 plot(allx,v);
 xlabel('x'); ylabel('v(x)');
 subplot(2,2,4)
 set(gca,'FontSize',[fs]);
 plot(allx,omega);
 xlabel('x'); ylabel('\omega(x)');
end

figure(1)
set(gca,'FontSize',[fs]);
%surfl(X,W,H);shading interp;colormap(gray)
surf(X,W,H,H);shading interp;colormap(jet)
xlabel('x'); ylabel('w'); zlabel('H');
view([12 46])
