clear all;clc;close all
fishery_main
close all

Jmin=cell2mat(dyn.Jo);
Uopt=cell2mat(dyn.Uo);
x=cell2mat(res.X);
u=cell2mat(res.U);
tiempo=linspace(0,200,length(x));
tiempo2=linspace(0,200,length(u));
plot(tiempo,x,'b');
title('Peces en el estanque');
xlabel('t [días]')
ylabel('Peces ')
grid on
figure;plot(tiempo2,u,'b');title('Tasa de pesca');
xlabel('t [días]')
ylabel('u [peces/día]')
grid on


figure
mesh(linspace(0,200,length(Jmin(1,:))),linspace(0,1000,length(Jmin(:,1))),Jmin)
view([0 90]);grid('on');colorbar;
xlabel('días');ylabel('peces'); title('costo mínimo')

figure
mesh(linspace(0,200,length(Uopt(1,:))),linspace(0,1000,length(Uopt(:,1))),Uopt)
view([0 90]);grid('on');colorbar;
xlabel('días');ylabel('peces'); title('control óptimo')