clc; clear; close; t=320; p=1; while p