public define plot_three()
{
   fancy_plot_unit("kev","ergs");
   multiplot([5,2,2,2]);
%%%
   mpane(1);
%%%

   variable scl1_0 = get_par("constant(1).factor");
   variable scl1_3 = get_par("constant(3).factor");
   variable scl2_0 = get_par("constant(5).factor");
   variable scl2_3 = get_par("constant(6).factor");
   variable scl3_0 = get_par("constant(8).factor");
   variable scl3_3 = get_par("constant(9).factor");

   notice_values([[1:9],[11:19],[21:23]],0.5,8;unit="kev");
   () = eval_counts;
   plot_unfold({[[1:3],[7:9]],[4:6]};dcol={15,5},decol={15,5},
               yrange={3.301e-11,7.999e-10},mcol={0,0},
               scale={[1/scl1_0,1/scl1_3,1]},no_reset=1);
   plot_unfold({[[11:13],[17:19]],[14:16]};dcol={15,5},decol={15,5},
               yrange={3.301e-11,7.999e-10},mcol={0,0},oplt=1,
               scale={[1/scl2_0,1/scl2_3,1]},no_reset=1);
   plot_unfold({[21,23],22};dcol={15,5},decol={15,5},
               yrange={3.301e-11,7.999e-10},mcol={0,0},oplt=1,
               scale={[1/scl3_0,1/scl3_3,1]},no_reset=1);

   notice_values([[1:9],[11:19],[21:23]],0.7,1.5,2.6,8;unit="kev");

   () = eval_counts;
   plot_unfold({[[1:3],[7:9]],[4:6]};dcol={1,4},decol={15,5},oplt=1,
               scale={[1/scl1_0,1/scl1_3,1]},no_reset=1);
   plot_unfold({[[11:13],[17:19]],[14:16]};dcol={1,4},decol={15,5},oplt=1,
               scale={[1/scl2_0,1/scl2_3,1]},no_reset=1);
   plot_unfold({[21,23],22};dcol={1,4},decol={15,5},oplt=1,
               scale={[1/scl3_0,1/scl3_3,1]},no_reset=1);

%%%
   mpane(2);
%%%
   notice_values([[1:9],[11:19]],0.5,8;unit="kev");
   () = eval_counts;
   plot_residuals({[[1:3],[7:9]],[4:6]};res=6,rcol={15,5},recol={15,5},
                   yrange={0.85001,1.14999},no_reset=1);

   notice_values([[1:9],[11:19],[21:23]],0.7,1.5,2.6,8;unit="kev");

   () = eval_counts;
   plot_residuals({[[1:3],[7:9]],[4:6]};res=6,rcol={1,4},recol={15,5},
                  yrange={0.85001,1.14999},oplt=1,no_reset=1);
%%%
   mpane(3);
%%%
   notice_values([[1:9],[11:19],[21:23]],0.5,8;unit="kev");
   () = eval_counts;
   plot_residuals({[[11:3],[17:19]],[14:16]};res=6,rcol={15,5},recol={15,5},
                   yrange={0.85001,1.14999},no_reset=1);

   notice_values([[1:9],[11:19],[21:23]],0.7,1.5,2.6,8;unit="kev");

   () = eval_counts;
   plot_residuals({[[11:13],[17:19]],[14:16]};res=6,rcol={1,4},recol={15,5},
                  yrange={0.85001,1.14999},oplt=1,no_reset=1);
%%%
   mpane(4);
%%%
   notice_values([[1:9],[11:19], [21:23]],0.5,8;unit="kev");
   () = eval_counts;
   plot_residuals({[21,23],22};res=6,rcol={15,5},recol={15,5},
                   yrange={0.85001,1.14999},no_reset=1);

   notice_values([[1:9],[11:19],[21:23]],0.7,1.5,2.6,8;unit="kev");

   () = eval_counts;
   plot_residuals({[21,23],22};res=6,rcol={1,4},recol={15,5},
                  yrange={0.85001,1.14999},oplt=1,no_reset=0);
}